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I  INTRODUCTION 

The  performance  of  a  naturally  fragmenting  weapon  depends  upon 
the  number  of  fragments  produced  by  the  detonation  of  the  high 
explosive  (HE),  their  angular  and  weight  distributions,  and  their 
velocities.  Calculation  of  these  quantities,  in  place  of  experimental 
firing  and  data  collection,  is  especially  attractive  when  configura¬ 
tions  and  explosives  for  new  or  improved  weapons  are  being  considered. 
The  computational  scheme  described  here  was  constructed  during  the 
past  year,  under  U.S.  Army  support1.  There  are  two  major  related 
parts,  the  computational  fluid  dynamics  to  produce  the  detonation  wave 
and  find  the  gas  and  metal  motion,  and  a  fragment  prediction  scheme  to 
get  the  number,  weight  distribution  and  average  velocity  of  fragments 
in  each  of  the  polar  zones  around  the  weapon.  The  computer  program, 
while  rather  long,  is  in  BASIC  (Beginners  All  Purpose  Symbolic 
Instruction  Code),  takes  very  little  time  to  run  on  a  large  computer, 
and  does  not  require  an  experienced  programmer. 

The  fluid  dynamics  are  discussed  in  Sec. II.  A  conventional 
Lagrangian  scheme  with  artificial  viscosity  Is  used  for  the  interior 
gas  dynamics.  Gas  grid  points  are  made  to  slide  along  the  metal 
boundary  by  essentially  the  method  in  the  HEMP  code2,  but  the  slide 
point  acceleration  formulas  are  different.  The  standard  subdivision 
of  the  metal  casing  into  cells  for  which  detailed  calculations  are 
made  is  not  used.  Instead,  the  metal  is  taken  to  be  a  set  of  mass 
points  whose  motion  is  a  boundary  condition  which  is  solved  for 
along  with  the  gas  flow.  This  avoids  the  inevitable  difficulties  which 

1  Fragment  Prediction  Method,  Work  Unit  No.  NOL-989/A,  for  U.S. 

Army  Materiel  Systems  Analysis  Agency,  Aberdeen  Proving  Ground. 

2  Mark  L.  Wilkins,  '‘Calculation  of  Elastic-Plastic  Flow,"  in  Methods 

of  Computational  Physics,  Vol.  3,  edited  by  B.  Alder,  S.  Fernbach,  and 
M.  Eotenberg,  Academic  Press,  New  York,  1961!.  - 
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arise  in  the  usual  Lagrangian  scheme  as  the  metal  becomes  thinner 
and  smaller  computational  time  steps  are  called  for.  The  mass  point 
idea  is  well  suited  to  problems  where  the  metal  casing  is  relatively 
thick  and  expansions  to  2  or  3  initial  radii  are  calculated.  Since 
there  is  a  large  saving  in  computing  time  if  the  computational  grid 
is  coarse,  some  attention  was  given  to  the  establishment  of  a 
reasonable  detonation  wave  in  a  coarse  grid,  and  to  the  effect  of 
grid  size  on  accuracy. 

The  calculated  results  are  to  be  used  for  lethal  area  studies 
which  are  generally  made  with  data  obtained  experimentally  in  arena 
tests3.  These  data  consist  of  numbers,  weight  distributions,  and 
average  velocities  of  fragments  in  5  degree  polar  zones  surrounding 
the  weapon.  Section  III  contains  the  scheme  for  producing  these 
data  analytically,  by  adding  the  calculated  contributions  of  all 
the  mass  points  to  each  polar  zone.  Fragment  weight  distributions 
assigned  at  the  outset  to  each  of  the  metal  mass  points  are  used 
here  in  conjunction  with  the  calculated  mass  point  motions.  Results 
of  computations  for  the  105  mm.  Ml,  and  155  mm,  M107,  projectiles, 
filled  with  military  grade  Composition  B  explosive  are  given  in 
Sec. IV.  There  is  also  provision  for  providing  these  results  in 
other  formats,  for  example,  as  a  punched  card  deck  which  can  be 
used  directly  in  the  JMEM  lethal  area  program. 

The  complete  BASIC  computer  program  appears  in  the  appendix, 
together  with  notes,  a  list  of  variables,  and  lists  of  input  and 
output  statements. 

The  scheme  can  be  used  for  other  fragmenting  systems,  e.g., 
bombs  and  warheads.  However,  in  many  cases  fragment  weight  dlstri- 
buttons  and  detonation  product  equation  of  state  data  needed  to 
make  the  calculation  will  not  be  available.  These  data  will  have 
to  be  assembled,  either  from  existing  test  data  for  various  types 
of  casings  and  explosives,  or  from  new  theoretical  and  experimental 
work.  The  program  can  be  used,  in  conjunction  with  the  arena  test 


3  Joint  Munitions  Effectiveness  Manual,  Test  Procedures  for  High 
Explosive  Munitions,  TH-61A1-3-7,  FM101-51-3,  NAVAIR00-130-ASR-2-1, 
FMFM5-2L,  12  Jun  1970. 
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data,  to  find  the  effects  cf  various  factors  on  fragmentation,  since  it 
provides  a  way  to  tell  what  parts  of  the  casing  the  fragments 
came  from, and  the  related  detonation  wave  impact  angle  and 
acceleration  history. 

A  computer  program  to  make  this  kind  of  calculation  was 
constructed  by  Lir.demann4,  who  used  simple  approximate  formulas 
for  the  casing  motion,  in  place  of  the  computational  fluid  dynamics. 
Detailed  hydrodynamic  calculations  have  also  been  made  wich  modern 
Lagrangian  and  Eulerian  programs  which  take  into  account  elastic- 
plastic  flow  in  the  metal2»s.  The  program  constructed  here 
(Appendix  A)  is  a  compromise.  It  has  enough  detail  to  produce 
useful  input  data  for  lethal  area  calculations.  At  the  same  time 
it  is  simple  and  fast  enough  to  be  operated,  on  a  routine  basis, 
by  weapon  designers  and  test  personnel  with  no  special  Interest  in 
computer  programming. 


4  Michael  J.  Lindemann,  "A  Computational  Method  for  Predicting 
from  Design  Parameters  the  Effective  Lethal  Area  of  Naturally 
Fragmenting  Weapons,”  Naval  Ordnance  Station,  Indian  Head, 
Maryland,  IHTR  295,  30  Jun  1969  • 

s  L.J.  Hageman  and  J.M.  Walsh,  ”HELP,  A  Multi-Material  Eulerian 
Program  For  Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two 
Space  Dimensions  and  Time,"  Ballistic  Research  Laboratories, 
Aberdeen  Proving  Ground,  Maryland,  BRL-CR39,  May  1971,  (VOL  I- 
726459;  VOL  11-726460). 
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II  FLUID  DYNAMICS 
A.  Input  and  Initialization 

The  program  is  designed  to  deal  with  axisymmetric  systems  of  the 
type  shown  in  Fig.l,  i.e.,  relatively  heavy  casings,  with  closed 
ends,  filled  with  HE.  In  v?ew  of  the  axial  symmetry  we  will  be 
concerned  with  half  the  cross  section.  In  Fig.l,  the  interior  has 
been  divided  into  cells,  while  the  masses  of  the  bands  of  metal 
indicated  by  the  broken  lines  have  been  assigned  to  mass  points 
located  on  the  interior  metal  boundary. 

Input  for  the  computation  consists  of  values  of  K1  and  LI  ,  the 
number  of  computational  grid  points  in  the  axial (Z)  and  radial  (R) 
directions,  dimensions  of  the  metal  casing,  material  densities,  the 
detonation  velocity  D  »  various  other  constants,  equation  of  state  data 
for  the  gaseous  detonation  products,  and  for  each  mass  point  a  frag¬ 
ment  weight  distribution  (numbers  of  fragments  per  gram  in  various 
weight  ranges),  or .preferably, a  parameter  value  from  which  this 
weight  distribution  can  be  calculated. 

Let  Z  and  R  be  the  fixed  grid  coordinates  and  U  and  V  the  velocity 
components  in  the  axial  and  radial  directions,  respectively.  Denote 
the  time  byT,  the  pressure  by  PI  ,  the  artificial  viscosity  byQl, 
the  density  by  p ,  the  density  of  the  undetonated  HE  byp0,  the  relative 
specific  volume  by Vl(-p0/p) ,  and  the  internal  energy  times  p0  by E2  • 

Also  let  2irWbe  the  mass  associated  with  a  grid  point,  and  ZirWlthe 
mass  associated  with  a  cell,  where  A  andAl  are  the  corresponding 
areas.  The  flow  variables  Z,R,U,V  W,A  will  be  located  at  the  grid 
points,  while  pi,Ql#vi,E2,  R1,21,W1  ,  will  be  located  at  the  cell 
centers.  The  quantities W  and  W1  (mass/2lT)  will  be  called  scaled 
masses. 

The  construction  of  the  computation  grid  proceeds  as  follows 
(see  Fig. 2):  Let  the  outer  metal  boundary  consist  of  the  grid  points 
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(O.L)and  (Kl  +  l.L)  for  L=0  to  Li,  i  outer  metal 

(K,L1H)  for  K*0  to  Kl+i..  )  boundary. 

These  points  will  be  used  only  in  the  calculation  of  the  scaled 
masses  WK>t  to  be  associated  with  the  metal  mass  points.  The  inner 
surface  of  the  metal  casing  will  contain  the  boundary  grid  points 
used  in  the  computation,  namely, 

(l.L)and  (Kl,L)  for  L-O  to  Ll-1 ,  i  inner  metal 

(K,L1)  for  K  =  i  to  Kl.  I  boundary. 

Then,  1 

a)  Insert  values  of  20t0  ,  Hij0  ,  HKl(0  .  HK1+1>0  . 

b)  Calculate  uniformly  spaced  interior  axis  points 

2„*  K.  *  «»..  -?«)  ^-DAKI-D- 

c)  Initialize  the  remaining  Z  by  putting 

2„>l=  Hh  #  for  K  =  0  to  Kl+i,  L-  O  to  Ll+1. 

d)  Read  in  RK1A  and  RK>LA+t  for  K*0  to  Kl  +  1. 

e)  Initialize  uniformly  spaced  R.  values  by  putting 

Rk,u=  RK,.t‘  L/Ll ,  for  K-O  to  «l+l.L*Oto  Ll-1 

It  is  assumed  above  that  the  end  walls  of  the  weapon  are  perpen¬ 
dicular  to  the  axis.  If  these  walls  are  curved,  special  grid  point 
values  must  be  read  in.  An  example  is  the  base  of  the  105  min 
projectile,  for  which  the  grid  is  shown  in  Fig. 3. 

Coordinates  of  the  interior  cell  centers, Z1K>U  andRlKtLare  gotten 
by  averaging  coordinates  of  the  four  cell  corners. 

In  order  to  prevent  large  distortions  of  the  mesh,  the  grid  lines 
extending  upward  from  the  line  L=L1~1  are  allowed  to  slide  along  the 
line  L*L1.  The  intersections  with  the  lineL=Ll,  are  denoted  by 
H3k,R3w.  The  notation  ZK  lA ,  R*,^  is  used  for  the  metal  mass  points. 

At  the  outset  H3K«HK(L4  and  R3H»RK>1.t .  At  the  ends,H3*«Zltil%RVR*.w  . 

x  throughout  the  calculation.  This  notation  and 
the  variable  locations  are  shown  in  Fig. 

There  is  provision  for  an  inert  compressible  material,  with  its 
own  density  and  equation  of  state,  to  occupy  the  space  between  the 
lines  K=1  and  KsK4>l,  where  K4  is  an  integer  specified  as  part  of  the 
initial  conditions.  This  takes  into  account  the  fuze  mechanism,  if 
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it  protrudes  into  the  HE  cavity.  If  the  fuze  does  not  protrude  into 
the  cavity,  K4=o. 

An  initial  density  R2K>U, either  the  inert  material  density  or  the 
solid  HE  density,  is  stored  for  each  cell. 

The  grid  point  locations  are  used  to  calculate  the  values  of  the 
scaled  masses  Wl„, associated  with  the  cells.  Pig. 5  shows  how  the 
volumes  are  calculated  by  decomposition  into  triangles.  For  example, 
to  get  the  volume  associated  with  a  cell,  using  the  notation  in  Pig.1*, 
let  the  vertices  of  one  triangle  be 

Y1  ■  Rk,i  Y2  -  Rh*j)L  Y3  =  . 

Then  calculate  the  scaled  volume  (volume/2  n)  swept  out  by  rotating  the 
triangle  about  the  axis  from 

A=[(X1Y2'X2  Y1)  +(X2  Y3-X3  Y2)  +  (X3  Yt-X1Y3)]/Z^ 

V3  =  IAKY1  ♦y2+'Y3 )/3  . 

A  similar  calculation  if  made  for  the  adjoining  triangle  in  the  cell. 
Then 

W1ku  =  (y/3+V4)  R2*>l  ^  (2) 

where  V3  and  V4  are  the  two  sealed  volumes  and  R2K,  is  the  initial 
density. 

The  scaled  masses  w,  ,  for  the  interior  grid  points  and  the  metal 
nass  points  are  calculated  in  a  similar  way.  The  subdivision  of  the 
various  areas  into  triangles  is  shown  in  Fig. 5.  Both  the  scaled 
masses  of  the  metal  mass  points  and  those  associated  with  the  Interior 
grid  points  are  denoted  by  The  scaled  messes  of  the  HE  associated 

with  the  gas  grid  points  on  the  boundary  are  called  W3-...  The 
quantities  wi„  ,  .w* ,  W3„  L  are  calculated  once,  at  the  beginning,  and 
then  used  throughout  the  calculation.  The  quantities  w/3^ are  not 
used  in  calculating  the  metal  mass  point  motions  because  the  boundary 
point  pressures  are  obtained  by  extrapolation  and  Interpolation.  They 
are,  however,  used  in  the  slide  routine  and  in  the  energy  check  made 
at  the  end  of  each  computation  cycle. 

The  velocity  components,  relative  specific  volumes, and  energies 
must  be  specified  for  the  initial  time.  Usually  one  sets  u“,*o, 
v*j(«o,  l .  Also  E2^t=o  for  K=i  to K4 (if Wo), and  EZ^Elp,, where  El 
is  the  energy  released  per  gram  of  HE,  for  K  >  K4. 
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The  detonation  can  be  started  on  the  left  boundary  of  the  HE, 
on  the  line  K  =  K4+L  The  initial  time  is  taken  as 

~  (  ^K4*l.O  ■*2K^4,  e>)/P 

the  time  for  the  detonation  front  to  traverse  a  half  cell.  There  is 
also  provision  for  starting  the  detonation  at  the  point  O). 

B.  Flow  Equations 

The  equations  governing  the  flow  can  be  written 
Velocity 

d  €  (d.  T  -  U  J  dR/dT-V  (3) 

Continuity 

p  AR  =  W  ,  pi-Ai  Rl=Wl  (U) 

Artificial  viscosity 

Q1  =  C3*p.  fAl/Vl)[(l/Vl)  dvi/dr]  (5) 

Equation  of  state  of  detonation  products 

Pi  =  Pl(£2,vi)  •  F  (6) 

P2  —  ■¥  Q1  (V) 

dE2/dT  +  P2  dvi/dT  =0  (3) 

Acceleration 

Interior  points 

du/d t  -  -  (jL/f>)  apst/az  y  dv/dr  -  -(i/p)  ap2/aR  ( ° ) 

Metal  mass  points 

2irwdo/dT=X,P5J  21TV/dv/dT=  A«R5  (10) 

In  (5)»C3  is  a  constant  and  a  quadratic  artificial  viscosity  has 
been  used.  Equation  (6)  contains  the  burn  fraction F,  the  fraction 
of  the  cell  traversed  by  the  detonation  front. 

Provision  is  made  for  detonation  product  equations  of  state  of 
the  form 

PI  =  C4 *  EZ  +C5.  (11) 


Here,  C4-  andcs  are  functions  of  the  relative  specific  volume  vi  . 
The  sound  speed  is  needed  in  the  time  step  calculation.  Since 


C  2  -0>Pi/dp)s  =  - ( vi2/^ )( a pi/svi)s 

OEz/avi)s  *  -  pi  , 


where  C2  is  the  sound  speed  squared  and  S  is  the  entropy,  we  have 
from  (11) 


NOLTR  74-77 

C2  =  -{\^)[dC5/dvl  -  Pl  C4  *(dC4/dvi)  El].  (12) 

Calculations  are  presently  being  made  with  the  JWL  equation  of 
state*,  which  has  the  form  (11),  with 

C*=  D^/V1  '  (13) 

C5-  D?Cl-C4/Fl)exp(-Fi/vl)  t  D4(t-C4/F2)exp(-F2Vl)j 

where  D2,  p5,  D4,  Fl  j  and  F2  are  constants. 

Let  K  and  t.  be  Lagrangian  coordinates,  or  grid  point  labels. 
Initially  m  the  z  and  r  directions,  respectively.  Then  all  the  flow 
variables,  including  *£  and  R  ,  are  functions  of  k,  L  »  andT.  Now 

3PZ  /3K  *»  (  3P2/3R)  2R/9K  +  (2P2/SZ)  2B/2K  , 

2*>2/dL  *  ( 2P2/3R)  3R/3L  «■  (0P2/9Z)  3Z/3L  , 

from  which, 

dPz/9%  =[(aP2/2K)  a?/5>L  -(ap2/aL)az/dK]  / 3(r,i)/3(*x)  , 

»P2/aZ  »[(3P2/3L)  3R/3K  -(3P2/3K)  dR/3L]/c>(R,*)/3(K,  L)  .  (1' 

Combining  ( ,  and  (9),  and  using  the  fact  that 

A  *  ~  9(R,?)  /2Ck,L)  , 

leads  to 

du/dr  *(R/W)[(aP2/,9L)3R/3K  -  (PP2/&K)  0R/3lJ  , 

{lb; 

dV/dT  =(R/w)[(aP2/5K)3f/aL  -  OP2/»L)d*/aK]  . 

The  equations  in  (16)  are  used  to  calculate  the  gas  grid  point 
accelerations.  The  motion  of  a  metal  mass  point  on  the  boundary  is 
found  from  (10),  where  2itwis  the  mass,  Aa  and  A*  are  projected  areas 
in  the  Z  and  R  directions,  and  P5  is  the  pressure  associated  with  the 
mass  point.  The  projected  areas  Aa  and  Apare  found  by  working  with 
lines  connecting  the  mass  points. 

*  Metal  Acceleration  of  Chemical  Explosives.  J.  W.  Kury,  H.  C.  Hornig, 
E.  L.  Lee,  J.  L.  McDonnel,  D.  L.  Ornellas,  M.  Finger,  F.  M.  Strange, 
and  M.  L.  Wilkins.  Proc.  Symp.  Detonation,  4th,  Office  of  Naval 
Research,  Rept.  ACR-126,  pp.  3-13*  U.  S.  Govt.  Printing  Office, 
Washington,  D.C.,  1965. 
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The  motion  of  a  slide  point  is  shown  in  Pig, 6.  Suppose  the 
slide  point  to  be  moved,  E3K,R3*  is  located  on  the  line  Joining 
the  mass  points  labeled  1  and  2  in  Fig. 6.  The  slide  point  is  first 
moved  along  the  line  Joining  the  old  grid  points  1  and  2,  to 
The  point  E4M>R4Hmay  extend  beyond  the  point  1  or  2.  This  is  done 
Just  before  the  new  grid  point  velocities  are  calculated.  At  the 
beginning  of  the  following  cycle,  after  the  interior  grid  points 
and  the  mass  points  have  been  moved  to  their  new  locations,  the 
new  location  of  the  slide  point  I3HfR3K  is  found  by  intersecting 
the  line  Joining  24«,R4Hand  with  the  new  boundary. 


C.  Main  Routine  -  Difference  Equations 

Figure  7  is  an  outline  of  a  computation  cycle.  At  the  beginning,  the 
cycle  number  n  is  advanced  by  setting  Relative  to  this  new 

cycle  number  we  have  the  previously  calculated  quantities  (see  Fig. 7) 


T2  = 

(17) 

H-l  n-1  _4* 

J  K«.L  » 

ViH'1  Qi"-**  P2n'X 

EZ*1”1  C2n~l  Al*"* 

Following  the  order  in  Fig. 7  put 

Th  a  Th_1  +  Tt  . 


Limited  Computation 

Assume  that  the  detonation  starts  on  the  surface  K»K4*ior  at  the 
point  k*K4+i,  L=o.  Let  K3  be  the  maximum  k  for  which  the  cell  center 
variables  Pl,Vi,Ql,  E2.C2  will  be  calculated  and  let 

Set  ~ 

K3  =  5  +  the  maximum  K  for  which  D-T-Z.  ,  +  Z>0  ) 

(l8) 

K3  =  min  {  K3  ,  Ki -1 }  •  I 


Here  D  is  the  detonation  velocity.  The  additional  five  k  lines  in  (18) 
are  needed  to  establish  a  detonation  wave  peak  pressure  close  to  the 
Chapman- Jouguet  (CJ)  pressure.  This  works  because  the  artificial 
viscosity  term  causes  the  solid  ahead  of  the  detonation  front  to  be 
artificially  compressed. 
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{♦>  ADVANCE  CYCLE  NUMBER  N  AND  TIME  T. 

CALCULATE  NEW  VALUES  IN  THE  ORDER  LISTED: 

1.  MAXIMUM  K  FOR  LIMITED  COMPUTATION  IF  THE  DETONATION  FRONT  HAS  NOT 
REACHED  K  =  K1. 

2.  POSITIONS  OF  INTERIOR  GRID  POINTS,  AXIS  POINTS,  AND  METAL  MASS 

*  1  POINTS  z£  L  Rj^  l. 

3.  POSITIONS  OF  SLIDE  POINTS  Z3jj,  R3g. 

4.  BURN  FRACTION  FK. 

5.  AREAS  AlJI  L.  RELATIVE  SPECIFIC  VOLUMES  Vlg!  L,  ARTIFICIAL 
VISCOSITIES  Q1*'1/2  •  PRESSURES  Pit!  .,  ENERGIES  E2S  .. SOUND  SPEEDS.  AND 
CELL  TIME  STEPS  T3$!  L. 

A  6.  TIME  STEPS  T1  AND  T2. 

7.  MOTION  OF  SLIDE  POINTS  TO  Z4g,  R4$!. 

8.  VELOCITIES  OF  INTERIOR  GRID  POINTS,  AXIS  POINTS.  AND  METAL  MASS 
POINTS  ugV*.  Vjj+L1/2. 

AT  SPECIFIED  TIMES,  PRINT  OUT  RESULTS  AND  MAKE  FRAGMENT  PREDICTION 
CALCULATIONS. 

MAKE  ENERGY  CHECK. 

■<  RETURN  TO  START. 


FIG.  7  FLUID  DYNAMICS-COMPUTATION  CYCLE 
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New  Positions 

The  new  positions  of  the  interior  grid  points,  axis  points,  and 
metal  mass  points  are  calculated  from  (3).  In  difference  form 


ZN  =  zH"4 


.  p 

+  T1  •  SJL 


o'*  ^R^1  +T1. 

kh,l  ^*,1.  vk,«-  •  ' 

After  all  the  ZtR  points  are  moved,  the  new  slide  points  Z3",R3* 
are  found  by  intersecting  the  lines  joining  (calculated 

during  the  previous  cycle)  and  2"  u.t ,  (see  Pig. 6)  with  the  r 


during  the  previous  cycle)  and  Z  t (see  Pig. 6)  with  the  new 
boundary.  To  do  this  the  appropriate  segment  of  the  boundary  must  be 
found  for  each  slide  point.  Suppose  the  equation  of  the  line  through 
Z < , and  ZK^.t  ,  is^ 

d  =  at  +bR  +  c  «o,  (20) 

If  the  coordinates  of  a  point  not  on  the  line  are  inserted  for  %  and  £ 
in  (18),  then  the  sign  of  d  will  depend  on  the  side  of  the  line  the 
point  is  on.  The  proper  boundary  segment  can  thus  be  found  by 
inserting  the  mass  point  coordinates  ,  successively,  into 

(18),  starting  with  J~ 1,  and  testing  for  a  sign  change.  Since 
Z4„  -Z„  .  may  be  zero,  it  is  convenient  to  put 

SS*  -  -  (Z^-2^.0  ,  (21) 

<*  =  +  S8K(R-RK,LW)  ,  (22) 

and  substitute,  successively  Z,iLi  and  R»  ,  for  Z  and  R .  After  the 
value  of  J  for  the  left  endpoint  of  the  proper  line  segment  is  found, 
Z3k and  R3K  are  obtained  from  the  intersection  of  the  two  lines,  i.e., 
by  solving  the  system 

S7  -  (Rj-44jtl  ”  ^  *”  )  >  (23) 

(R-  R»i(u)  -  S7  (Z-2TntU)=0  ,  (24) 

<Z-  z*.,.-*)  +  SeH(S 

The  variable  is  now  set  equal  to  j,  to  identify  the  first  mass 
point  to  the  left  of  the  slide  point  Z3k,R3k. 

Burn  Fraction 

The  term  burn  fraction,  denoted  here  by  F ,  refers  to  the  fraction 
of  an  HE  cell  considered  to  be  converted  to  detonation  products.  If 
the  detonation  front  is  taken  to  be  perpendicular  to  the  axis,  F  is  a 
function  of  K  only.  Then,  with  D  the  detonation  velocity,  we  take 
(see  Ref. 2) 


15 


7 


t:.  **■  5f#V2“?  4 


aarwa^Si 


NOLTR  7^-77 


F*  =  max  {(D«T-ZKo4-?)/(ZK^tl-ZK>#)  ,  (l-VlK  J/(l  -  Vl„)  } 

(25) 

Fk  *  min  {  i,  P.J  • 

Here  H  -  Z„^  at  T=o,  the  initiation  plane,  and  Vl  is  the  relative 
specific  volume  at  the  CJ  state,  a  constant  for  the  explosive. 

To  initiate  at  a  point  let  the  burn  fraction  be  the  ratio  of  the 
time  the  detonation  front  has  been  in  the  cell  to  the  time  for  the 
detonation  front  to  traverse  the  cell.  Let  be  the  coordinates 
of  the  initiation  point  at  T=o  and  let 

G«.u=[(rk  .u-R)*  +(Zk.,-2)2.T\  (26) 

Then  take 

FK>t=  max  -  GK>J  ,  (1"  ~  vlc»)  ,  ^ 

Fv.u  *  min  {  >  FU|t.T|  . 

At  the  present  time  there  is  provision  in  the  program  for  plane 
detonation  only  (F*F\,).  Initiation  at  one  or  more  points  can 
be  done  by  making  appropriate  changes  in  the  burn  fraction  and 
equation  of  state  routines. 

Area,  Specific  Volume,  Pressure,  Energy,  and  Sound  Speed 

The  area,  relative  specific  volume,  pressure,  energy,  sound  speed 
and  cell  time  step  are  all  calculated  for  a  single  cell  before  moving 
on  to  the  next  cell.  To  get  the  area  Al„.  the  cell  is  divided  into 
two  triangles  and  the  triangle  routine  (1)  is  used.  Note  that  for  the 
uppermost  cells  (see  Fig. 4)  the  corners  are 

(ZS^.RS^),  (Z3k,R3K).  Let  V3  and  V4  be  the 
scaled  volumes  of  the  two  triangles,  obtained  with  (1).  Then 

V1K,L  =  (V3  +V4)-  R2K,L/wiK  t  ,  (2S) 

The  artificial  viscosity  (see  (5))  is  now  calculated  from 

<<r = c  ■  r  ^  -  c  )/c <  -  > ,  ( 29 ) 


where 


VI  -  r c  Vl  1  -  <* )  /T1J  ■  z/Cv C  *  vC). 


If  v<\  >o,then  =  o. 
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Arrays  in  storage  are  for  a  single  time  step  only.  The  values  of 
Al"**  .Vi;;  (  pi"'t*  ,  and  E2**^  are  stored,  temporarily,  as  unsubscripted 
quantities,  before  the  new  values  are  found,  for  use  in  (29),  (30), 
and  in  the  pressure,  energy  and  cell  time  step  calculations. 

The  pressure  Pi  and  the  energy  E2  are  found  by  solving  (6)-(8) 
simultaneously.  When  the  JWL  equation  is  used,  C4  and  C5  are 
calculated  from  (13),  using  .  The  following  are  then  solved 
by  iteration,  starting  with  the  first  approximation 

pC=  +C5,  (31) 


<*  N-l  r.  N  M-i  \  ,  H- 1  /» 1  /  N  H-t  . 

E2kA=E2Kjt  -[(nM  +PI.J/2  +Q1^  J(V1„-  V1K>1  ). 


H  H  H 

The  new  values  of  PI  ,  ,  Vl„tL  with  derivatives  of  C+  and 

C5  from  (13),  are  inserted  into  (12)  to  get  the  square  of  the  sound 
speed 
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Time  Step 

The  time  step  Is  determined  by  numfrical  stability  criteria,  with 
the  form  taken  directly  from  the  HEMP  iode2.  This  is  a  composite  of 
two  criteria,  one  for  the  shock  regions,  the  other  requiring  that  a 
signal  pass  only  part  way  through  a  cell  in  one  time  step.  The  minimum 
of  the  time  steps  found  for  the  individual  ceils  is  taken  as  the  new 


time  step. 

To  find  the  cell  time  step  we  first  get  the  smallest  diagonal  from 

I  (33) 

SS^ijnin  {  S3,  S+}]  #  (34) 

For  the  uppermost  cells,  where  L=Li-t,  23K  and  R3k  are  used  in  place 


*^K,ct  an<* 


Then  the  cell  time  step  T3  is  given  by 


ss/(3«t.  e*D  , 

where  8  =  0  lfvs»o,  (35) 

0=2'  C5  SS-V<1  if  vqiO 
The  quantities  C3  and  vq(=vi/vi)  appear  in  (29). 

The  new  time  step  is  then  found  from 

Tl"**  min  {  min  T3K(L,  Tl*  E8 }  (36) 

where  E©  is  a  constant  (1.1  is  now  being  used,  i.e.,  the  time  step  is 
not  allowed  to  jump  more  than  10?  in  a  single  time  step).  Also  set 

T2N*1/i*  _  ^  Tl***  +•  Tl")/2  .  (37) 

The  time  step  T1  will  be  used  in  the  next  computation  cycle,  but 
T2H+t/*  will  be  used  in  this  cycle,  to  get  the  velocities. 

Slide  Point  Motion 

To  move  the  slide  point  along  the  line  joining  points  1  and  2  in 
Fig. 6,  let  be  the  velocity  along  the  line.  Then 


dG^/dT  *  -(i/p)[9P2/c>2  cos  oi  +  9P2/3R.  sin  d] 


(38) 
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du7/dr  -  dG<\  /dT  cosot  j 

(39) 

dvydr  -  dG^/dT  sino<  ) 

(40) 

dE/dT*U7  ,  dR/dT  =  V7, 

where  d  is  the  angle  with  the  E  axis  shown  in  Fig. 6  and  u7,V7are  the 
velocity  components.  To  find  3P2/3Z  and  9P2/3Rin  (38),  we  can  use 
(15).  Since  the  line  L=Li  is  being  held  fixed,  3P2/3L -O.  Also,  since 
there  are  only  two  quarter  cell»  associated  with  the  point 
equation  (4)  must  be  replaced  by 

pAR  =  2  W3  .  (i»i) 

Then,  using  (15),  (38),  and  (4l), 

CdG^/d-r)^  =[R3K/(2  WJK)ll)]  (  3P2/aK)(|  W[(SZ/aL)H  tt  sin  cL  -(3R^L)K  tl  cos  cl]  ( 


where 


(.sn/3 k)„„  =  -  P2l.,„., 

(43) 

(  31/31).,.,  =  H3."  - 

(44) 

(3#/a).,l,  -  rs;  - 

(45) 

cos  <*  *  (V  2*)/  [(2*-  2,)‘  +(R*“  Rif  ]  , 

(46) 

sin  d  *  (^-«J/r(2z-2t)4+(R*'Rx)*]1/l 

(47) 

The  subscripts  1  and  2  in  (46)  and  (47)  refer  to  the  labeled  mass 
point  locations  in  Fig. 6.  We  then  get  V7K ,  Z*„and  using  (39) 


and  (4C)  with 


m-i/x  N-  t/i 

U  7  -  U 

'x  tc,  V.l-1 


V7. 


24- 


*♦1 


-  ?3« 


+  (dca/dT;swCos  ot-T2 


N-4/Z 

*  V- 


fCdG^/dT) 
H+Vx 

+  U7k 
+-V7. 


N*i/X 


KjlA  sin  «<  T2 

•  Tl 

•  T1 


(48) 


(49) 


Note  that  here,  as  in  the  HEMP  code  ,  the  old  velocities  at  the  point 
(k,  ul-i)  are  used  to  calculate  cnK  and  viK  from  (39)  and  (40). 
Velocity-Interior  Grid  Points 

New  velocities  of  the  interior  grid  points  not  on  the  axis  are 
found  from  a  difference  form  of  (16),  viz„ 
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<r  *  ^  +  t2  cm/d«)K>jd«/*i-\L] t  (so 

<;l/l  ’  C"  +  T2  (<. /%,  J  [^P^K)K  l(&?/9L)KA-(aP2/^^?/dk)K  t  ]  ^  (51) 
where 

(?P2/*K)Kl  =  (  P2"K(t  4-  P2^t.,  -  P2l.t^  -  P2;.1>t.O  /2  ,  (52) 

CSP2/3L)^  =(P2:.t>  4  P<c  -  P^K-t. c-1  -  ^-,.0/2  ,  (53) 

(  3R/(9K)Kjt  -  (  Rn4<,t  ~  > 

(SP/^L)^  =  (<,»*.  ~  <..-0/2  , 

(54) 

3  (Ci..  -2L.0/2  , 

(92/^L)K,t  =  <<*>,-  Hh.c-0/2  . 

Velocity-Interior  Axis  Points 

For  the  Interior  axis  points  (9)  will  be  used  directly.  Note 
first  that  on  the  axis  dv/dT«0,so  that  £P2/dR«0,  and  from  (14) 

SP2/5Z  - OP2/c>L)/(df/aL)  (axis  points).  (55) 

The  density p  needed  in  (9)  is  given  by 

p  *  /bi  ,  (56) 

where  81  is  the  sum  of  the  scaled  volumes  of  the  two  quarter  cells  at 
the  point  (k,o)  .  To  get  these  volumes  the  centers  2iHO  ,  Ri*>0  of  the 
adjoining  cells  are  found,  first,  by  averaging  the  H  and  R  coordinates 
of  the  corners.  Then  the  scaled  volumes  of  the  two  quarter  cells  are 
found  by  subdivision  into  triangles  (see  Fig. 8)  with  the  same  quarter 
cell  routines  that  are  used  in  the  initialization.  The  velocity 
components  are  now 

O  *  <7*  -  T2  (81  /vO  2  ( PC  -  P2".-^)/CZ.„,. -H^)  ,  ( 57  ) 

v„7«o.  (58) 

Veloclty-Hetal  Mass  Points 

The  equations  in  (10)  are  used  to  calculate  the  velocities  of  the 
metal  mas3  points.  The  notation  is  shown  in  Fig. 8.  The  pressures  at 
the  slide  points,  denoted  by  PS,  ,  are  found  first,  by  extrapolation  with 
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P4=(P2,A1.1  +  P2^>tt.1)/2 

P5  =(P2,>tl.2  +  (59) 

P3y  =  (3P4-P3)/2 

The  subscripted  variable  P5r  is  different  fromPS,  a  distinction 
allowed  in  BASIC  programming. 

The  pressure  at  a  mass  point  is  found  by  locating  the  nearest 
slide  points  on  both  sides  and  interpolating  with  respect  to  distance 
from  the  slide  points.  The  adjacent  slide  points  are  located  by 
essentially  the  same  method  used  to  find  2 3 *  and  R3K  (Eq.(20)).  The 
equation  for  the  line  through  a  slide  point  can  be  written  (Eq.(21)) 

+  S8,^R-ftT<ul^)=0.  (60) 

Here  S8, ,  the  negative  reciprocal  of  the  slope  of  the  slide  line,  was 
previously  calculated  with  (21)  and  stored.  Now, starting  with  j*!, 
evaluate 

^*7  —  C Zn,la."’  Z7,i.4.-t)  *  S8,  (61) 

for  successive  values  of  J,  until  the  sign  of  G,+l  is  different  from 

the  sign  of  G,  .  Then  for  the  values  of  J  and  j+1  where  the  signs  differ, 

compute  ,  ,  ,,*Ai 

+(R..„-R5,)*]  *  {62; 

08=r(*«.u-*3r.,)’ 

and  P5  the  pressure  at  the  mass  point  by 

P5=  P3y +  07  (P3,*,  -P3,)/(D7  +  D8).  (63) 

Assume  that  the  mass  associated  with  the  mass  point 
is  uniformly  distributed  over  a  band  whose  cross  section  consists  of 
the  pieces  between  the  point  and  the  midpoints  of  the  line  segments 
joining  the  point  and  the  adjacent  points  (Fig. 8).  Then,  in  (10), 
take  the  projected  areas 

=  ( iT/4)[(RK.itW  +  R*t„)1  -  ( ]  ;  j 

" ^^)[(3  RhiU  +  +  (3  R,^,  +  Zk,m)  .  I 

Now,  using  (10), 

lC*  sU«.'u?  +T2  A,  P5/(8W#(J  , 

(65) 

v:r  *  cr  +  t2  ^ps/cb^)  , 
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where  _  _ 

A,*  (”4// r)A,  >  A„  *  <4/w)A.  . 


Similar  formulas,  with  appropriate  subscripts, are  used  to  get  the 
velocities  of  the  metal  mass  points  located  on  the  lines  k*  l  and  K*wi. 
Energy  Check 

The  total  energy  in  the  system  is  summed  at  the  end  of  every  cycle 
Denote  the  total  internal  energy  by  If.  and  the  total  kinetic  energy 
by  K£  .  The  total  energy  trT  is  found  from 


where 


E,*1E  +  KE 


IE 


62«..  /«* 


•».w 


(66) 

(67) 


KE  = 


IT 

Kll  U*0 


+7rH  f  +wwu-^+nO} 

L«© 

KA-t 

+  TrH  W3h..JU7*+V7*)- 


I 

a 


(68) 
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FIG.  8  VELOCITY  CALCULATIONS 


0 


[ZO,  ROj 

FIG.  9  WEAPON  POSITION  AND  POLAR  ANGLE 


NOLTR  7^-77 


III  COMPUTATION  OP  ARENA  DATA 

Suppose  that  the  weapon  is  located  in  the  center  of  an  arena  of 
radius  oi  and  center  at  ?o,RO  in  the  coordinate  system  used  for  the 
fluid  dynamic  computation.  Put  the  weapon  nose  and  the  zero  polar 
angle  at  the  left,  as  in  Fig.9»  and  divide  surrounding  space  into 
5  degree  polar  zones,  measured  from  20, RO.  The  total  weight,  total 
number,  average  velocity,  and  weight  distribution  of  fragments  in 
each  of  the  polar  zones  are  needed  for  lethal  area  studies. 

The  total  weights  and  average  velocities  of  the  fragments  in  the 
various  polar  zones  are  obtained  from  the  fluid  dynamics  computation. 
Assume  that  at  a  specified  time,  say  the  time  for  the  center  of  the 
weapon  to  expand  to  twice  its  initial  radius,  the  casing  has  broken 
into  fragments  and  burst,  acceleration  of  the  fragments  has  ceased, 
and  the  directions  of  fragment  motion  are  constant.  In  the  cases  for 
which  fluid  dynamic  computations  were  made,  the  fragment  directions 
were  established  early,  and  then  changed  very  little  with  time.  The 
velocities  increased  slowly  with  increasing  radius  after  case  expan¬ 
sions  to  1.5  or  2  radii.  Hence,  the  velocities  could  vary  by  a  few 
percent  depending  upon  the  time  chosen.  The  results  of  the  fluid 
dynamics  computations  for  the  first  cycle  after  the  specified  time 
are  used. 

The  metal  mass  points  are  numbered  in  order  starting  with  1  at 
K=  2 ,  L-o  and  ending  with  K1+2UL  at  k=ki,  l=0.  Let  t^s  be  the  mass 
point  number,  and  let  ?MS .  ,  and  VM,  be  the  corresponding 

coordinates  and  velocity  components,  respectively,  calculated  in  the 
fluid  dynamics  section  of  the  program.  The  midpoints  of  the  line 
segments  joining  the  mass  points  are  then 

+  ,  rms«-v*  ~(^»r  +  KMSfi.)/z  (69) 

Also,  let  j3ns  be  the  fragment  direction  angle  for  the  mass  point  M5 

i  •  e  •  >  t  1  \ 

Pm*  “  arctan  , 

24 


(70) 


i 

! 

I 

& 

& 


m 


% 


1 


and,  for  the  fragment  directions  from  the  midpoints , take 

~(£m3  *  •  (71) 

The  corresponding  polar  angle  ,the  *:»gle  between  the  horizontal 

and  a  line  drawn  from  the  point  where  the  fragment  direction  intersects 
the  circle  to  the  point  20, Ro  (see  Fig. 9),  is  obtained  by  solving, 

simultaneously  t  2 

C?6-ZO)2  *  (R4-RO)  -  oi 

(f?4  -  «?„,*»/*)  -  tan  Pns+tft  "  ° 

for  24  and  R4,  and  taking 

arctan  [(24 -zo)/(R4- RO)]  .  (73) 

The  scaled  mass  (=wKtin  the  fluid  dynamics  program)  associated 
with  each  mass  point  M5  is  distributed  uniformly  over  the  full  and 
partial  polar  zones  between  the  polar  angles  and 

Denote  the  contribution  of  the  mass  point  nr  to  the  J2th  polar  zone 
by  W7J2M#.  After  the  scaled  masses  of  all  the  mass  points  are 
distributed  over  the  appropriate  polar  zones,  the  total  fragment  mass 
and  average  fragment  velocity  in  each  polar  zone  are  calculated  from 

KUlLl 

w07l=E  W7n,«« 

■» 


m 


1 


6 


total  fragment  mass  in  polar  zone  J2«2*r  wa,» 
average  fragment  velocity  in  polar  zone  s z 


where 


+V^f 


and 


(76) 


is  the  velocity  of  the  mass  point  ns . 

The  numbers  of  fragments  in  various  weight  ranges  in  the  individual 
polar  zones  are  obtained  from  the  directions  calculated  with  the  fluid 
dynamics  program,  together  with  weight  distributions  which  are  specified 
for  each  of  the  mass  points. 

To  every  mass  point  M5  there  is  assigned  a  value  of  the  parameter 
M ,  called  the  average  weight  of  fragments  weighing  more  than  one 

grain.  It  is  assumed  that  the  weight  distribution  of  fragments  from 
the  mass  point,  the  number  of  fragments  per  unit  weight  in  the 
different  weight  ranges,  depends  only  on  this  parameter  and  is  com¬ 
pletely  determined  once is  specified.  Let  the  total  number  of 
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weight  groups  be  J«  .  Specify  JR  weights  M*.*, . ,Ma> _ My,  »  80  that 

Mt1!  grain,  thejth  weight  group  (j  <jq)  contains  fragments  with  weights 
between  M,  and  ,  and  the  last  group  contains  all  the  fragments 
weighing  more  than  M„  .  For  each  mass  point  MS  read  in  or,  pre¬ 
ferably,  calculate  from  a  formula  the  weight  distribution  associated 
with  Mine  •  This  will  consist  of  the  quantities  NM#y  the  number  of 
fragments  of  weight  greater  thanM, ,  in  each  of  the  JR  weight  ranges, 
for  the  mass  point  Mf .  Also  read  in  or  calculate  ,  the  weight 

fractions  of  fragments  weighing  more  thanM,,  for  each  mass  point. 

Let  N7,4,  be  the  number  of  fragments  of  weight  greater  than  M, 
in  the  J4th  polar  zone.  This  number  is  accumulated  with 


MVs2irL  W  ,  (77) 

while  the  weights  from  each  mass  point  ms  are  being  put  into  the  polar 
zones  between  the  polar  angles  and  AR„g .  At  the  same 

time  the  scaled  weight  of  fragments  in  each  weight  group,  in  each 
polar  zone,  is  accumulated  with 

mill 

W5T4 .7  *  W2H»,»  W\Mr  .  (78) 

Let or  J*«W.  The  average  weight  of  fragments  in  the  j  th 
weight  group  in  thejfrth  polar  zone  is  then  zn  (ws,4,  -  ws,to 
for  j<JR,and  2ir  W5,W1/N7W  W  for  J«J* l.  The  total  number  of  fragments  of 
weight  greater  than  one  grain  in  the  Jfrth  polar  zone  is 

)  "*■  •  (79) 

The  average  weight  of  fragments  weighing  more  than  one  grain,  in  the 

J*th  polar  zone  is  2irw8,4/M8,4  (see  (7*0). 

The  value  of  R  which  must  be  assigned  to  each  mass  point  depends 
upon  the  casing  diameter  and  wall  thickness,  the  explosive  composition 
and  density,  the  type  and  treatment  of  the  steel,  and  the  impact  angle 
of  the  detonation  front.  For  the  present  work,  fragmentation  data 
from  arena  tests  of  105  and  155  mm  shells,  with  standard  projectile 
steel  casings  and  military  grade  Composition  B  explosive  fillings, 
were  used  to  construct  a  plot  (Fig. 10)  of  R  vs.  the  parameter  where 

X  « 1  dT/Cl  ♦  2  C/M)  ,  (  i"-)**  .  ( 80 ) 

Here,  t  is  the  wall  thickness  (in.),  d*  is  the  inside  diameter  (in.), 
and  C/i M  is  the  ratio  of  the  explosive  mass  to  the  metal  mass  at  the 
weapon  section  being  considered.  The  parameter  X  was  proposed  by 
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Magis7.  The  examination  of  che  arena  test  data  showed  that  when  the 
fuze  was  at  the  nose  end  of  the  projectile,  the  fragmentation  from 
the  side  wall  close  to  the  base  end,  beyond  the  pert  forming  the  main 
beam  spray,  was  substantially  finer  than  the  fragmentation  indicated 
by  the  value  of  X.  This  may  be  due  to  the  higher  pressure  in  the 
section,  produced  by  the  reflection  of  the  detonation  wave  at  the 
base  of  the  projectile.  A  second  curve  in  Pig. 10  was  drawn  for  use 
with  this  section.  Also,  in  the  arena  test  data,  the  value  of  R  for 
the  base  corner  was  less  than  one  fourth  the  largest  value  of  ft  in 
the  sections  forming  the  main  beam  spray. 

For  the  sample  calculations  that  were  made,  the  values  of  ft  for 
the  mass  points  on  the  side  wall  were  taken  from  Fig.]0,  and  theft  for 
the  mass  points  on  the  base  corner  and  the  base  end  were  set  equal  to 
those  observed  in  the  arena  test  data.  The  values  of  ft  for  the  mass 
points  on  the  line  K*iwere  arbitrarily  set  at  200  grains.  No  attempt 
was  made  to  model  the  fuze. 

Figure  10  was  constructed  with  a  limited  amount  of  work  involving 
both  available  arena  data  and  fluid  dynamic  computations.  The  arena 
data  was  used  to  get  the  weight  distributions  (number  of  fragments 
per  unit  weight  in  the  various  weight  ranges)  in  the  different  polar 
zones,  while  the  fluid  dynamics  results  were  used  to  find  the  mass 
points,  and  the  corresponding  X  values,  that  had  contributed  to  these 
polar  zones.  This  kind  of  analysis  is  felt  to  be  worth  applying  to 
other  existing  arena  data  for  a  wide  range  of  shells,  bombs,  and 
warheads.  It  has  the  advantage  that,  in  many  cases,  the  weight 
distributions  can  be  related  to  particular  sections  of  the  weapon 
where  the  value  of  the  parameter  X  (Eq.  (80)),  the  detonation  wave 
impact  angle,  and  the  acceleration  history  are  known,  from  the  fluid 
dynamics  computation. 

Some  attempts  were  made  to  correlate  the  cylinder  fragmentation 
results  of  Magis7,  and  a  weight  distribution  fitted  to  the  Magis  data*. 


7  S.F.  Magis,  Naval  Weapons  Laboratory,  Dahlgren,  Virginia. 

•  H.M.  Sternberg,  "Fragment  Weight  Distributions  from  Naturally 
Fragmenting  Cylinders  Loaded  with  Various  Explosives,"  U.S.  Naval 
Ordnance  Laboratory,  White  Oak,  Silver  Spring,  Maryland,  NOLTR  73-83* 
12  Oct  1973. 
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with  arena  test  data  for  the  standard  105  and  155  mra  projectiles.  The 
agreement,  overall,  was  poor.  A  likely  reason  for  the  poor  agreement 
is  that  the  cylinder  fragments,  collected  in  sawdust  pits,  were 
subject  to  secondary  breakup.  Another  possibility  follows  from  the 
fact  that  the  cylinder  data  are  averages  of  weight  distributions  from 
various  sections  of  relatively  short  uncapped  cylinders.  These 
averages  may  not  be  adequate  to  represent  weight  distributions  from 
sections  of  a  projectile.  The  sample  size  may  also  be  a  factor.  There 
were  very  few  standard  projectile  steel  experiments  in  the  Magis  work7. 

The  Mott  formula9  was  used  for  the  weight  distributions.  This  can 
be  written 

Nfm)  =fi/2,u)  exp  [-  (™//*) ]  ,  (81) 

where  m  is  the  fragment  size,  jUis  a  parameter,  and  N(m)  is  the  number 

of  fragments  per  unit  weight,  of  weight  greater  than  m.  With  ju in  grains, 

M(i)  is  the  number  per  grain  of  fragments  weighing  more  than  one  grain. 

Hence,  the  parameters  Fi  and yu  are  related  by 

f\-ZjuL  exp[(i//41*].  (82) 

To  calculate  the  Mott  distribution  (8l),  starting  withFtjyuis 

needed  from  (82).  In  the  computer  program  this  was  calculated  with 

the  fit  .  _  N  _ 

JJ.-ZO+  0.43  <.G,(m-5o)  J  \ 

^4  =  20  +  (7/lS  )  (R-5-O)  J  SOS  M<2SO  ' 

Figure  11  contains  plots  of  the  Mott  distributions  N(m)  vs.  m17* 
for  various  values  of  the  parameter  jn  .  Note  that  by  the  choice  ofFl, 
one  member  of  this  family  of  lines  is  assigned  to  each  mass  point  ins'. 

The  family  forms  an  envelope  with  some  interesting  properties  (see 
Ref. (8). 

The  input  quantities  wzniJ  the  weight  fractions  of  fragments 
weighing  more  than  M,  (grains)  were  calculated  from 

W(w)s  [(l/2)(m^u)  +  l]  .  (8*1) 


9  See  R.  W.  Gurney  and  J.  N.  Sarmousakis,  "The  Mass  Distribution  of 
Fragments  from  Bombs,  Shells,  and  Grenades,"  Ballistic  Research 
Laboratories,  Aberdeen  Proving  Ground,  Maryland,  BRL  Report  No.  M8, 
7  Feb  19^. 
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IV  RESULTS 

About  12  computations  were  made  with  the  fluid  dynamics  section 
of  the  program  in  order  to  test  the  effect  of  the  grid  size  on  the 
calculated  velocities.  These  were  done  for  metal  cylinders  40  cm 
long  and  10  cm  inside  diameter,  with  uncapped  ends,  filled  with 
Grade  A  Composition  B  explosive,  and  plane  detonated  at  one  end. 
Explosive  to  metal  masses  (C/M)  between  0.1  and  2.0  were  tried.  The 
JWL  equation  of  state,  with  the  explosive  constants  for  Grade  A 
Composition  B  taken  from  Ref.  6,  was  used.  Calculations  were  made 
for  20x3,  40x6,  and  80x12  grids,  e.g.,  in  the  20x3  grid  the  cells 
were  2  cm  in  the  axial  direction  and  5/3  cm  in  the  radial  direction. 
Typical  machine  times  for  the  computation  of  expansions  to  2-3  radii, 
with  the  CDC  6400  computer,  were  45  seconds  for  the  20x3  grid, 

3  minutes  for  the  40x6  grid  and  15  minutes  for  the  80x12  grid. 

Figure  12  shows  the  effect  of  cell  size  on  the  calculated  ve.lscity, 
which  is  given  in  terms  of  the  Gurney  value,  viz., 

Gurney  value  =  (U*+v*)v*  [( l+o.5 

These  computations  were  made  for  C/M=0.2,  a  value  like  that  encountered 
in  projectiles.  The  plots  in  Fig. 12  are  for  various  expansions, R/Rp  , 
where  R.  is  the  initial  inside  radius  (=5  cm  here).  Note  in  the 
figure  that  the  velocities  calculated  with  the  2  cm  long  cell 
(20x3  grid)  are  close  to  those  obtained  with  the  finer  grids.  As  the 
C/M  ratio  is  made  larger,  corresponding  to  thinner  walls  and  more 
rapid  expansion,  the  effect  of  cell  size  becomes  greater.  For 
C/M*2.0,  wall  velocities  calculated  with  the  20x3  grid  are  about  3? 
lower  than  those  calculated  with  the  80x12  grid. 

Sample  weapon  calculations  were  made  for  the  105  mm.  Ml,  and 
the  155  mm,  M107,  projectiles,  with  standard  projectile  steel  casings, 
filled  with  military  grade  Composition  B  explosive  (RDX/TNT/wax, 
59.4/39.6/1.0).  The  JWL  equation  (13)  with  the  following  values  of 
the  constants  was  used: 
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D  =  0.7839  cm/usec  1.634  grams/cc, 

£i=0.08l4  mb-cc/cc  (-  El  fV) , 

02®  0.30  ,03®  4.9055  ,04=  0.058, 

FI  *  4.2,  ,  F2®  0.9  • 

The  above  values,  worked  out  for  a  different  RDX  containing  explosive, 
are  believed  to  be  adequate  for  military  grade  Composition  B.  The 
density  ut>ed  (1.634  g/cc)  may  be  a  bit  low,  but  this  should  not  be 
significant.  The  maximum  reported  experimental  density  for  cast 
Composition  B  is  1.68  g/cc. 

Figure  13  shows  the  calculated  directions  of  the  various  sections 
of  the  105  mm  projectile  at  two  different  times.  Note  that  a  large 
fraction  of  the  metal,  usually  called  the  beam  spray,  appears  in 
the  20th  polar  zone  (95-100°),  and  that  the  directions  are  fairly 
well  stabilized  after  about  60  microseconds  from  the  initiation 
time. 

The  input  and  the  calculated  weight,  velocity  and  angular 
distributions  of  fragments  from  the  105  mm.  Ml,  and  the  1G5  mm, 

M107, shells  are  listed  in  Tables  I  through  XII.  Tables  I  and  II 
give  the  C/M  ratios,  the  inside  radii,  the  scale  factors  X ,  and 
the  values  cf  M  assigned  to  the  mass  points  (see  Sec. III).  The 
initial  grid  and  the  calculated  results  are  Tables  III-VII  for  the 
105  mm  shell  and  in  Tables  V1II-XII  for  the  155  mm  shell.  It  was 
pointed  out  in  Sec. Ill  that  the  calculated  velocities  can  vary  by 
a  few  percent  depending  upon  the  assumption  made  about  when  the 
acceleration  ceases.  For  each  5  degree  polar  zone,  the  calculated 
numbers  of  fragments  in  the  various  weight  ranges  are  listed  in 
Tables  VII  and  XII. 
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TABU  !  f/M  RATIOS,  INSIDE  RADII  (CM), AND  SCALE  FACTORS  X  (IN4^  FOR  1C5  ,AM,  Ml,  AND 
155  MM,  Ml 07  SHELLS. 


105  MM  SHELL,  Ml.  HE  LOADING  DENSITY=1 .634  G3AMS/CC 
C/«  A«0  SCALE  FACTOR 
X«ST*DI*«U/3I/<1*2*C/m> 


K 

1 

2 

3 

4 

5 

6 

7 

8 

9 

to 

!1 

12 

13 

14 

15 

16 
IT 
16 

19 

20 
21 


Z (K*0) 

IRICM) 

0 

2.35 

1.755 

2.775 

3.51 

3»2 

5.265 

3.525 

7.02 

3.65 

8.775 

4.025 

10. S3 

4.2 

t 2.285 

4.215 

14.04 

*>.23 

15.79a 

4.15 

17.55 

UC7 

19.305 

3.965 

21.06 

3.86 

22.815 

3.755 

24.57 

3.65 

26.325 

3.55 

28.08 

3.45 

29,635 

3.36 

31.59 

3.27 

32.9 

3.93 

35.1 

2.55 

C/M 

CHI <lM**4/3 

.123762 

.58116 

.15489 

.575323 

.187499 

.£64955 

.23135 

•504983 

.283879 

.443403 

.315046 

•411023 

.351816 

.37573 

.346186 

.388644 

. £40754 

.398612 

.325526 

,412235 

.310954 

,425875 

.273072 

.483807 

.24122* 

.542416 

.21493 

.598T15 

.192151 

.655058 

.173139 

.708544 

.15631 

•761675 

.137018 

.044104 

.120679 

.925896 

.11263* 

•915206 

7.13284E >2 

1.03668 

155  MM  SHELL,  M 107 .  HE  LOADING  DENSITY  =  i  .634  GRAMS/CC 

C/M  AMO  SCALE  FACTOR 
K»ST«0I»*  Cl/3) / (1 ♦2#C/M) 


K 

Z(Kt0) 

1R (CM) 

1 

4 

2.6 

2 

6,6295 

3.23 

3 

9.259 

3.75 

4 

11.8885 

4.22 

5 

14.518 

4.65 

6 

17.1475 

5 

7 

19.777 

5.35 

8 

22.4065 

5.6 

9 

25.036 

5,81 

10 

27.6655 

6 

11 

30.295 

6.12 

12 

32.9245 

6.12 

13 

3S.554 

6.05 

14 

38.1835 

5,92 

1$ 

40.813 

5.7 

16 

43.4425 

5.45 

17 

46.072 

5.17 

18 

48.7015 

4.9 

19 

51.331 

4.61 

20 

53.9605 

4.35 

21 

56.59 

4.1 

C/M 

CHIUM**4/3> 

.152479 

•336315 

.206232 

.513635 

.241187 

r522 

.265091 

.545188 

,2B9l 12 

.557504 

.306671 

.570282 

.33861 

.549701 

.33072 7 

.602287 

.32352 

.650773 

.31588 

.700387 

.345283 

,641166 

.357498 

.612645 

,336248 

•653534 

.301286 

.730476 

.252666 

.862817 

.209226 

1,01484 

.171087 

1.18497 

•14184 

1.34689 

.107616 

1,63807 

.100756 

1,61515 

9.79935E-2 

1,53253 
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TABU  I!  M  VS  M5.  AVERAGE  WEIGHT  OF  FRAGMENTS  WEIGHING  MORE  THAN  ONE  GRAIN 
ASSIGNED  TO  MASS  POINTS 


NCH.TR  74-77 


TABLE  Hi  105  MM  SH&L  Ml  (NOSE  PLUG).  INITIAL  GRID,  INCLUDING  METAL  OUTSt  BOUNDARY 


POINTS 

K 

Z  (CM) 

*  (CM) 

L 

K 

Z  (CM) 

1  (CM) 

• 

•1.1* 

• 

t 

11 

17. ** 

2.711)1 

1 

t 

t 

t 

It 

It.SIS 

2.1*311 

t 

1.7S* 

• 

t 

11 

tl.tt 

f .57111 

} 

J.S1 

t 

t 

1* 

tt.ti* 

t.Stlll 

* 

liltl 

t 

t 

1* 

t*.*7 

t.tllll 

» 

7.tt 

t 

t 

It 

tt.lt* 

t.ltttl 

* 

1.77* 

• 

t 

17 

tt.tt 

t.l 

7 

lt.»l 

8 

t 

1* 

tt.ti* 

t.t* 

• 

It.tt* 

• 

t 

It 

31.  *t 

t.l* 

* 

lt.O* 

• 

1 

tt 

It.t 

t.ttlll 

1* 

l*.7t* 

t 

t 

tt 

It.tt 

1.7 

11 

>7.»* 

• 

t 

tt 

It.t 

1.7 

It 

lt.lt* 

t 

1 

• 

•1.1* 

1 

1) 

tl.tt 

• 

1 

1 

• 

t.l* 

1* 

tt.ll* 

t 

1 

t 

1.7*5 

t.775 

IS 

10.17 

t 

1 

1 

1.51 

l.t 

I* 

tt.m 

t 

1 

t 

s.tts 

l.St* 

17 

tt.tt 

• 

1 

* 

7.tt 

l.t* 

I* 

tt.ti* 

• 

1 

t 

t.77* 

*.815 

1* 

11. *t 

t 

1 

7 

lt.Sl 

*.t 

It 

st.t 

• 

1 

• 

It.tt* 

*•215 

tt 

1S.I 

t 

1 

• 

It.t* 

«.tl 

tt 

it.t 

• 

1 

It 

15.7*5 

♦  .1* 

• 

•3.1* 

1 

1 

11 

17. *5 

♦.•7 

i 

t 

•7*1)11 

1 

It 

lt.lt* 

s.tts 

t 

!.»** 

.t  tt 

1 

11 

tl.tt 

l.tt 

t 

J.*» 

i.ttttt 

1 

1* 

tt.tis 

).»** 

♦ 

s.tts 

l.|7* 

1 

IS 

t*.s» 

l.*5 

• 

7.tl 

;.t*)i) 

1 

It 

2t.lt* 

1.55 

» 

•  .77* 

l.Stlt7 

1 

17 

tt.tt 

l.*5 

7 

It.*] 

;.* 

1 

It 

t*.tl» 

».lt 

• 

tt.tt* 

!.*•» 

1 

It 

11.** 

1.27 

• 

lt.lt 

l.tl 

1 

tt 

12.1 

l.tt 

!• 

1S.7I* 

1. Mill 

1 

tl 

11  tf 

i.ss 

11 

17.** 

l.l*t*7 

1 

tt 

It.l 

!.*5 

It 

lt.lt* 

1.111*7 

* 

1 

•3.1* 

l.fi 

1) 

tl.lt 

1 .1*1*7 

t 

1 

• 

»  H 

It 

tt.ti* 

1.1SU7 

t 

t 

1.75* 

♦  .t* 

IS 

tt.*7 

l.ttt*7 

t 

1 

3.51 

♦.*!> 

U 

tt.M* 

1.1*111 

t 

t 

s.tts 

».** 

17 

tl.tt 

1.1* 

* 

* 

7.tt 

*.»» 

10 

tt.ti* 

l.It 

t 

t 

t.77* 

5.1*5 

l* 

11. St 

l.tt 

* 

7 

tt.ti 

5.1 

tt 

It.* 

l.ttttl 

t 

• 

It.tt* 

S.3S5 

tl 

It.l 

.** 

* 

t 

l*.t* 

*.J' 

tt 

It.t 

.** 

t 

It 

15.7*5 

5.115 

t 

.1.1* 

I 

t 

11 

17.S5 

s.tt 

1 

t 

i.ssttr 

* 

It 

lt.lt* 

s.tts 

t 

l.VJS 

i.t* 

t 

11 

21. It 

*.t» 

) 

l.»l 

t.issii 

* 

It 

tt.ti* 

*.t7 

« 

s.tts 

t.ss 

* 

15 

t».S7 

*.*7 

s 

7.tt 

I. tot* 7 

* 

It 

2t.lt* 

5.27 

t 

t.  77* 

t.ttlll 

* 

17 

tt.ot 

S.t7 

7 

II. *1 

t.t 

* 

It 

tt.ti* 

5.115 

t 

ll.tt* 

t.t: 

* 

1* 

11. St 

*.♦ 

t 

It.tt 

t.tt 

t 

tt 

It.t 

l.t 

It 

l*.7t* 

t.7ttsr 

* 

tl 

11.17 

*.ts 

* 

tt 

It.* 

*.*7 
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TABLE  V  105  MM  SHELL  Ml  (NOSE  PLUG).  COMPOSITION  B  FILL.  CALCULATED  WEIGHT 

CONTRIBUTIONS  OF  MASS  POINTS  TO  5  DEGREE  POLAR  ZONES,  120  MICROSECONDS 
AFTER  INITIATION 


NABS  NT 

NOLAN  ZONE 

NT  If  Z0NC (9NA1N9I 

1 

1 

239,944 

t 

1 

1911,59 

S 

1 

214.976 

9 

2 

917,499 

9 

9 

917.4SS 

9 

A 

917.499 

9 

9 

966.966 

A 

5 

136.983 

A 

9 

2076.15 

A 

1 

2076.15 

A 

9 

2076.15 

A 

9 

2076.15 

A 

10 

2076.15 

A 

11 

997.443 

9 

11 

1013.75 

9 

12 

1951.12 

9 

19 

1951.12 

9 

1A 

1951.12 

9 

19 

37.7109 

9 

19 

5122.24 

A 

16 

2367.76 

7 

16 

7457.11 

9 

16 

1604.97 

• 

IT 

5761.69 

* 

17 

3677.1 

* 

19 

3242.31 

10 

16 

5229.01 

10 

If 

1776.54 

11 

19 

7129.69 

12 

If 

2940.14 

12 

28 

4340.15 

19 

29 

7340.87 

1A 

20 

7469.93 

19 

20 

7996.45 

19 

20 

9574.79 

17 

20 

9112.04 

if 

20 

9627.23 

If 

29 

10110.2 

29 

20 

10620.2 

21 

20 

11444.9 

22 

20 

996,709 

22 

21 

7421.35 

22 

22 

3019.29 

29 

22 

2261.93 

29 

29 

3913.33 

29 

24 

3913.33 

29 

29 

1376.77 

2A 

29 

3192.34 

2A 

26 

4994.11 

2A 

27 

4994.11 

2A 

29 

997.99 

29 

29 

2910.34 

29 

2f 

3254,1 

29 

90 

1561.76 

29 

90 

327.257 

29 

91 

629.291 

29 

92 

929.2(1 

29 

99 

929.291 

29 

34 

207,234 

2T 

94 

30.0731 

27 

99 

49.249 

27 

96 

40.249 

2? 

97 

0 
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TABLE  VI  105  MM  SHELL  Ml  (NOSE  PLUG).  COMPOSITION  6  FILL.  CALCULATED  FRAGMENT 
WEIGHT,  PERCENT  OF  TOTAL  METAL  WEIGHT,  AND  AVERAGE  VELOCITY  IN  EACH  5 
DEGREE  POUR  ZONE,  120  MICROSECONDS  AFTER  INITIATION 


POUR  ZONE 

MET  WT  (GRAINS) 

PCT  BY  WT 

AV  Va  'FT/SEC) 

1 

2365.07 

1.15783 

1412.25 

2 

917,435 

.449158 

1629,61 

3 

927.485 

.449158 

1629.61 

♦ 

9i  7 .485 

.449156 

1629.81 

5 

995.049 

.46713 

1596.46 

6 

2076.15 

1.01639 

1405.39 

7 

2076.15 

1.01639 

1405.39 

8 

2076.15 

1.01639 

1405.39 

9 

2076.15 

1.01639 

1405.39 

10 

20*6.15 

1*01639 

1405.39 

11 

2011*19 

.994585 

2694.58 

12 

1951.12 

.955177 

3963.04 

13 

1951.12 

.955177 

3963.04 

14 

1951.12 

.955177 

3963.04 

15 

5159.95 

2,52608 

3808.16 

18 

11329. B 

5.54657 

4009.29 

J7 

9639.78 

4.7187 

4480.45 

18 

8471o32 

4.14717 

4760.44 

19 

118*5.4 

5.79895 

4801.33 

20 

87618.4 

42.8939 

4014.62 

21 

7421.35 

3.63315 

3398.29 

22 

5281,19 

2.58543 

3116.71 

23 

3913.33 

1.96683 

2740.86 

24 

3813.33 

1.86683 

2740.86 

25 

4568.11 

2.23633 

2136.49 

26 

4994.11 

2.44489 

1876.01 

27 

4994.11 

2.44489 

1876.01 

28 

3599.33 

1.76158 

2326.3/ 

29 

3254.1 

1.59306 

2496.82 

30 

1889.02 

.924775 

2704.8 

31 

529.261 

.306057 

3697.35 

32 

629.261 

.308057 

3697.35 

33 

529,261 

.308057 

3697.35 

34 

238.807 

.116909 

4263.84 

35 

46.248 

2.26409E-2 

8065.03 

36 
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TABLE  VII  105  MM  SHELL,  Ml  (NOSE  PLUG).  CALCULATED  NUMBERS  OF  FRAGMENTS  IN  VARIOUS 
WEIGHT  RANGES  IN  EACH  5  DEGREE  POLAR  ZONE,  120  MICROSECONDS  AFTER 
INITIATION 


40.  Of  7843M6NTS  19  »ot«*  20965 
#7.  awmi  in  a*«i  4S 


fOUtZONC 

1 

t 

1 

* 

r 

■ 

* 

it 

n 

i« 

u 

i* 

i* 

1* 

IT 

1* 

1* 

I* 

tl 

ti 

tl 

I* 

II 

11 

IT 

I* 

I* 

1* 

11 

II 

31 

1* 

is 

i* 

•>t 

i 

i 

i 

* 

9 

0 

T 

• 

* 

19 

11 

II 

11 

l* 

IS 

l* 

IT 

!• 

1* 

It 

11 

11 

tl 

1* 

IS 

IS 

IT 

IS 

IS 

IS 

11 

31 

11 

IS 

IS 

IS 


l  - 1 

2  •  S 

5  -  10 

I*  -  15 

90k*8  tOHt 

150  -  250 

07.  TK*N  | 

.5051*1 

.919152 

.905511 

.47951 

1 

1.13157 

2.4*149 

.1*5977 

.104404 

.374501 

.241411 

2 

.41*97* 

.9*2725 

.195*77 

.104404 

.374501 

•2*1011 

1 

.438*74 

.902725 

.1*5*77 

.104404 

.17*501 

.201011 

4 

.418*74 

.902725 

.21)1*1 

.401101 

.404**1 

.141103 

5 

.511*14 

1.84*9 

I.1SSS* 

2.11151 

2.09129 

1.4212 

0 

1.52111 

2.19*22 

1.1055* 

2.11141 

2.0*129 

1.4212 

7 

1.52111 

2.1*022 

1.10559 

7.11151 

2.0*129 

1.4212 

* 

1.52111 

2.1*022 

i.issss 

2.11141 

2.0*129 

1.4212 

* 

1.52111 

2.1*021 

1. 10559 

2.11141 

2.9*129 

1.4212 

1* 

1.52111 

2.1*821 

1.I2SII 

2.14940 

2.02505 

1.17*71 

11 

1.47151 

2.1214* 

1.09539 

1.90410 

1.90514 

1.11501 

12 

1.42*5 

2.88*11 

1.0951* 

1.99410 

1.90514 

1.31501 

11 

1.42*5 

2.08*11 

i.o*si* 

1.9*430 

1.9*51* 

1.31501 

14 

1.42*5 

2.05011 

1.11 750 

5.01-.1S 

5.20110 

3.7471* 

15 

1.87214 

5.1*117 

11.7147 

19.01* 

18*2907 

12.102* 

10 

9.1*047 

10.7511 

II.IM) 

31.5107 

29.0112 

18.1922 

17 

0.4*289 

0.01*5* 

It. C0*T 

40.0710 

41.0410 

24.924 

18 

0.73221 

3.9*82 

IT. IZSO 

02.1108 

55.1040 

31,7022 

1* 

9.59132 

5.7117* 

74.0103 

111.227 

124.12* 

01.1720 

28 

01.7791 

04.22*5 

0. SOOTS 

15.1147 

14.5411 

9.559*3 

21 

0.400*9 

0.555*1 

0. 101S7 

10.4913 

18.3839 

0.8027 

22 

4.5971* 

4.004«t 

4.40352 

7.8092 

7.497*0 

4.91194 

tl 

3.119* 

1.10*17 

S. 40152 

7.0092 

7.*97*0 

4.9119* 

I* 

1.119* 

3.10*17 

I*. 0041 

19.102* 

12.0191 

19.220* 

IS 

2.01998 

1.777*5 

10.0057 

50.7100 

45.0*00 

15.7128 

20 

2.51011 

.8**411 

10.0097 

SO. 71*0 

45.8*00 

25.7125 

17 

2.51*11 

.08*411 

15.0244 

24.47 

20.8480 

12.1078 

18 

2.0*07 

1.47010 

•.*1524 

10.5119 

14.0727 

8.49701 

2* 

2.000 

1.023*8 

S. 19702 

9.00704 

9.10724 

*.***78 

1* 

1.50718 

.9*585 

1.10242 

2.21441 

2.8*04 

1.29108 

11 

.553141 

.417*1* 

1.10242 

2.21*51 

2.0*8* 

1.2*108 

It 

.551141 

.417*1* 

1.10242 

2.21*41 

2.0*84 

1.2*108 

11 

.55314} 

.41701* 

.494271 

.8*8120 

.777177 

.440*58 

1* 

.20**2 

.15820 

9.572216-1 

.10*25 

.1509*9 

9.508*2242 

15 

4.00517C-2 

1.804*16' 

••57221C-I 

.10*25 

.150549 

9.50B02E-2 

10 

4.805176*1 

1.804916' 

IS  -  IS 

24  -  SO 

50  -  10* 

too  -  IS* 

.970451 

1.52127 

1.05029 

.945080 

.379572 

.590151 

.042529 

.1740*0 

.179572 

.49*151 

.0*2529 

.1740*0 

.179572 

.49*151 

.442529 

,1740*0 

.405942 

.7*2304 

.787*00 

.4*70* 

1.90055 

2.00311 

2.80792 

1.45432 

1.90055 

2.80111 

2.50792 

1.40412 

1.90955 

2.00311 

2.80792 

1.45412 

1.90955 

2.00311 

2.80792 

1.45412 

1.90055 

2.00111 

2.80792 

1.45*32 

1.90095 

2.77154 

2.72000 

1.4185 

1.94999 

2.09009 

2.03881 

1.37013 

1.04999 

2.09*09 

2.01881 

1.17011 

1.04999 

2.09*04 

2.01881 

1.37013 

5.1995 

7.51010 

7.30135 

3.70571 

10.2745 

22.510* 

20.430 

9.82177 

13.0071 

30.1040 

24.1*50 

10.1090 

10.5105 

10.027 

25.5000 

*.1*1*8 

41.5151 

49.4951 

35.52*1 

11.177* 

IOT.911 

140. 0*T 

111.589 

*1.1710 

12.7171 

17,2591 

15.1841 

7.042*8 

9.04971 

12.2021 

10.8009 

9.011*1 

0.51441 

8.80010 

7.8031* 

1.018* 

0.53441 

4.95010 

7.8031* 

3.018* 

20.020 

22.0*44 

13.14*1 

4.20102 

24.0954 

29.414* 

10.47*2 

4.021*4 

28.9954 

29.418* 

10.47*2 

4.023*4 

14.0459 

10.5298 

11.022* 

1.821*2 

11.1117 

11.1512 

9.0700 

1.0207* 

0.2922 

7.49705 

5.50*71 

2.0*51* 

1.05312 

2.4*401 

1.05*71 

.082108 

1.05112 

2.09*01 

1.05*71 

.082188 

1.05112 

2.0*401 

1.05*71 

.002108 

.027100 

.794914 

.02797# 

.258801 

.121497 

.151945 

.12101* 

•85*112 

.121497 

.isms 

.12101* 

.0501)2 

40 
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TABLE  VIII  155  MM  SHELL,  M107  (M51A5  FUZE).  INITIAL  GRID,  INaUDING  METAL  OUTER 
BOUNDARY  POINTS 


9 


K 

Z  (CM) 

R  (CM) 

L 

K 

Z  (CM) 

R  (CM) 

• 

•3 

0 

2 

11 

30.209 

4.00 

1 

» 

0 

2 

12 

32.0249 

4.0* 

2 

•.*299 

• 

2 

13 

35.554 

9.03333 

3 

9.290 

0 

2 

19 

30.1*35 

3.949*7 

* 

11.9095 

0 

2 

15 

40.013 

3.0 

9 

<9.519 

0 

2 

19 

93.9425 

3.93333 

* 

17.1475 

0 

2 

IT 

49.072 

3.9940/ 

r 

10.TTT 

• 

2 

10 

40.7015 

3.249*7 

• 

22.9099 

• 

2 

10 

51.331 

3.07333 

* 

29.039 

0 

2 

20 

53.9*05 

2.9 

10 

2T.999S 

0 

2 

21 

59.5* 

2.73333 

ti 

30.299 

9 

2 

22 

90.95 

2.73333 

iz 

32.9299 

0 

3 

0 

•3 

2 

13 

35.599 

0 

3 

1 

4 

2.* 

19 

30. 1035 

0 

3 

2 

6. *2*5 

3.23 

IS 

90.013 

0 

3 

3 

9.259 

3.75 

19 

93.9925 

0 

3 

9 

11.0095 

9.22 

JT 

99.072 

0 

3 

5 

14.510 

4.95 

1* 

99.7015 

0 

3 

* 

17.1475 

5 

1* 

91.331 

0 

l 

7 

19.777 

5.35 

20 

93.9*05 

0 

3 

• 

22.4095 

5.9 

II 

59.99 

0 

3 

9 

25.039 

5.01 

II 

00.95 

0 

3 

10 

27.9*55 

9 

• 

•3 

.****»/ 

3 

11 

30.295 

9.12 

1 

9 

.39*** J 

3 

li 

32.9245 

9.12 

t 

9. *2*9 

1.079*7 

3 

13 

35.954 

4.05 

3 

9.259 

1.25 

3 

19 

30.1*39 

5.92 

* 

11.0005 

1 .909*7 

3 

IS 

40. *13 

5.7 

S 

19.910 

t.SS 

3 

19 

43.4425 

5.45 

* 

17.1975 

1.599*7 

3 

17 

49.072 

5.17 

T 

19.777 

1.70333 

3 

10 

40.7015 

4.9 

• 

22.9095 

1.0***7 

3 

10 

51.331 

4.41 

* 

2S.039 

1.939*7 

3 

20 

53.9*05 

9.35 

It 

27.4*55 

2 

3 

21 

99.59 

4.1 

u 

30.295 

2.09 

3 

22 

90.45 

4.1 

12 

32.0295 

2.09 

9 

0 

•3 

2.3 

13 

35.5S9 

2.019*7 

9 

1 

4 

4 

1« 

30.1*35 

1.97333 

* 

2 

9.92*5 

4.5* 

IS 

90.013 

1.9 

• 

3 

9.259 

5.12 

t* 

93.9925 

1.01*67 

9 

9 

11.9*05 

5.94 

IT 

99.072 

1.72333 

9 

5 

14.519 

9.1 

IS 

90.7015 

1.93333 

9 

* 

17.1475 

9.90 

1* 

51.33; 

1.53*97 

9 

7 

19.777 

9.0 

it 

53.9*05 

1 .95 

9 

• 

22.4*95 

7.15 

21 

59.59 

1.399*7 

9 

0 

25.039 

7.45 

It 

90.95 

1.399*7 

9 

10 

27.9*55 

7.73 

• 

•3 

1.33333 

9 

11 

30.2*5 

7.75 

1 

9 

1.73333 

9 

12 

32.9245 

7.7 

2 

*.*299 

2.15333 

13 

35.554 

7.7 

3 

9.259 

2.5 

9 

19 

34.1*35 

7.7 

• 

11.0009 

2.01333 

9 

15 

99.013 

7.7 

9 

<9.910 

3.1 

9 

19 

93.9925 

7.7 

t 

17.1975 

3.33333 

9 

17 

49.072 

7.7 

T 

19.777 

3.59*97 

9 

10 

40.7015 

7.’ 

• 

22.9095 

3.73333 

9 

It 

51.331 

7.9 

* 

25.035 

3.07333 

9 

20 

53.9*05 

7.42 

10 

27.9*59 

9 

9 

21 

59.59 

7.25 

A.  « 
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TABLE  IX  155  MM  SHELL,  Ml 07  (M51A5  FUZE).  COMPOSITION  B  FILL.  CALCULATED  SCALED 
MASSES,AND  POLAR  ANGLES  121  MICROSECONDS  AFTER  INITIATION 


m  (MASS/2  ir,  GRAINS) 
4.05502 
32.4402 
64.8804 
226.854 
109.105 
125.50/ 
143.915 
160.392 

174.13 

183.51b 

203.406 

224.10/ 

240*752 

233.52b 

227.164 

234,781 

251.239 

276,638 

305.338 

335.43 

367.824 

414.098 

504.858 

446.744 

211.511 

56.5234 


49(H54l/2) 

0 

3.37265 

25.2768 

54.7868 

73.9444 

82.9255 

83,9218 

85.0536 

86,6071 

89.0314 

90.7568 

91.1007 

91.9013 

94.527 

97.9178 

100.446 

102.447 
103.568 
103. B74 
103.914 
103.911 
1C3.122 
113.303 
148,9 
173. B42 
177.269 
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TABLE  X  155  MM  SHELL,  M107  (M51A5  FUZE).  COMPOSITION  B  FILL,  CALCULATED  WEIGHT 
CONTRIBUTIONS  OF  MASS  POINTS  TO  5  DEGREE  POUR  ZONES,  121  MICROSECONDS 
AFTER  INITIATION 

MSS  FT  POL*«  ZONE  * T  IV  70V£<GRA1NS> 


1 

1 

3^2.84 

i 

1 

3142.72 

3 

1 

466,972 

3 

2 

1434.76 

3 

3 

1434.76 

3 

4 

1434.76 

3 

5 

!  <*34,76 

3 

6 

79.4232 

♦ 

b 

3517.5 

♦ 

7 

3723.73 

♦ 

8 

3723.73 

4 

9 

3723.73 

4 

10 

3723.73 

4 

11 

3564.96 

5 

11 

117.628 

S 

12 

2758.65 

S 

13 

2758.65 

S 

14 

2758.65 

s 

15 

2176.22 

1 

IS 

1445.23 

6 

16 

6845.27 

6 

17 

3965.28 

7 

17 

13942.1 

» 

17 

14802.6 

• 

10 

735.806 

9 

18 

16869.2 

10 

18 

17778.6 

11 

18 

11062.2 

11 

19 

8643.27 

12 

19 

21710.9 

13 

19 

23323.4 

14 

19 

22623.3 

15 

19 

3163.35 

IS 

20 

18844.1 

14 

20 

1 9685.1 

14 

21 

3859.78 

17 

21 

24339.4 

14 

21 

26799.9 

19 

21 

29580.3 

20 

21 

32495.5 

21 

21 

35633.8 

22 

21 

40116.7 

23 

21 

10808.9 

23 

22 

29779.3 

23 

23 

19000.7 

24 

23 

2069.65 

24 

24 

6096.17 

24 

25 

6096.17 

24 

26 

6096,17 

24 

27 

6096.17 

24 

28 

6095.17 

24 

29 

6096.17 

24 

30 

4633.16 

25 

30 

991.056 

25 

31 

4091.29 

25 

32 

4091,29 

25 

33 

4091,29 

25 

34 

4091.29 

25 

35 

3143,64 

24 

35 

1950.23 

24 

36 

3625.6 

27 

36 

694.48 

27 

37 

0 

TOT  MET4L 

«Ta  81.1606  l„S( 

568264.  98AIN5) 
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TABLE  XI  155  MM  SHELL,  M107  (M51A5  FUZE).  COMPOSITION  B  FILL.  CALCULATED  FRAGMENT 
WEIGHT,  PERCENT  OF  TOTAL  METAL  WEIGHT,  AND  AVERAGE  VELOCITY  IN  EACH  5 
DEGREE  POLAR  ZONE,  121  MICROSECONDS  AFTER  INITIATION 


POLAR  ZONE 

MET  WT  (GRAINS) 

aCf  8Y  WT 

AV  VEL  (FT/ 

1 

4002.53 

.704343 

1159.48 

Z 

1434.76 

.252481 

1304.8 

3 

1434.76 

.252431 

1304.8 

♦ 

1434.76 

.252491 

1304.8 

5 

1434.76 

.252491 

1304.8 

6 

3597.02 

.632984 

1226.37 

7 

3723. 73 

.655282 

1224.6 

8 

3723.73 

.655292 

1224.6 

9 

3723.73 

.653282 

1224.6 

10 

3723.73 

.656232 

1224.6 

11 

3682.59 

.645041 

1341.57 

12 

2758.65 

.485o53 

48B6.54 

13 

2758,65 

.463453 

4886.54 

1 A 

2758.65 

.486453 

4886.54 

15 

3621.46 

.637294 

4694.05 

16 

6845.27 

1.20459 

4404.2 

ir 

32612.9 

5.73904 

4168.55 

18 

46445.9 

3.17327 

4385.27 

19 

79464.3 

13.9837 

4405,34 

20 

37729.3 

6.63938 

4810.68 

21 

203634. 

33.8345 

3594.75 

22 

28779.3 

5.06442 

2245.96 

23 

21078.4 

3.70926 

2233.1 

2* 

6096.1 7 

1 .07277 

2115.02 

25 

6096,17 

1.07277 

2115.02 

26 

6096,17 

1.07277 

21 15.02 

27 

6096,17 

1.07277 

2115.02 

28 

6096.17 

1.07277 

2115.02 

29 

6096.17 

1.07277 

2115,02 

30 

5615.02 

.9881 

2206.7 

31 

4091,29 

.719962 

2639.32 

32 

4091.29 

.719962 

2639,32 

33 

4091.29 

.719962 

2639.32 

34 

4091.29 

.719962 

2639,32 

35 

4993.87 

.878754 

3372.22 

36 

4310. uB 

.753465 

4615.4 

1> 


•  .- 

.'65' 


F  ‘  # 


S  V  « 


;  v? 


p  •  V 


4  ; 

jr  ci 


l 

i 


jg.  ^ 


I  ^ 
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TABLE  XII  155  MM  SHELL,  M107  (M51A5  FUZE).  CALCULATED  NUMBERS  OF  FRAGMENTS  IN 
VARIOUS  WEIGHT  RANGES  IN  EACH  5  DEGREE  POLAR  ZONE,  121  MICROSECONDS 
AFTER  INITIATION 


«o(  90  M*a«c«rs  is  *ow*»  toots 

if*  MOMS  In  n<u Jn$ 


1  •  t 

f  •! 

5  -10 

IS  •  1 

•«*49»4 

t.*o*fi 

!.*H0/ 

l.ts 

.  J06*t>* 

,?0*tf64 

.s«bTjr 

•411134 

«  30**^*1 

•S»S7)7 

•4)113* 

.J«44»>1 

.sssrj/ 

•41(110 

•JU4444 

.iisnr 

•  4J  ft 94 

•  *s»  7.6 

4.4f*>N# 

I .*• T1S 

/.»*?•»% 

4.7*479 

4.4609/ 

3.1 f t40 

t./tVM 

*.**•4/ 

J.ttr** 

i.MfM 

t.H.N 

4.46*9/ 

J.lt»4fc 

f.*4fj* 

••*•4/4 

4.4609/ 

3.1*740 

f.OfOjt 

*.//v/4 

♦.WS*. 

).l»*4i 

9,  HU. 

4.1#»7J 

-,04/*%. 

t.tLOSf 

4.1  n 

4.L  *44 

f.TlSSI 

t'34J40 

4.0*044 

t.tOOSf 

Mkhl 

N. JO604 

6.16*64 

3.44*11 

5,fm7 

US4.< 

4.f *0l( 

2/. 66*4 

40.M1O 

.  *.>#:  * 

f 4 i 7*04 

ie.it** 

6J.744# 

*o.«r*j 

33.01 #i 

Jf.7>J» 

*>».#I6* 

60.011 

*t.flJ9 

14.V6J4 

JO. •/»* 

)0.«30> 

ft. 14*1 

M.  IPli 

4M4«| 

IT.it/ 

II. Of 

n.ofct 

»4.»;as 

9*  1*1# 

)f • 1 3*1 

'!*•' rs 

?1.bgk 

ll.*4*» 

If. 31» 

U.S»  *0 

ll.SSni 

19.044* 

It*  633 

If. 4176 

t i.ssas 

IM44S 

1I.S33 

1I.S1TS 

II  .MSS 

1 f .*44* 

1?. 333 

If.SI ts 

II. t  ns 

;t.#44* 

7 4.133 

12.4t»t 

ft. 46*0 

\*  .*44* 

)#.*33 

10.41J6 

1 7*#44| 

16*31? J 

1S.4SSS 

MtHJ 

S.(U*?J 

6.0ft 

4.41JST 

I.ITSU 

6.  Off 

*.4L  JS7* 

3.47SSJ 

k.«kM 

6.0tf 

4*vi)«: 

l.l’SSI 

6.0ft 

4.0) JOT 

l.rosiy 

*.7# 

4.6111 1 

*.4J JOT 

l,MSM 

4.70*69 

♦.44T«» 

3.14011 

IS  -  IS 

41  •  SO 

S5  -  100 

ISO  -  > 

MISS! 

(.INS) 

r.HJoi 

1. 0J40f 

.S»JSM 

.41/174 

1 . 0047# 

•  MISS) 

.SUSTJ 

••lists 

l •sssts 

•  MJM1 

NOCOT* 
I.8S4  * 

i.MJ?* 

S.MllO 

I.Mlli 

s.isto* 

4.7040# 

».mn 
i.im* 
u.iit* 
so.  our 

IMlTf 

UM1« 

♦o.m  i 

10.0307 

u.ojor 

14.0307 

IM)t> 

\?'**or 

u.«3lf 

u.mi 

I.MM 

r  mi 
’••MS 
I.ItTH 
S.1SS77 


PfLM  IONC 

I 
t 
S 
• 
s 
s 
7 
• 

• 

IS 

II 
II 
13 
If 
IS 
I* 
l» 

IS 

IS 

IS 

II 

If 

II 

l« 

IS 

IS 

if 

if 

is 

is 

ii 

si 

>3 

3* 

31 

IS 


IS#  -  III 

ST.  T«4* 

l.4)S0f 

4.10000 

•SSS4SS 

1.S0SSI 

•ssooss 

I.SOSS. 

•SSS46S 

1 .sossi 

.4*4*46 

i.soss: 

0.SSS4 

3.OS007 

IrSlOOS 

3.70717 

f.#306? 

l.tSMI 

*.#3«’tS 

3.7S7I7 

I.Dkk) 

3.7071/ 

f.##fS 

3. 70S* 

f. *7/140 

«.SHO 

t. 1/440 

|.S0«S 

t.ifost 

I.5S40 

I.SSirff 

3.060)0 

S.SIIII 

4.L  JC’77 

•S. *6*6 

Si.iirs 

34.004 J 

SI .SOSf 

ss.uir 

M.4«m4 

♦  0.6*10 

SS.4MS 

IDV.173 

(1.0604 

3S.OOSS 

IS.SS*} 

11.131? 

4.04 

S.JSS7S 

♦  .04 

S.3SS7S 

4.S4 

S.3S07* 

4.04 

s.issts 

4.04 

S.ISStf 

4.94 

4.0*007 

4.01061 

3.30700 

3.0047) 

3.167*0 

3.00073 

3.10709 

?.Oht3 

3.307*0 

4.00073 

3.07440 

'.SS300 

3.134*4 

4.S0I0 

«  15  + 


I 
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Appendix  A  -  The  Fragment  Prediction  Code 

Table  AI  is  a  complete  list  of  the  computer  program,  in 
BASIC.  A  set  of  notes,  which  explain  what  the  various  statements 
do  ,is  m  Table  All.  Table  AIII  is  a  list  of  variables,  together 
with  the  quantities  needed  in  the  dimension  statements. 

Figures  A1-A3  illustrate  descriptions  in  the  notes.  The  particular 
weapon  treated  in  the  program  list  (Table  AI)  is  the  105  mm.  Ml 
shell.  For  other  items,  which  do  not  have  a  curved  base,  the 
input  is  somewhat  simpler  and  some  of  the  statements  can  be 
removed.  Tahle  AIV  is  a  list  of  the  input  statements.  To  run 
the  program  for  a  particular  weapon  one  proceeds  through  this  list, 
providing  the  necessary  data  at  the  listed  statement  numbers  and 
removing  cr  bypassing  any  statements  that  do  not  apply.  Parts  of 
the  program,  for  example,  the  initialization  to  get  values  of  the 
scale  factor  X  ,  or  the  fluid  dynamics,  can  be  run  by  inserting 
appropriate  bypass  or  exit  statements.  Output  statements  are 
listed  in  Table  AV.  Since  all  the  statements  in  BASIC  are  numbered, 
any  of  the  output  can  be  bypassed  by  inserting  a  GO  TO  statement 
at  the  appropriate  place. 

It  will  be  seen  in  Tables  AI  and  AV  that  there  are  provisions 
for  output  in  formats  need'-  by  lethal  area  programs.  For  the  JMEM 
lethal  area  program,  the  ne c'jsary  calculated  quantities  are  put 
on  tape,  from  which  a  c  'eck  In  the  appropriate  format  is  prepared. 
The  FORTRAN  program  which  determines  how  this  card  deck  is  made  is 
listed  in  Table  AVI.  For  the  AMSA4  lethal  area  program,  the 
necessary  input  quantities  are  listed.  Provisions  for  getting  a 
card  deck  for  this  program,  and  for  machine  plotting  of  any  of  the 
output,  could  be  added  by  putting  the  needed  material  on  tape  and 
using  auxiliary  FORTRAN  routines. 
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TABLE  A-l  BASIC  PROGRAM  LIST-FRAGMENT  PREDICTION  CODE 


01000 
C  *  v  1  0 
01020 
01030 
01040 
01050 
01060 
01070 

oioeo 

010*0 

01100 

onio 

01120 

01130 

01140 

01150 

01160 

01170 

01180 

01190 

01200 

01210 

01220 

01230 

01240 

01250 

01270 

012BO 

01290 

01300 

01305 

01310 

01320 

01330 

01340 

01350 

01360 

01370 

01380 

01390 

01400 

01410 

01420 

01430 

01440 

01450 

01460 

01470 

01480 

01490 

01500 

01510 

01520 

01530 

01540 

01550 


PRINTiloS  MM= 

P9«3. 14159265 

P6*2*P9 

Q9»P6#15.4185 

P7»5*P9/1B0 

P8»P6/454 

REM  *0.  OF  STERADUNS  IN  POLAR  ZONE 
FOR  J6*l  TO  18 

Y8«J6)sP6*(C0S( <J6-l)PP7)-C0S(J6*P7l) 

Y8(37-J6)*rUiJ6> 

NEXT  J6 

J9=10 

T7»40 

T4«150 

Tl*0.3 

E8»l.l 

N4*20 

N5»7 

N6*200 

N7*l 

K4*0 

R2«4 

R4=7.84 

C3*2 

E7*lE-l0 

01=20 

R0»0 

Z0320 

Ll*3 

*1=21 

REM— ———BRING  IN  EQUATION  OF  STATE  CONSTANTS 
G0SU8  08*50 

PRINT5K1==KI *rLl“5Ll *EK4*=K4 
PRINTS J9»sj9 

PRlNTST7«ST7fST4*5T4,STlwST1.5E8*=E8.SN5*SN5 

PRINT5N6*=N6.SE7*=E7 

PRINTSR2«5R2*SR4*SR4 

PRINTSR0=5R0  »?20*£20  »S01 sSDlrFTS 

DIM  N3< 37 *101*97(37*101 *Ml (27) *M(10> »N(27»10> 

OIM  *6 (271  * A9 ( 27) « 33 (27) 

DIM  *(7(37*31  )»*(2(27»  10)  *85(37, 10)  *88(37)  *89 (37) 

OIM  N8(37) *YB( 37) «Y9 (37) *M6( 37* 10) *06(37*10) 

OIM  R(23»7) *Z ( 23 1  7 ) »R1 (20*5) *Z1 (20*5) 

OIM  R2 (23*7) 

OIM  81 (20*5) *8(23*7) 

DIM  83(23*7) 

OIM  U(23»7) »V (23*7 ). VI (20*5) *V5(20*5) 

OIM  01 (20*5) *T3 (20 *5) *02 (20.5) *P2(20»5) .22(20*5) 

OIM  Y (36) »Yl (36) »F (20)  * U7 (21 ) » V7 (21 ) 

DIM  SB(21)*R3(21)*Z3(21)»P3(21) 

OIM  R4(21)*24(21)*Al(20*5)*X9(2l) 

REM— — — — — —  i nPUT  ZIK1.C)  *Z(0*0)  «zn*0>*Z*Kl*l*9> 

Z(K1, 01*35.1 
Z (0*  0) *«3. 15 
Z(1*0)«0 
ZCKl*l.0)*3b,9 
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01560  REM— - - - - - - - INITIALIZE  Z(K,0>  — - — ' 

01570  FOR  K«1  TO  K1 

01580  Z(K»0»«Z<1,0>*(Z<K1»0)-Z(1*0>  >»<K-1)/<<1-1) 

01590  NEXT  K 

01600  REM - - - INPUT  SPECIAL  Z  UiO)  VALUES-———--' 

01610  Z (20*0) *32,9 

01620  REM  —————— —INPUT  R »K,L1  >  »R  >  — — — — — 

01630  GOSUS  09960 

01640  REM - - INPUT  SPECIAL  R <K »Ll > ♦» <K»L1 *1 >  VALUES 

01650  GOSUS  1 00  TO 

01660  REM - — - - - - - INITIALIZE  Z  *K,L»  »R<K*L> 

01670  FOR  K*o  TO  KIM 
01680  FOR  L*1  TO  L1M 
01690  Z (K,L) &Z (Kf 0) 

01700  NEXT  L 

01710  FOR  L* 1  TO  Ll-1 


01720  R(K,L)«R(K»Li; »l/li 
01730  NEXT  L 
01740  NEXT  K 

01750  REM— - - - - INPUT  SPECIAL  Z<K1,L>  VALUES 

01760  Z  (K1  »2 )  »34,t>2 
01770  i (K1.L1 1*33.87 
01780  (K1'L1«1)3Z(K1.L1> 

01790  REM— —————WEIGHT  DISTRIBUTION  INPUT— 

1800  GOSUB  2l2o0 

01910  REM— —PRINT  CASE  DIMENSION  INPUT— - - 

01520  PRINT  =L=»=65«iZ5t5R= 

01930  for  l*o  tu  lim 
01940  FOR  KaO  TO  KX#1 
01850  PRINTL.K.Z(K.L) «R(<*L> 


01960  NEXT  K 
,11970  NEXT  L 
J 1 975  REM 

01980  REM— —————— INITIALIZE  R3  U>  »Z3  PO  — — — 

01990  GOSUB  110*0 

01900  REM— —————SET  Z  for  lighting——— 

01910  Z*Z (Kfc«l » 0) 

01920  REM - - INITIALIZE  ME  AND  FUZE  DENSITIES—— 

01930  G0SU9  07190 

01960  REM— - - - — CDORO  OF  ZONE  CENTERS, ZONE  MASSES 


01970  FOR  K*1  TO  *1-1 
01980  FOR  L«0  TO  Ll-1 

01990  Zl  <K*L)«(Z<K,L>*Z<<M*L>*Z(K*l»L*l>*Z<<.LM>  )/A 

02000  Ri (K»L)a(R(K»L)^RIK^l*L)^R(K*ltL*i>*RIK,L*l) )/♦ 

92010  GOillU  07300 

02020  GOSUS  07390 

02030  Ml (K»L) *N5 

02040  NEXT  L 

02050  NEXT  K 

02060  rem— —interior  grid  point  masses——— 

02070  FOR  K*z  to  Kl-1 
02080  FOR  L« 1  TO  Ll-1 
02090  Gl»Zl  «<-l»L-D 
02100  H1*R1 <<-l»L-l> 

02110  52»Z1 (K.L-ll 
02120  »*2»Rl  (K.L-l) 
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02130  G3»Zi  «iU 
1)2 1  40  H3a«i  K,L> 

02550  G4=zi 
02160  H4=R1(«-1»L) 

02170  30S08  07390 
02180  H(*.U*WS 
02190  NEAT  L 
02200  NEXT  K 

02210  RE* - HASSES  OF  INTERIOR  AXIS  POINTS— 

02220  FOR  <*2  TO  M-l 
02230  L*0 
02240  G0SU8  08090 
02250  S0SJ3  07J90 
022N0  «(K*Dsw5 
02270  G0SU8  08180 
02280  GOSOd  07390 
02240  W(*«L)*W<*»L)425 
02300  NEXT  k 

02310  R£8— — HASS  CORRECTION  AT  FJZE-HE  BDuNORr— — — — 

02320  IF  K4»i)  THEN  02370 

0233*1  K*X*M 

02340  FOR  L*0  TO  Ll-1 

02350  *MK*Ll>(R2/R3«ll*«MK«U/2 

02360  NEXT  L 

02370  REN 

02380  REN  —————— ——BAS  MASSES  ASSOC  HIT*  NETAL  POINTS- 

02390  FOR  X=2  TO  <*i-l 

02400  L»'  1 

*32410  G0SJ8  08270 

02420  GOSOd  0 739J 

02430  »S 

02443  GOSUS  98360 

02*50  GOSOs  37390 

Ql*- CO  •t3«<»L)«rA3(8.»Ll  ♦HS 

02*70  NEXT  K 

v'i4H0  L«L1 

02*93  CF  K4«0  THEN  02520 
Oeb.03  *{**♦♦  1 

025)  P  83 (KtL)  =  ( l *R2/R31 *43 (H»L) 72 

3252;  FOH  tci  TO  Ll-I 

02533  «*1 

02540  GOSOd  08090 

02550  GOSOd  97390 

0l560  «ilK»L)*25 

02579  GOSOd  08360 

02‘80  GOSOd  07390 

025V0  4J(K»L)»R31<^LJ48S 

02600  NEXT  L 

026 1 J  cOR  L«1  TO  Ll-1 

02579  *»*1 

326.1)  GOSOd  08180 

02543  GOSOd  J7390 

02650  8J«X»l:«<»3 

02660  G0SJ5  08270 

02670  J3S03  07399 

1)2680  *5<X»ul  s*3(ft«L)  *46, 
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0269.1 
02700 
02>  10 
02720 
02730 
Ot  7*0 
02750 
02760 
02770 
0’780 
02790 
02000 
02910 
02920 
02930 
029*0 
02950 
02960 
02970 
02980 
02990 
02900 
02910 
62320 
02930 
029*0 
02950 
02360 
02970 
02990 
02990 
03000 
53010 
93'JeO 
03039 
030*0 
03050 
04060 
03o  70 
03080 
03090 
03100 
03110 
03120 
031  aO 
031*0 
03133 
03160 
031  70 
03199 
03190 
03200 
03210 
03220 
03230 
032*0 


NEXT  L 

5*1 

L*0 

GQS08  08090 
30SJ9  17390 
*»3(K.L)*85 
L*L1 

GOSOS  08360 

G0SO8  07390 

*3<5tL)*»5 

5*51 

L*0 

GOSOS  08180 
GOSOS  07390 
83<5«U=*5 
L*Ll 

G9S08  09270 
GOSOS  07390 
*3(K»L>«*5 

GEM - - - - - METAL  NASS  POINTS— 

FOG  L«0  TO  U1 
G2(0»L>*8* 

G2'51fU»«* 

NEXT  L 

F 08  5*1  TO  51-1 
G2<K.LI >»** 

NEXT  K 

K*0 

L*0 

GOSJG  07>i0 
90S J9  0739* 

9(l»0>«<MlvO)«93*9* 

FOG  L*1  TO  U-l 
5*0 

30S03  077J0 
SOSOS  07390 
*(l«L)*tftl*l.)«23«M 
GOS 0907820 
G0S09  37399 
«m«U)sw()  .}.>♦* 3»«* 

NEXT  L 
AmQ 

L*U 

G0SO8  1TJ00 
60SJ9  07390 
«  ( 1  *L1  >  *3*  <  1  »L.  I  -««3**>* 

GOSOS  37320 

?0S08  07390 

*11  tLl  !*aMl»Ui)*#3*NA 

5*1 

GOSOS  08000 
GOSOS  07390 
tf(l*Ll>3«U»Ll>*«S 
FOG  5*2  TO  51-* 

L*L1 

GOSOS  0791O 
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03250  G0SU8  07390 
03260  4<K»L)*W<K»L>«45 
03270  G0SU8  98000 
03280  G0SU3  07390 
03290  4<KtL>*W<K»L>*45 

03300  next  k 

03310  KsKl 
03320  L*L1 
03330  G0SU8  07300 
03340  GOSUd  07390 
03350  «MK*L)*«MK»L)*45 
03360  S05U8  07910 
03370  G0SU3  07390 
03380  4<K»U=4(K»L)*45 
03399  G0SJ3  07620 
03400  G0Sl)3  07390 
03410  4 (*»l) «W CK,L) ♦US 
03420  *«K1 

03430  FOR  L*1  TO  Ll-1 
03440  GOS'JB  07730 
03450  G0SU8  07390 
03460  M(KtL)«M(K*L)«45 
03470  G0SU8  07820 
03480  GOSUS  07390 
03490  4(K.L)«W<K*L)**5 
03500  NEXT  L 
03510  L=0 
03520  G0SU8  07730 
03530  G0SU307390 
03540  4(<»L)*M(K*U«45 

03550  «E9 - — -SPECIAL  BASE  EN3  METAL  MASSES-—-— -v*. 

03560  30SU3  09119 
3570  REM 

03610  8EM - PRINT  ZONE  CENTERS  AND  MASSES- - — 

03620  REM— ——PUT  INPUT  PRINT  OUT  BYPASS  HERE—.., 
03630  PRINT=K=tSL=»=9l  (K.L!=»=Z1  <K*l)E*=41  tK»L)s 
03640  FOR  *=1  TO  <1-1 
03650  FO«  L*0  TO  Ll-1 

03660  PRINT  <,L*Ri <K»L> tZl «K»U t*l <K«L> 

03670  NEXT  L 
03680  NEXT  K 
03690  PRINT 

03700  PRINT  2K£*SLs«:R(K*L)s*sZ(K*L)~*sW(K*L»s 

03710  FOR  K*1  TO  K1 

03720  FOR  L*0  TO  LI 

03730  PRINT  <«LtR(XfL>  »/ (K»L) t 4 (K«L) 

03740  NEXT  L 
03750  NEXT  K 

03760  REM— ———C/M  AND  SCALE  FACTOR—————— 

03770  PRINT5C/M  and  SCALE  FACTORr 
03790  PRlNT=X*ST«OI*»ll/3)/(l*2*C/M)i 
03790  PRINT 

03800  PRINT=K=.5Z (K.O) =»=IR<CM) E.HC/M=f =Cr!I ( lN«»A/3)r 
03910  FOR  X*1 TO  Kl 

03820  Y4«R3/R4/(  (R{K»L1*1)  /RU*L1>)«»2.1> 

03930  T5*<R<X.H*i).R(XfUn»R<X*LU»6(I/3)/tl»2»r*)».36355 
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039*0  PRINT  <»ZU»0>  «R(K»L1)  »Y4»Y5 
03850  N£XT  .( 

03960  REN————. INITIALIZE  N#T— — <— — 

03870  NiO 

03889  T«(Z1  <<4*1*0)-ZK4*1,0>  >/!) 

03890  «xm— — — initialize  yi#E2»for  k«i  to  K4— — — 

03909  FOR  C*0  TO  Ll-1 
039 i 0  FOR  <»1  TO  K4 
03C20  VllKiU«l 
03930  E2(KtL)>0 
33940  NEXT  K 

03950  REM  ——————— INITIALIZE  Yl»E2tF0R  K*<4*1  TO  Kt-1  — 

03950  FOR  <*<4*1  To  <1-1 
03970  VI(K*D*I 
03980  £2 (K«L) »£1 
03909  NEXT  < 

04003  NCX7  L 

04010  REM— — — — INITIALIZE  J»v— — — — — — 

04020  FOR  <*1  TO  <1 
04030  FOR  L«0  TO  LI 

04040  ur.«u*o 

04050  Y(X«U>0 
04080  NEXT  C 
04070  NEXT  K 
4090  REM 

04093  REM  —INITIALIZE  INTERMEDIATE  SlIOE  PTS  Z4<<> ♦*(»<> 
04100  FOR  K«2  TO  K1 
04110  K9(K}*K.l 
04120  k4«)«R(i<»LI) 

04130  Z4(K)«Z(K«i.l) 


0*140  NEXT  < 

4145  01«01*30.48 

04150  R£M - —ENERGY  CNEC<— — — 

04160  GOSu’8  11400 

04170  REM—  — — FLJI3  OYNAMli'S  -  MAIN  ROUTINE 
04180  NaN«l 


04190  IF  T>T4  THEN  11970 
04200  IF  M>N6  THEN  11970 
04210  T*T*T1 

«0220  REM— - — — <3  FOR  LIMITCO  COMPUTATION 

04233  FOR  KX1-1  '0  1  STEP  -1 

042*0  iF<O*T-Z(K.6>*Z>>OTMEN0t26O 

04250  80  TO  04260 

04260  M3*K 

04270  *C*1 

J4280  NEXT  K 

0*290  <3»<3*5 

W4300  IF  <3X1*1  THEN  94320 
04310  80  TO  04330 
14320  <3«K1-1 
04330  REM 

04340  PRlNT5N«sNt5T»sTt5<3«SM3 


04350  REM— —————— NEW  POSITIONS 

04360  REM— — — r.RtZ  FOR  <3«1*I  — 


04370  FOR  <«1  TO  <3 
04380  FOR  L*0  TO  Cl 
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04390  Z<«* L)*2<K.U*T1»U{K.U 
04400  R<K»L)*R(K»L)*n*V<KtL> 

04410  NEXT  L 
04420  NEXT  K 

04430  IF  K3<«1-1  THEN  04490 

04440  REM— —————— -Q«Z  FOR  K»K1 

04450  FOR  L*0  TO  LI 

04460  Z<*ltD=2<Kl,L>*Tl*U(Kl«L» 

04470  R(K1  »U«R(X1»L)4T1»V<KI,L) 

04480  NEXT  L 
04490  REM 

04500  R£M— - - - —n£w  POSITIONS  OF  SLIDE  PTS*Z3<K>  ,R3<K>  — 

04510  G0SJ9  10200 

04520  REM  ——CELL  CENTER  VARIABLES———. 

04530  FOR  K*1  TO  *3 
04540  FOR  L*0  TO  Ll-1 
04550  P=P2 (K»L)-0I(K»L) 

04560  V5 (X*L) =V1 (K»L) 

04570  REM- - -——CELL  CORNERS  FOR  VI  (K* Ll-1)  CALC— — - 

04580  IF  LsLl-1  THEN  10660 
04590  50SU3  07300 
04600  G0S08  07390 

04610  VI (K«L1*(V34V4)4R2(K«L)/N1 (K«L) 

04620  A1 (K.L)=A3«A4 
04630  V6«(V1 1K*L)-V5(K,L) } 

04640  V7«V6/T1 

04650  Vg*(Vl (X»L)+V5(K«L) 1/2 
04660  v9*v7/v« 

04670  IF  V9>0  THEN  04700 

04690  31  (X »L 1 sC3*C3*R2 (KiL)*(A3*A4) *V9*V9/V8 


04690  GO  TO  04720 
04700  OKK.LlsO 
04710  REM 

04720  REM - —————BURN  FRACTION—————.—. 

04730  IF  MKl-i)*i  THEN  04760 
04740  IF  L=0  THEN  11740 

04750  REM - EG  OP  STATE  ROJTINE  CALC  P2*E2tC2 

04760  G0SU3  08450 

04770  REM— — .... — .— ...TIME  STEP— — — 
04780  REM - - - —TIME  STEP  FDR  EACH  CELL.T3U.L)-' 


04790  IF  L=L1-1  T1£n  04830 

04800  S3*<Zf«^l.L4l>-Z<K»L)l**2*<ROUl,L4l)-R(K#L>>*»2 
04810  S4*(Z(<*L*l)-Z(K*ltL))#»2*(R(x.L*l)-R(K*l»L) )*»2 
04920  GO  TO  04930 

04930  S3*(Z3(K4l)-Z{<,Ll)*#24(R3(K41».RtKtLn**2 
04940  S4* ( 23 (K ) -Z (K*l t L) ) **2* <R3«XJ-R 1X41 tL)l*»2 

04950  IF  S3<S4  THEN  04980 
04960  S5*S4 
04970  GO  TO  04890 
04980  S5»S3 
0**990  S5*S0R  { S5) 

04900  IF  V9>o  THEN04930 
04910  d«2*C3*S5*V9 
04920  GO  TO  04940 
04930  3*0 

04940  T3 (X»L ) *S5/3/S3R  (C2  (X *L ) *8*8 1 
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0*950  NEXT  L 
0*960  NEXT  K 

0*970  REM - - - NE*  THE  STEMS  T1  ,T2 

0*960  T6«Tt 
0*990  TS«T3(1,0) 

06000  F OR  K» 1  TO  K3 

05010  F 09  L*0  TO  Ll-1 

05020  IF  T3l*tL)-T5<0  THEN  050*0 

05030  30  TO  05050 

050*0  T5*T3(<#L> 

03050  NEXT  L 
05060  NEXT  K 

03070  IF  T5>£8*T1  THEN  05100 


05080  T1*T5 
05090  50  TO  05110 
05100  T1«E8*T6 
05110  T2*(Tl*T6>/2 

05120  REN—— ........... — q*,**  FOR  SLIOINQ— — 

05130  GO  TO  11090 

051*0  REM - VELOCITY  OF  INTERIOR  POINTS 

05150  FOR  K*2  TO  *3 
05160  FOR  L*1  TO  Ll-1 

05170  REM— ————— —CYCLONE  ACCEl  FoRHJLAS- 

03180  GO  TO  10830 
05350  NEXT  L 
05360  NEXT  K 

05370  REM— —————VELOCITY  OF  AXIS  INTERIOR  POINTS 

05380  FOR  K=2  TO  Ml 


05390  R1 IK-1 •0)*(R(K-l(0)«R(<*0)*R(K»l)»R(K-l«n )/* 

05*00  Z1 (M-l»0)B(Z(K-l*0)*Z(<t0)*Z(Ktl)*Z(M-l til)/* 

05*10  NEXT  K 
05*20  L»0 

05*30  FOR  K=2  TO  *1-1 
05**0  G0SU8  08090 
05*50  G0S03  07390 
03*60  31»V3»V* 

05*70  G0SU8  08180 
03*80  GOSUS  07390 
05*90  8l«81*V3*V* 

05500  82»Z(M*l»L)-Z(X-ltL) 

05510  G(Mt0>aU(K.0>-T2*bl*2*(P21KtL)-P2<K-l tL) )/M(KtL)/32 


05520  VHX»01*0 
05530  NEXT  < 

055*0  REM— ———VELOCITY  OF  METAL  MASS  <»0INTS— 

05550  R£M— — — — ....... u<l  1 0 1 1 V  ( 1 1 0 )  — — — 

05560  SlsRCltl 1/2 

05570  0 ( 1 »0) =U( 1 *0) -72*81*31* ( 3*P2 C 1 tO)-P2 (2t 0) 1 /*/8| 1 *01 
05580  V ( 1 tO>  *0 

05590  REM— — — — — jtV  FOR  <*1*L»1  TO  Ll-1—— 
05600  K*1 

05610  FOR  L»1  TO  Ll-1 


05620  P**<P2(KtL)*P2(«tL-l)»/2 
05630  P3»«P2{KtltL>*P2CK*l,L-l))/2 
056*0  P5«(3*P*-P3)/2 
05650  81«CR(MtL)*R(KtL"in **2/* 
OS660  82«<R(<,U*R<K.L*1>>**2/* 


i  T?^- t^v?*^  *>  T>???*^''V^'  '*^-r  s-1-;  i?  ■  f-.  %x:  -» *^'7, «;  »j!?,^T>5.-^,*^t>  rr;^’^  -,vi^  i  *>  ■Jn?i-  *■  's'-o^v^  i*».  c;-m=t^xtt'  ■«*-fTT  ^"v 
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05670  83*<3<*R<K.L>*R(K»L-l) ) • (Z (K„L$ «Z (K*L-1 > ) /8 
03680  34*(3*R(KtL>*R<K.L*l) J • (Z <K.L*1 > -Z (K.L) >/8 
05690  tMK»L)»U(K«w)-T2*<32-81)*P5/2/«ICKiL> 

05700  V(X*L>*V<KtL>*T24«33*84)*P5/M(K«L> 

05710  NEXT  L 

05720  REM— — — — — — — U»V  FOR  K*1»L*L1  — 

05730  K«1 
05740  L*l1 

05750  83»<R(<*U*R«K.L-l>)*»2/4 
05760  34* (R(K*L)*R(K»1«L) ) **2/4 
05770  85*T2/2/«MK»U> 

05780  J(K#L)*U<K«LM35»<33-B4)4P2<1*L1-1> 

05790  33*  <3*R (K»L) ♦R(K*l«L))*<Z(K*l«L)"ZCK*L))/8 
05900  94*<3*R<X*L>*R<*<*L-1)  > * < Z (XtL) -Z <K»L-1 >  >/8 
05910  V(l.Ln=V{l#Ll)*T2*<93*B4)»P2{l.Ll-l)/iCK.L) 

03920  R£M— —————— P3(l)  tP3 (K3*  1 ) -- 

03933  ?3  ( 1 )  *P2  <  1  «L1»1 ) 

05940  P3<K3»1 )*P2(K3»Ll-ll 
05950  R£M— — — — — U*V  FOR  L*Ll*K*2  TO  *3———— — ....  1 

o596o  rem. — .Extrapolation  for  P3<k>  at  slide  pts— — — —  ■ 

05970  L*L1 

05990  FOR  K*2  TO  K3 

05990  P4*(P2(K«L.l»*P2CK-ltL-l))/2  \ 

05900  P3»(P2(K.L-2)4P2«K-l.L-2))/2 

05910  P5*P4*(P4-P3)/2 

05920  P3<K)aPS 

05930  NEXT  K 

05940  FOR  K*2  TO  K3  | 

05950  J*1  4 

05960  05*SGNlZ<K,Ll)*S8U)*R(K»Ll)«S»(J»*RCJ.Ll*n-ZCJ#Ll-I)) 

05970  D6*SGN(Z(K,Ll)4S8U*l)»R<KtLl)«S8U*l)»R(J*ltLl-l>*ZU*ltLl-l)  J 
05980  IF  J5O06  THEN  06050 
05990  05*06 
06000  J*J»1 

06010  IF  J>K1-1  THEN  06030 
06020  GO  TO  05970 

06030  GO  TO  09750  : 

0&040  GO  TO  11970 

06050  07*S8R((Z(KtU)-Z3(jn#*2*(R(XtU-R3U))»*2>  ' 

06060  D8*SQR<<Z(KtL>*Z3U«m*»2MR(K*l.>*RJ(J*l>>**Z> 

06070  P5*P3 ( J) *07* (P3(J«1).P3(J))/(D7«DB) 

06080  83* (R ( K-l»L)*R(X»L))4*2 

06090  84* (R (4*1 »L> ♦R(KtL) ) *#2 

06100  85*T2/9/*(KfU 

06119  0 (K tL)*U(K»L)«95»( 33-84 )»P5 

00120  REM 

06130  REM 

06140  83*  CZ (<*L)-Z(K-l»L))/2 
06150  34* (Z<<»1»L)-Z(K»L) 1/2 
06160  85* (3*R (K»U) ♦R(K-ltL))/4 
06170  86* (3*R (X*L) ♦RIK-lfL) )/4 
06180  37*  (B3*85^ 34486)  *T2/rf (<#L> 

06199  V(K«L)*V(K*L) ♦37*P5 
06200  NEXT  K 

06210  IF  X3*4l-1  THEN  06230 
06220  GO  TO  06350 


**&^£uibjS ottss^aa 


ma mmum 
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TABLE  A-l  (CONT.) 


•U#V  FOR  K»KlfL»0* 


TO  Ll-1- 


06230  REM 

06240  REM— ................ 

062^0  «s*l 
06260  L*0 
06270  8l*R (K *  1 ) /2 
Of'ESO  REM 

06290  f»5*«3»O2{Kl-l,0l-P2Ca-2.0n/2 
06300  U(K«L)*U(K»L) ♦T2*61 *9 1*65/2/4 (<tL> 

06310  V(K»U*0 

06320  REM— — — — — FOR  K»K1  »L«1 
06330  FOR  L«1  TO  Ll-1 
06340  32*(R(<*L)'R(K«L+1) ) **2/4 
06350  93*lR«tU-  R(K«L-1)  )*«2/4 
06360  P4*(P2(K-1*L>*P2<K-1,L-1)  > /2 
06370  P3*(P2(K-2,L)*<»2<A-2tL-l>>/2 
06380  85e(3*P4-P3)/2 

06390  U (K» U>»U«K.U* (32-83) *T2»BS/6(<,L)/2 
06400  83«(3*R(K.U4R<KtL*l))»(2(K,U)-2(K,L*l))/8 
06410  34*<3*R(K»L) ♦R(KtL-l) ) • (Z <K t L-l ) -Z <K *L ) ) /8 
06420  V (K»L)«V(K«l) « (33*84) *85»T2/rf (<»L) 

06430  NEXT  L 

06440  REM— —————— Ufy  FOR  K*K1  •  L*L1  —————————— 

06450  L«Ll 

06460  33*  IR(*'-l*L)*R(KtL))  #*2/4 
06470  B4b(R(K«L>*R(K»L-1)  >##2/4 
06430  B5BT2/2/M(KtL> 

06490  J  (M« U  *'J  (K.L)*  (33-94)  #*2 (K-i *L-1 )  *85 
06500  33B(3«R(K*L)«R(K-l*L))*(Z<KtU-Z(K-l»L)  >/8 
06510  34*(3#R{K,L)*R<«*L-l>)»(Zl«»L-l>-Z(KtL)>/8 
06520  7  CKtL)  *V (K tL; * (33*84) #P2 (K-l *L-1 ) *T2s B (K*L) 

06530  REM 
06540  REM 

0or-50  REM— — — — r(.0B  VARIABLE  PRINT  ROUTINE——— 

06560  IF  N*1  TMENU6590 

06570  IF  T>«T7#N?*THEN  06590 

06580  SO  TO  07130 

06590  IF  N<N4  THEN  06610 

6600  SOSUB  20020 

06610  PRINT  =N*£N»=T«5Tf £T1*5TJ *iT2«rT2 
06620  PRINT 

06630  REM  FOR  F(.o«l  VARIABLE  PRINT  BYPASS 

06640  REM— —so  TO  END  OF  FLO*)  VARIABLE  ROUTINE— 

06650  REM 

06660  FOR  K*1  TO  K3 

06670  PRINTs5*sKtcF (X) *sF (<) «iMV£LBsSQR(U(<»Ll > ##2*V (Kt LI) *62) *32808 
06680  PRlNT5L5»5Z5f5RS»sJ5,rVS  ° 

06690  FOR  L«0  TO  LI 

06700  PRINT  L*/(^«D »R(K»L) »U(K»L) *V(K*L) 

06710  NEXT  L 
06720  PRINT 

06730  PRlNTSLr»sPl*<31£tr31s»5VlrtSE2s 
06740  FOR  L*0  TO  Ll-1 

06750  PRINT  L»P2(M.UtOUK»L)*Vl(K»L)»E2(K,U 
06760  NEXT  L 

06770  PRINT5«=#sLs»Sr3<KfUSt=C2{K#L)S 
06780  FOR  L*OTO  Ll-l 
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0&790  PRlNTK*LtT3<K*L) *C2(<«L1 
06900  NEXTL 
06310  NEXT  K 

06320  IF  K3<>KI-1  THEN  06920 
06330  K=Ki 
06940  PRINT=<s=K 
00350  PRlNT=Li»iZi*5REt=J=*SV= 

00969  FOR  L»0  TO  LI 

06970  PRINT  L*Z<K*L> *RIK»L) • J<X*L> *V(K*L) 

06980  NEXT  L 
06890  PRINT 
06900  PRINT 
06910  PRINT 

06920  PRlNT=KS*=«3=*sZ3=*=P3= 

06930  FOR  K*1  TO  K3*l 

Q6940  PRINT  K*R3<K) *Z3(M *P3<K> 

06950  NEXT  K 
06960  PRINT 

06970  PRlNT5<=»iLs»rVEL < FT/SEC) H 

06980  Kal 

00990  60SU8  11920 

07000  IF  K3<K1-1  THEN  07030 

07010  KsKl 

07020  G0SU8  11920 

07030  REH 

07040  PRINT  =Ki«=K9 (<) =*iW4 (K) =»=Z4 {<) =.=S8 (K)= 

07050  FOR  K»1 TO  K3*l 

97060  PRINT  K,K9(K> *R4<K1 »Z4<K)  .SlMK) 

07070  NEXT  K 

07080  REH - ———BO  OF  FLO*  VARIABLE  PRINT  ROUTINE- 

07090  IF  N»1  THEN  07130 
07100  N7»N7 ♦ 1 

07110  IF  N7>=N5  THEN  11970 

07120  REH - - - - - - -ENER6V  CHECK . . 

07130  G0SU3  1 14J0 

07140  REH - - - RETURN  TO  St  ART  OF  MAIN  ROUTINE - - 

07150  REH 
07163  GO  TO  04180 
C  7 1 7 0  REH 
07180  REH 

7199  REM  — - - - INITIALIZE  HE  AND  FUZE  CELL  DENSITIES' 

07200  FOR  L»0  TO  t-1 
07210  FOR  K»l  TO  K4 
07?20  R2<K«L)»R2 
07230  NCXT  K 
07?40  FOR  K»K4*  1  7!)  <1 
07250  R2(K»LJ*R3 
07260  NEXT  K 
07J70  NEXT  L 
0'280  RETURN 

07290  REH - - - - -SUBROUTINES  FOR  MASS* VOL JME* AREA 

0  73 jO  61JZ(K*L> 

07310  H1«R(K*L) 

07320  G?«Z <K*l *L) 

07330  H2sR<KM*L; 

07340  G3»Z (K*l *L*i ) 
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TABLE  A-l  (CONT.) 


% 


% 


If 


i 


1 


07350  H3»R<K*l»L»lJ 
07360  G*«Z(6.L*1> 

07370  H**R (K*L*1 ) 

07380  RETURN 
07390  X1»G1 
07*00  V 1 *H1 
07*10  X2*G2 
07*20  Y2»H2 

07*30  IF  K<Kl/2  THEN  07*70 

07**0  X3«G3 

07*50  V3*rt3 

07*60  QO  TO  07*90 

07*70  X3*G* 

07*80  Y3«H* 

07*90  G0SU8  07650 
07500  V3«V2 
07510  *3sA 
07520  83»8 

07530  IF  K<Kl/2  THEN  07570 
075*0  X2*G* 

07550  Y2«H* 

07560  GO  TO  07590 
07570  XJ»G3 
07580  Y1 «H3 
07590  G0SU8  07650 
07600  V*aV2 
07610  A**A 
07620  **»* 

07630  *5«*3 ♦•»* 

076*0  RETURN 
07650  Afc  »X1 •Y2-X2*Y1 
07660  A7«X2*Y3-X3»Y2 
07670  A8oX3*Y1-K1*Y3 
07680  AsA8S(A6«A7«A8)/2 
07690  R»«Yl*Y2*Y3)/3 
07700  V2«A*R 
07710  8aV2«H2(K,U) 

07720  RETURN 
07730  Gl»Z<K#U 
077*0  H1«R(K*U 
07750  G2*Z(K«1'L) 

07760  H2«R(64*»U 
07770  G3«(ZO'«l«LM>*G2>/2 
07780  H3*(RrK4l»L*l)*H2)/2 
07790  G*a(Z(KtL»l)«Gl)/2 
07900  H*s(RtK»L*lJ*Ml)/2 
07910  RETURN 
07*20  Gl«Z(KtU 
07930  H1«R(K»L) 

079*0  G2>Z(K«lfU 
07950  H2»R (K*l » L) 

07960  G3*  <Z  <K»1 »L"1 ) »G2) /2 
07970  H3*(R(K»l»t“l) *Hi) /2 
07980  G*a(Z(X«L-l)«Gl)/2 
07990  H*a<R<K»L-ll*Hl)/2 
07900  RETURN 
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TABLE  A-i  (CONT.) 


07910  31*ZK,U 
07920  H1«R<K,U 
07930  S2jZ<KiU*1) 

0  79*0  H2*R(K,L»1) 

U7950  33*  < Z  <K-1 »L*1 ) *62) /2 
07960  H3*<R(K-l»L*l)*H2)/2 
07970  G4» (Z (K-l »L>  ^Gl ) /2 
07980  H4*<R<K-l*U*Hl>/2 
07990  RETURN 
08000  Gl«Z(K«U 
0B010  Nl»R(l\*U 
06020  G2*Z ( K t L* 1 ) 

08030  H2=R(K*L*1> 

0»040  33*(Z(K»l*L4l)*G2>/2 
06050  H3«<R<K*l*L*lJ*H2>/2 
06060  G4*(ZtK»l»L)*Gl)/2 
0»070  H4*(R<K*1»L)  ♦Hll/2 
06080  RETURN 
08090  G1*Z(X,U 
08100  H1»R(X»|.) 

06)10  G2*(ZK*2«U*Gl>/2 
08120  H2«(R(K-l.L>»rtl>/2 
08130  G3«Z1(<«L> 

081*0  H3=R1  U.L) 

08150  G4*(Z(«.L*l}*Gl)/2 
08160  H4*(R(K?i.»l)*Hl)/2 
Obi  70  RETURN 
06)50  G1»Z iK»L) 

08190  H1*R<K.L> 

08200  GP* (Z(X*L*l)*Gl>/2 
06210  Hpfc (R(KtL+l)*Hl)/2 
08220  G3*Z1 <K-i»L> 

08230  H3*R1 (K-l.L) 

08240  G**(Z<K-l»l!*Gl)/2 
08250  H4»(R(K-1»U)*H1)/2 
08260  RETURN 
08270  G1«Z(K«U 
06280  H1*R(K»L) 

06290  Gc*(Z(K«L-l)*el)/2 
08300  M2«<R(X,L-l,4hl)/2 
0»310  G3*Z1 (K-i»L-l) 

08320  H3*R1 K*l«L-l> 

08330  G4*<Z(K-1.U *Gl>/2 
08340  R4*<Rt«-l#L>*Hl >/2 
06350  RETURN 
08360  31  *Z<K«U 

06370  hi*r;k»l) 

08330  G2»(Z<K»L-lWQl)/2 
08390  H2e(R(X,L-l)*Hl)/2 
08400  33«Z1 (KfL-1) 

08*10  H3»R1 (K«LW3 1 
08*20  54» (Z CK*1 *L1 *G1 1 
08430  H4*(R{K*ltL)*HI)/2 
08440  RETURN 

08450  RE«  — — *•—  EO  Of  STATE  ROUTINE  CALC  »2tE2*C2- 

08*60  IE  N>0  THEN  08490 
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08470  305UB  08760 
08480  00  TO  08750 
08490  IF  K<*K4  THEN  10160 
06500  N2*0 
08510  HsVlIK.U 
08520  C4*02/H 

08530  C5*D3*  ( ’.-C4/F1 )  <*£XP  I-FlSH)  *1)4*  (1-C4/F2)  *£XP  <-F2*6> 
08540  E6*£2 (K»L) 

0»550  F0«  N2«o  TO  1 
0856.0  Pl*C4«t2<K.L)*C5 
C85T3  *P  (K) 

08580  P3»Pl*antSfC< 

08590  £2<K«L>»E6-: f P.P1 5 /?*ca <K»U  )<*V6 
06610  NEXT  N2 
08620  P2(;;«u«p3 
08630  F3«-02/H/H 

08640  F4»-03»«F3/Kl4Fl»Jl-C4/Fin*EXP(-Fl«H) 

08650  F4«F*-i;6*<F3/F2*F2»<l-C4/F?>>t>EXP<-F2*H> 

08661  C£<K.L)*-IH*H/'»JMK«U  }«{-PlfcC4*E2(X»L>*F3*F4) 

08670  IF  C2(K.L)>0  THEN  08750 

08680  PRINT=N*=N'5K«SK*=t.«rU=®I*SPl»sVl  (K.L)sEVl  <K»U 
08690  PRlNTi£2(K*L)a=£2(<fk.»  «I£6«sE6 
08700  PRINTER <K*L* *59 <K*U  »=*  (<»U  *SZ  (K.LI 
06710  PRINTER <K«1 »L> *SR I **1 #l> »SZ <*♦! »L>*SZ  OU1 »L> 

08720  POINTS*  <**1  »LM  )*ERU*'..M.*I>*5Z<K*i»LM  I  *SZ«KM#U*i> 
08730  PRINTER <X*L* 1 J *5* <K»l* 11 »=Z  (K*L«1  )■=£  (K.LM  ) 

08740  PRlNT=R3fK)«=R3(K) t=Z3(X)«=Z3(K) 

08750  RETURN 
08760  Os. 7839 
08770  02**30 
0B780  R3*l .634 
08790  03*4.9055 
0B89C  04*. 058 
08810  El*. 0814 
08820  Fls4.2 
08930  F2*0.9 
08840  83*1»0. 731304 
08950  PRINTS J6L  EONs 

08860  PR’NTSO*=0.=RHD(0>*sR3.=02==02»=03s=03.=04s=04 
08870  PRlNTE£l*EElE(HB«R«CC/CC)i*EFlsEFl*EF2sEF2 
08880  PRINT  EQ3*1-V1 IC J) *533 
08890  PRINT 
08900  RETURN 

08910  REM— ——BASE  CORNER  MASS  CM4N3E-— ------ 

08920  X1*(Z(K1-1»L1)«Z(K1«U)  1/2 
08930  Yl*<R<Kl-WLl)*R<Xl»Ll»>/2 
08940  X2*X1 

08950  V2*(R(Kl-l.Ll*l)4HIKltUlfi)>/2 
08960  X3* 1X1 (Kl+ltLl+l) 1/2 
08970  Y3«(Y2^RIKlM»Ll»l))/2 
08980  80SUB  07650 
08990  6(KI-1.L1I»«(K1-1»UI)*6 
09000  tf(Kl.U>**<Kl»Ll>-6 
09010  Xl*(Z(Kl#LlUZ(KltLl-l>)/2 
09020  Yl*(R(Kl*Ll)«R(Kl»Ul*l))/2 
09030  X2*(Z(KlM*il)*Z{Kl»!»Ll«l>)/2 
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TABLE  A-l  (CONT.) 

(•■#043  Y2* ( W (Kl ♦ 1  #L1  1 ♦R<Kl*l /2 
(>9o50  X3»{X<J*Z(M*l,LUl>)/2 

09(160  Y3*(Y2*«(K1»1  *U>1) )  /2 

ovnro  Gosud  o7=50 
UYOttO  *UM.U-na**Kl,LW)*rf 
UY090  rt(Kl«Ll 
09100  RETURN 

09110  RE*-- - - — 3  A  St  EMU  MASStS-105  *M. 

09120  K*K1 

09130  L*1 

09140  G0SU9  07920 

09150  GOSUet  07390 

09lt>0  *MK»L)  *W5 

091  70  G0SU8  (*•'  730 

09180  33=36.9 

09190  H3»2.37 

09?o0  GOi-Ua  07390 

093J0  *  <K»L)*W0V«L1'*W5 

09220  L=? 

09330  G1*Z(K»L) 

09240  81=R(K«U 

0*250  62*36.9 

092f,o  H2*3.2b 

09270  GOSUa  07390 

093*10  W(K«U=W5 

09290  64=<ZK,U*Z(K.L*l>)/2 

09300  H4*(Rt<*Ll»RtK*L*l))/2 

09310  63=36.9 

09320  H3*4.30 

09330  GOSJH  07390 

09340  *<KtU«*<!\»L)4WS 

09350  L*3 

C93h0  G1*Z(K,L) 

09370  HI*R(h,l) 

09380  62=36.9 
09390  H2=4.55 
09400  GOSJS  07390 
09410  «MK«t)  *m5 
09420  64=33,35 
09430  *4*2.8 
09440  63*34,55 
09450  H3=4,95 
094b0  60SU3  07.390 
09470  rf<5»U*W<5»U**5 
09480  K=*l-1 
09490  U*L1 
09500  G1*Z(K«L) 

09510  H1=R (K*L> 

09520  32*32,9 
09530  *2*5.2 
09540  faOSUa  0739y 
09050  «M**U*rfS 
095&0  G4*32.2 
0*570  H4=3.15 
09580  33*32.5 
09590  *3*5.25 
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0*600  GOSU9  073*0 
0*61  0 

0*620  KsKl-2 
0*630  L*Ll 
0*640  Gl *Z (K • L ) 

0*650  H! «R (K  »L 1 
0*660  G2*Z(K»L*l* 

0*670  H2*H(H.LM> 

0*660  G0SU3  073*0 
0*6*0 

0*700  sosua  07*10 
0*710  aosua  o7i*o 

0*720  « (fttL) >M (K*L> *«5 
0*730  RtTuRN 

0*740  REM— EXIT  IF  J>K1*1  1*  MASS  PT  PRESSURE  ROUTINE— 
0*750  PRlNT=aORY  pressjke  ROUTINES 

0*760  PRINT=JsrJ*sHo=K»rR<<»Li)«s«<K»tl> ♦=4(<»U)«SZtK»Ul) 
0*770  PRINT  =J=»5Sft< J1=»=H< J«L1-1 »=.=Z{J*L1-1>=«5065 
0*780  POR  Jal  TO  M 

0*7*0  36»Z(K,L11*S8(  J)*H<K.L1>-SB«Jj*P<  J*L  1-1  > -Z  UtLl-U 
0*900  WRINTJ.SBIJ)  .M(J#U«1>  »Z(J*Ll-l>  »06 
0*810  NEXT  », 

0*920  GO  TO  11*70 

0*930  RE* - *E*P  P*E  I T£PA V ION t At.T  METMOOI— — « 

0*940  F5«F(«) 

0*950  E6s£2<<.D-lP*31  <6tL)l*V5 
0*960  P3*(C4*E6-*C5)«F5 
0*970  E2(K.L)*Eb-,5«(P3-P)»V6 
0*9«0  Pl«(C4«E2(K*L)*C5)*F5 
0*9*0  P2(KfUePl*«l  (<»U 
0**00  GO  TO  08630 

0*510  REM-—— ———INPUT  RK»L1>  •9t<»Ll*l)  STEP  i- 
0**20  FOR  *30  TO  <1*1 
0**30  :«£A0  «(KfL11  .R(KtLlM) 

0**40  NEXT  K 
0*550  RETURN 

0**60  RE**  — — ——INPUT  R(k*LD  *R(K*L1»1>  STEP  ?•••■ 
0**70  READ  R(0»L11  .R«0*Ll«l> 

0**80  FOR  Kal  TO  K1  STEP  2 
0**90  READ  R(K«L1>  »R<KtU*)> 

1UOOO  NEXT  k 

10010  READ  R(K1*1.L1>»RH1*1»L1*1> 

10020  FOR  9*2  TO  *1-1  STEP? 

10030  R<K»Ll  >*<R(K*t*Ll  >  *R(X*1»U)  >/2 
1  U 040  R(A.Ll4l)»{R1K*l»Ll*n*R(K-l*Ll*l>)/2 
1  0  050  NEXT  K 
10060  RETURN 

10070  REM— .-.——INPUT  SPECIAL  *  <*tLl »  *R  1  VALUES— 

10080  R«20*L1 1*3.09 

1U0*0  R(20«L1*D«S».2 

U'100  R(2l«Ll )«2.55 

10110  R(21*U*I)«S,05 

10120  Rc22.U)=R(21.l:! 

10; 30  RETURN 

10140  REM— — EQ  OF  STATE  FOR  FUZE  CELLS—"* - - 

lolbG  P 3*P2(<»L) 
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TABLE  A-t  (COMT.) 

10160  »2!««D>»3t/tVl 

10170  £2(K,L)sE2(K*L)*((P3*P2(K«tl)/2)4>V6 
10180  C2<K#L»*#07/R2<KtL>/»Vl(KtL>##6> 

10190  80  TO  08670 

loroo  rex— ..............  ........  .slide  routine— 

10210  L*L1-1 
10220  IF  N»1  THEN  10310 
10230  FOR  K.l  TO  K3«i 
10240  R3<K}«R(KtLll 
10250  23 (K) *Z (K  ,Ll > 

10260  NEXT  K 
10270  00  TO  10680 
10280  REX 
10290  REX 
10300  REX 

10310  FOR  K*i  TO  *3 
10320  REX 
10330  «EX 
103*0  REX 

10350  S8(K)«»(24(K).2(K«L) )/(R*lK)*R(KtL) 1 
10360  REX 
10370  J*1 

10389  05*Z ( JtLl  J  ♦  IRl-JiLl  )»R  (K*L>  )«S8|K)»Z  (K«Ll 
10399  06*2 C J«2 VL1  > •(R(J*I til )*R<K*U)*S8(K>«ZfK*L> 

10*00  IF  06*0  THEN  11540 

10*10  IF  SON (06 ) «»56N (05)  7 HEN  10*60 

10420 

10430  IF  J*Kl  THEN  10520 
10**0  05*06 
10450  60  TO  10390 

10*60  S7*(R(J«l«U)«R(JtUn/<Z<j«ltLl)*Z(JtU>> 

10*70  23(K‘aZ<*tLJ  »S8  HO  »S7«Z  ( Jtf.1 )  -S8  (K)  •  (R I  J,L1  >  -R(K#t  11 

10*80  23<F)*23lK»/|ltS8lK»«S7l 

10*9o  R3(K>*R?wtU)-$7*Z(J#Ll)*S7«Z3(K> 

10500  K9 (K ) *j 

10519  60  70  10570 

10520  RRIHTSSLIOE  ROUTINE, J*=Jt=K*=K 
10530  80  TO  11970 

105*9  RJ (F 1 *R (W*i *Lt ) 

10950  73(F)«Z(\;*ltU> 

10569  K9<K>ij 

10570  REX 
10520  NEXT  K 
10590  R3<l)*R(ltLl> 

10600  73(l)*Z(ltLll 
10610  »3(K3tl)*R(R34ltLl> 

10620  Z3(K3«l)*2(K3ti,t.l) 

10630  S8<1)*1 

106*0  S8(Klt«*l 
1065C  RETURN 

1 0660  REX  ———CELL  CORNERS  FOR  L«H-l— — 

10679  61*2 (KtL) 

10680  Hl*R(KtL> 

10690  02*Z (K«l »L  J 
10700  H?»R (xtl *L! 

10710  63*Z3 (K*l ) 


NOLTR  74-77 


TABLE  A-!  (CONT.) 


10720  H3c*3u»U 
10730  04>Z3(<) 

10740  M4aR3,<) 

107S0  90  TO  04600 
10760  REM 
10770  REM 
10780  REM 
10790  REM 
10800  REM 
10810  REM 

10820  REM - - - CYCLONE  ACCEL  FORMUL4S— — 

10830  Hla(P2(K»L)*R2<KiL-l)‘»R2(K-l*D-R2!K-l*L*'U>/2 

10840  M2* (P2 (Kt L)  +R2<K~1  tL) *?2 (K«l *L**1 ) -P2 (Kif'ln/2 

10850  81*(R(K*1»L)-R<K-1#L) 1/2 

10860  IF  L<M  -f  THEN  10900 

10870  92* (R*  <K.L-l)>/2 

10380  84*'23..  .  «L-l)>/2 

10890  00  TO  109c* 

10900  82*<R(K»L*1)-R(K*L«1 J  > /2 
10910  84*(Z(<tL«l)-Z(K«L-ll)/2 
10920  83*(Z(K»ltL>>Z(K«icL) *  V2 

10939  94*(Z(KtL*!)-Z(K»L-ll  )/2 

10940  JS*T2*R<K*L)/M<K«U 

11 000  U(K«L)*U(K«L) *35* (91*H2»82*H1 } 

11010  rf(«tL)*V(K,L>*85*(64*M1.93«M2j 
11020  90  TO  05350 

11030  REM— ——INITIALIZE  SLIOE  PTS,ZJ(K)  »R3(K)  — — 
11040  FOR  K*1  TC  Kl 

11059  R3(K)a')(r.»LD 

11060  Z3(K)a2<ft*Ll> 

11070  NEXT  K 

hi 080  RETURN 

11090  REM— —••-CALC  R4*Z4  FOR  SLI0IN9-— 

11100  IF  K9U3>*0  THEN  11120 
11110  9C  TO  11130 
11120  K9(K3|aK3-l 
11130  FOR  Ka2  TO  «3 
11140  R8*R(K9(K1 ,L1) 

11150  Z6aZ(K9(K)«Ll) 

11160  R60R(K9IKJ *1,L1> 

11170  Z9*Z (M9(M)*l»Ll) 

11180  09*$QR( ?Z9-Z8) **2* (R9-R8 J  **2 J 
11190  $9«(R9-R8)/09 
11200  C9a(Z9**Z8)  /09 
11210  S5*P2(<tLl-i>-P2(K-l,Ll-l) 

11220  G6»Z3(K)»Z(K«L1«1) 

11230  97*R3'.<0«R(K»Ll«l) 

11240  98*R3 (8) /2/W3 ( <9L1 ) 

11250  G9*G8*35*<S6*S9-87*C9) 

11260  90  TO  11340 

11270  35**2 1  <U1*L1~1)* 1*3 <K«1)»R(K»L1*1)) 

11280  96*P2<<*U»D*<R<KtLl»ll-R3<*Ul>  > 

11290  S7*R2U-l*Ll-ii«MZ3<X-l>*Z<KtLl-I>> 

11300  S8*P2<<.U«  l )«MZ(KfLl-l)-Z3 !*♦!>> 

11310  REM 

11320  F9»P3«MAl  <K»1 1 Ll*l  >/Vi  (K»ltLl*l)«Al  («»11<-1) /Vl  (KtLl-1)  )/4 
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11339  S9*S9/2/F9 

11340  J7<K>=J(*»Ll-l)*T2*C9*39 
11353  y7<*)*V<*»Ll-l)*T2*S9*39 
1 1 350  Z4U)*Z3;k)4J7(K)*T1 
11370  R4<K)s93(KJ*v7<k)#T1 

1 13»9  ME* T  K 
1  i  390  SO  TO  05140 

11490  a£N  ————————— ENERGY  CHECK-"— 

11410  £3*0 
11420  £4*0 
11430  £5*0 

f 1 440  FO«  K*1  TO  Kl-1 
11450  F 09  L«0  TO  Ll-1 

11450  E3«£3*€2(K.L)«2*«»9»(<!  <<«L> /32 (<*L> 

11470  MEAT  L 

11499  NEXT  K 

11490  FOB  <*2  TO  K3 
11599  EOa  L*4  TO  Ll-1 

11510  E4*£4^09«i»(K.U»<U(X.L1»*2*V«K»L)**2) 

11320  NEXT  L 
11530  NEXT  « 

11540  «*1 

1 1559  F09  L*9  TO  Ll-1 
11550  S05JB  11710 
11570  NEXT  l 
11590  L*Ll 
11590  FO*  K*1  TO  K3 

11500  £4*£44!>99*3<<.L»*<jf  (K»*»2*v7«<>»«2) 

11510  E5»E5439*n(K.l> •<U(X»U»»2*v(K.L)**2l 
11520  NEXT  K 

11539  IF  K3«l-1  THEN  11530 

11540  K«K1 

11559  FOB  L»f!  TO  LI 
115*0  SOSJB  11 710 
11570  NEXT  L 
11590  E«£3«E4*£5 

11590  PSlNT5N*sN*5IE3AS*sE3*S<E<JAS*i£4tsK£NET«s£f 

1 1 ?  9  3  aETuHN 

11710  £4*£4*P9*H3  (<».)  •«J(K»L)  **2*  V  C<»L) 4#2> 

11729  £5*£5*»9*«<X»L)4<'HK*L)»*2*y (K*L»##2l 
11739  SET jaw 

U740  hem— ——309.9  fraction  iflane  detonation) - 

11759  IF  F (K ) 3 1  THEM  1.900 
11759  IF  0*T<(Z(5»9)-ZJ7HEN  11999 
11770  IF  K*K4*1  THEw  11370 

11730  F<X)«O*WUt0)*Z)/(ZU4l»Ll»-ZIK*0J ) 

11790  F(<)*F(K)/1.5 

11500  04* ( 1 -V 1 (K«0) )/J3 
11910  IF  M5XJ4  TMEV  11930 
11929  C*0  TO  11940 
11919  F  xi  -05 

11440  %'  F (n; >1  THEN  1 1970 
l 1 950  IF  FlKXO  THEN  11H90 
11550  30  TO  11900 
1  1579  FX»*1 
11430  30  !0  11900 
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11990  F{K)*0 
11900  R£M 
11910  50  TO  04760 

11920  REM— — CALC  MASS  ?T  VEL  £0R  K'  *  ..'«Kl  •  — 
11930  FO*  L*0  TO  Ll-1 

1194C  *»RlNT  <«LtSa3<  J(<»U  »*2*V  1K»L>  ••21«32»09 
11950  NEXT  l 
11960  RETURN 
11970  60  TO  30000 


20000  R£M«— — -Fpa&mEvT  PREDICTION  SCmEME— — 
20010  R£M 

20020  9£M— ...............  I6MT  PER  POLAR  ZONE 

20060  RE"«— -———SET  M7  (  J2*M51  «0 


20050  EOR  J2*l  TO  37 
20060  FOR  M5*l  TO  2*L1«*1 
20070  47  ( J2«  **5)  so 
20090  NEXT  M5 

20099  NEXT  j 2 

20100  R£M  — — — — ANQtE  A9<MS*l/2>  — 
20110  PRINT5M55  =6=»”Lf«=45=»'«9(M5*l/2)= 
20120  R£M 

20130  R£M 
20140  M£M 
20150  REM 
20160  R£M 

20170  FOR  MS* 1  TO  tl 
20190  K*1 
20190  L«MS-1 

20200  R5*lRlM.C)«R(K.L*in/2 
20210  Z5*(2(<tLl*ZJM.L4l) 1/2 

20229  J5*U(<.L«15 

20230  V5*0U,L*ll 

2o240  031  M3) *S3R{U(K»LJ #*24V;<»L> **2> 
20250  SOSU9  2U90 
2026C  NEXT  MS 

20270  FOR  M5*L1M  TO  Li^M-l 
20260  L*L 1 
20240  «*M5-L1 

20300  R5*1R«<<U»*R(K*1»L)) /2 
20310  Z53<Z(<«LWK*ltU  >/2 
20320  U5*U(9M»U 
20330  V5*V (X* i »L  > 

20  34  0  33  (*5>  «S  J»  <U  IK »L>**2*V  »<,d  **2) 

2>1350  GOSU9  21340 
2.1363  NEXT  MS 

20370  *OR  ms«L1*<1  TO  l1*2*«1-1 

20390  <**1 

2u390  Us2*Ll*Xl-M5 

20409  R5*<R<X»U*R<*»L-1  J 1/2 

20410  2S*UK«cW(K»l-li5/2 
20420  U5«J(*.L-1> 

20430  W5*V1«,l-H 

2j440  *3<M5l  *S3R  (UttWU  #*2*V  l<,L>  **21 
20450  30SJ9  21340 
20460  N£XT  xj 
20479  49  f  0 >  *0 
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20*03  46(2*U««l)»tf(<lfO> 

20*90  Q3(2*LWKl>»SQR< J (<1 *0) »*2*V (Kl *0) »»2) 

20500  A9(2*LW<1)*180 

20510  REM— —SET  N7»«5tM6t96»N3»N8»ZER0— — — 

20520  FOR  J6*l  TO  37 

20530  F0»  J*1  TO  J9 

205*0  N7<J6«J)«tf 

20550  R5(J6.J)*0 

20560  M6(J6«J)*0 

20570  06(J6»J)>0 

20580  N3 { J6«  J>  *0 

20590  ME X T  J 

20600  N8U6>«0 

20610  ME*T  J6 

20620  PRINT 

20630  PRINT=MASS  PT  =  »=P0|.A9  Z0NE=tiwT  IN  ZONE  (GRAINS)  5 

206*0  REM— ——TOTAL  FRAG  «T  AND  AVE.  VEL  IN  POLAR  ZONES 

20650  FOR  M5sl  TO  2*L1*M 

20660  IF  A91M5-l)>A9«M5)  THEN  20700 

20670  C7sA9<M5*l) 

20600  C0*A9(M5I 
20690  SO  TO  20720 
20700  C7«A9(M5) 

2u710  C8=A9(M5-ll 
20720  J2»INT (C7/5) *1 
20730  J3alNT(C0/5)*l 
2u7*0  IF  J2«J3  THEN  20030 
2o  759  J*»J3-J2-1 
20760  86«C3-C7 

20770  M7(J2.M51»(J2*S-C7)«26(M51/00 
20700  m7U3«M5>«(C0-(J3*IJ*51*#6{M0)/06 
20790  FOR  N9*l  TO  J* 

20500  *7U2«N9«M5)s5**6(MS)/06 
20910  NEXT  N9 
20920  SO  TO  200*0 
20930  *7<J2,M5)*«6<M5> 

209*0  FOR  Jb» J2  TO  J3 
20950  FOR  j*i  TO  J9 

20960  N7U6  J)*N7<J6«JMN(M5»J1»07(J6»M5>*P0 
20970  M5  C J6» J) **5 ( J6t J) *22 (M»t J) *M7 ( J6«M5) 

20980  NEXT  J 

20990  PRINT  M5«J6»0?(J6*NS>#39 
20900  NEXT  J6 
20910  NEXT  H5 

20920  REM  ———————TOTAL  WEIGHT— — — — — 

20930  T8»0 

209*0  FOR  M5«l  TO  2»Ll*Xl 
2o950  T3»T9*m6(M5) 

20969  NEXT  M5 

20970  T9«T8»»8 
20980  BRINT 

2«990  PRlNTiTOT  METAL  *T«=T9=L3 l =T9*7000=GRAINS) 5 
21000  PRINT 

21010  PRINT s*0LAR  ZONEi«5MET  WT5»=PCT  BY  WT=*=AV  VELs 
21023  FOR  J2*l  TO  36 
21030  *8(J2)«0 
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21040  89(J2)*0 

21050  FQR  «5»1  TO  2*tl*Kl 

21060  *8<J2)««SU2)»*7  <J2#M5> 

210  70  89  ( J2 )  xw9 ( U2 )  ♦  *7  ( U2 • M5 )  ( M5 ) 

21060  IF  J2<>36  THEN  21110 
21099  «8(U2)**8(J2)*87(37.M5) 

21  ICO  *9  ( J2)  stf9 ( J2)  *87  (37tM5)«03(H5) 

21 1 1 C  NEXT  M5 

21120  r(U2>*89<J2>»32808/89< J2> 

21130  PRINT  J2 *88 ( J2 1 *i9» «9  <  J2 ) #99/T9/70 • Y  < J2) 

21140  NEXT  J2 

21150  R£M— ...... ......... ...mo.  OF  FRAGMENTS  IN  POLAR  ZONES 

21  loO  G0SJ8  21710 
21170  G0SU3  22270 
21160  G0SU3  22480 
21190  RETURN 


21195  R£M— —.—END  OF  FRAGMENT  PREDICTION  SUBROUTINE 
21200  REM - — - - - - —HEIGHT  DISTRIBUTION  INPUT  — 


21210  FO«  J=1  TO  J9 
21220  READ  m< J) 

21230  NEXT  J 

21240  FOR  M5«l TO  2«l1*<1 
21250  READ  Mi <M5) 

21260  IF  Ml (MSI <*50  THEN  21290 
21270  J6*  20* (//1S)*(M1 <M5>-5?> 

21280  GO  TO  21300 

21290  U6*20*0.^366*<M1 (MS). 50) 

21300  FOR  J*1  TO  J9 
21310  XI«M(J)/J6 
21320  X2«SQRU1) 

21330  X3«EXP(-X2) 

21340  N (M$« j) s7* 70925*X3/U6 
21350  «2(M5.J)*X3*(0.5«X1*X2*1) 

21360  NEXT  J 
21370  NEXT  M5 
21360  RETURN 

21390  REM— -——ANGLE  A9(M5*l/2) 

21400  81*N(K«l1*U5 
21410  32*U(KiL)MV5 
21420  33*U(K.L)#U5 
21430  34*V<K.L>*V5 

21440  35»SaR (3l»81 ♦32*32»63*93*84*d4 ) -83*84 

21450  IF  ABS <851 <E7  THEN  21600 

21460  A4*(31«32)/85 

21470  A5»A4* (R0-R5) *25-20 

21480  *6*A4* A4 

21490  A7«A5«A5 

21500  R6*S0R { A6*A7- i A6+ 1 ) • ( A7-0 1 *01 ) ) 

21510  R6«R0*(R6-A4*A5)/<A6*1) 

21520  26*2 5* A4# 1R6-R5) 

21522  IF  M5<ll*l  THEN  21525 

21523  IF  M5X1*H  THEN  21525 

21524  GO  TO  21330 

21525  IF  (Z6-Z0)MZ5«Z0)<0  THEN  21600 
21530  IF  A3S<Z8-Z0)<E7  THEN  21580 
21540  A8*ATnuR6-«0)  /(Z6-Z0)  M180/P9 
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21550  IF  Z(>-Z0<0  THEM  21650 
21560 

21570  30  TO  21660 

21580  A9<M5>«90 

21590  SO  TO  21660 

21600  IF  Z5-Z0<0  THEM  21630 

21610  A9(MS)*10O 

21620  60  TO  21660 

21630  A9 (MS) »0 

21640  SO  TO  21660 

21650  A9(MS)«.A0 

21660  RE* 

21670  46 (M5) *4 (K»L) 

21680  PRINT  H5*K,U*46(M5) *A9(H5) 

21690  RETURN 
21700  REH 

21710  REM- - - - —MO.  OF  FRAGMENTS  IN  POLAR  ZONES 

21720  PRINTsM5=tsM9ARs 
21730  FOR  M5«l  TO  2*L1>*1 
21740  PRINT  M5«Ml(M5> 

21750  NEXT  MS 

21760  PRINT 

21770  FOR  J6»l  TO  36 

21780  FOR  JalTO  J9 

21790  IF  J*J9  THEN  21850 

21900  N3( J6t J)*N7 ( J6* J)-M7 (J6. J*l) 

21916  Xl«45(J6tJ)-H5(J6t J*l) 

21920  IF  J6«>36  THEN  21840 

21930  M3(J6t J)*.M3U6.J)*M7(37,J)«N7(37tJ*I) 

21940  GO  TO  21890 
21950  M3(J6. J)*N7(J6#J) 

21960  Xl«45(J6tJ) 

21970  IF  J6<>36  THEN  21890 
21980  N3(J6.Jl»N3(J6.J)*M7(37t J> 

21890  M6(J6«J)«X1*P6/N3(J6*U) >15.4324 
21900  N8 ( J6) sN8 < J6) >N3(JbtJ) 

21910  NEXT  J 

21920  Y9< J6)sN81J6)/Y8<JSJ 

21930  NEXT  J6 

21940  PRINT 

21950  FOR  J6«l TO  36 

21960  FOR  JslTO  J9 

21970  96(J6tU)sN3(J6tJ)/M8(J6) 

21990  NEXT  J 
21990  NEXT  J6 

22000  PRINTS  NO.  OF  FRAGMENTS  IN  POLAR  ZONESs 

22010  PRINT24T.  RANGES  IN  GRAINSs 

22020  PRINT 

22030  U«1 

22040  G0SU9  22140 

22050  PRINT 

22060  J*5 

22070  GOSU8  22140 

22090  PRINT 

22090  PRINTS  POLAR  ZONES *M (J9-1 )  =-sM ( J9) isGT,  THAN2M(J9) 
22100  FOR  U6«l  TO  36 
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22110  PHInT  J6*N3 ( J6* JV-1 ) .N3( J6.J9) 

22120  NEXT  JO 
22130  RE f URN 

2t  >>0  PE'  — — PRINT  NO.  OP  ."RAGS  IN  POLAR  ZONES*  — .. 

22150  MIIMIJMI 
22160  X (2, «M{ j«2) 

22170  X(?)«0*(JO) 

22160  X(4)aM>'J«4) 

22190  PRINTSPZS»H(J>S-SX(1>  * X ( 1 > E*~X (2) *X (2) S*iX (3)  * X  ( 2;  £-=X  (4 ) 

22200  POR  J6»l  TO  36 

22210  POR  M2*J  TO  JO 

22220  N3  < J6.M2) *N7 { J6*N2) *N7 ( J6«*2*D 

22230  NEXT  M2 

22240  PRINT  J6t N3 ' J6« J) » N3 ( J6« J»1 ) .N3 ( J6t J«2) *M3(J6» J«2) 

22250  NEXT  J6 
22260  RETURN 

22270  REM— — ................  pm  *(T  LETHAL  AREA  INPUT—————— 

2228C  PRINT 

22290  PRINTS  LETHAL  AREA  INPUTS 
22300  PRlNTSHaAV.  NT  IN  3R4IN53 

22310  PRINTS  0*PRACTI0N  OP  NO.  IN  ZONE  >  i  GRAINS 
22320  POP  J6*l  TO  36  STEP  2 

22330  PRINT  ! J6-1 ) •5:-=J6«5* =-=» =-S»  J6*5S»i ( J6« 1 ) *5 
22340  PRINT  SMi«SUS.S*i»sMS*sOi 
22350  POR  J»1  TO  J9 

22360  PRINT  N6(J6*J) «Q&( J6.J) «S*£*H6(J6*1»J> i06( 

22370  NEXT  J 
22380  PRINT 
22390  NEXT  J6 
22400  PRINT 
22410  PRINT 

22420  PRINT  SAV  NT »A V  NT  In  GRAINS  OP  PRAGNENTS>1  GRAINS 

22430  PRINT  S  POLAR  ZONEStSlNIT.  VEL. (PPS)i» SpRAG/STERS.SAN  NTIGRATNSs 

22440  FOR  J6*l  TO  36 

22450  PRINT(  J6-U*5S-=  J6*5 •  V  ( J6)  t T9  ( J6)  tM6 ( J6) /NB  (U6)  *Q9 
22460  NEXT  J6 
22470  RETURN 

22480  REH  — — RARE  JnEN  DECK  FOR  LETHAL  AREA  INPUT  .............. 

22490  Ti (O)sASS(U(ltO) >*32808 

22500  VI (36) aA8S (U(Kl*0) >*32808 

22510  FOR  J6*l  TO  35 

22520  VI U6»*(V (J6)*V( J6*l) ) /2 

22530  NEXT  J6 

22540  RESTORE  PILE(TAPEI) 

22550  POR  J6»  1  TO  35 

22560  PRINT  PILE  (TAPED  ( J6-D  *5»  J6*5#  C  J5-.5)  *5» Y1  (U6-D  *T1  ( J61  »Y  ( Jfc)  »J9 

22570  POR  J»l TQJ9 

22580  PRINT  PILE CTAPE1 )M6 ( J6» J) I 

22590  NEXT  J 

22500  PRINT  PILE(TAPEl) 

22610  POR  J«1  TO  J9 

22520  PRINT  PILE (TAPEl ) N3 ( J6. J) I 

22530  NEXT  J 

22640  PRINT  PlLE(TAPEl) 

22650  NEXT  J6 
22660  RETURN 
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22670  DATA  3*3,05*2.35*3.83*3,2*4,65 

22680  DATA  3. 85* 5.07,4,2.5.3.4.23,5.37 

2269Q  DATA  4. 07,5*26*3. 86*5. 27, 3. 65, 5. 27, 3, 45, 5. 27 

22700  DATA  3.27.5.4*2,97.5.17,2.97.4,57 

22710  DATA  1,2*5*10*15*25*50*100*150*250 

22720  DATA  200 *200 .200 *  103 *  105, 100 .75.55 ,40*30 

22730  OATA  35*35,45*50*70*95*110*130*145*160 

22740  DATA  190*65*65*20*35*45*45 

30000  ENO 
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Statement 

Number 


1000 


1010-1050 


1060-1100 


1110-1300 


1310 


1320-1370 


1380-1500 


1510-1780 


m 


NOLTR  74-77 


Table  All 


Notes  on  Fragment  Prediction  Code  List 


Notes 


INPUT  AND  INITIALIZATION 
Identification 


Fixed  constants,  it,  2i r,  2*  *  15.4185,  5if/l80,  2ir/454. 


Calculate  and  store  the  number  of  steradians  in  each 


5°  polar  zone,  for  use  in  the  fragment  prediction 


scheme  (see  Fig. 9).  The  area  of  the  curved  surface 
of  a  spherical  segment  of  height  h,  of  unit  radius, 
is  2irh,  or  2ir(l-cosa),  where  a  is  the  polar  angle. 

The  surface  area,  or  number  of  steradians,  between  the 
polar  angles  ax  and  a2  is  then  2ir(coso2-cosuj ) . 


These  are  constants  for  the  computation,  defined  in 
Table  AI. 


Equation  of  state  constants  are  gotten  from  the  equation 
of  state  subroutine  (8450-8900). 


Constants  used  in  the  computation  are  printed  out. 


Dimension  statements  for  subscripted  variables.  The 
dimensions  and  descriptions  are  listed  in  Table  AIII. 

GRID  POINT  INPUT  AND  INITIALIZATION 
Dimensions  in  cm,  taken  from  drawings  of  the  weapon 
casing  and  fuze,  are  inserted  here.  At  the  same  time, 
the  values  of  the  interior  grid  points  are  obtained, 
either  by  uniform  spacing  or  by  reading  in  special 
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values.  The  example  in  the  program  list  (Table  AIV)  Is 
for  the  105  mm  projectile,  which  has  a  curved  base 
(see  Pig. 3)*  Special  grid  point  entries  were  made  in 
order  to  get  more  uniform  cell  sines  in  the  base  region, 
and  to  meet  the  requirement  that  each  cell  have  four 
corners.  These  special  entries  are  not  needed  if  the 
base  walls  are  essentially  perpendicular  to  the  axis, 
as  in  the  155  mm  and  5"/54  cal.  projectiles.  Input  and 
initialisation  for  the  example,  the  105  mm  projectile, 
were  done  as  follows: 

Dimensions,  in  cm,  were  obtained  from  a  drawing  of  the 
weapon  casing.  The  fuze  was  replaced  with  a  steel  plug 
in  the  nose.  Also  plane  detonation  of  the  HE,  starting 
at  the  nose  was  assumed. 

1510-1550:  Enter  metal  boundary  points  on  the  axis 
(see  Fig. 1) . 

1560-1590:  Initialize  h«0  values  by  putting  uniform 
spacing  between  Z, „  and  2„  _ . 

1600-1610 :  Input  special  coordinates  for  the  axis  points 
I63O:  Read  in  rh,«  .R*,.**,  (O*  Ll+l)  in  cm,  taken  from 
the  drawing,  assuming  HK  4  =  2,^  “  Z»(0  .  This  can  be  done, 
one  K  line  at  a  time,  in  9910.  To  reduce  the  number  of 
data  entries  one  can  read  in  every  other  R  value  between 
K=1  and  K=K1,  and  have  the  program  produce  the  inter¬ 
mediate  points.  This  is  done  in  9960.  In  the  example 
(Table  AI)  9960  was  used  and  the  data  appear  in  22670- 
22700. 

1650 :  Special  R«iW  values,  replacing  those 

already  inserted,  are  put  in  here,  to  make  the  cells  more 
uniform  in  the  base  region.  This  was  done  in  10070. 
1660-17^0:  Grid  point  coordinates  are  calculated,  using 
uniform  spacing. 
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Number 

1790 

1820-1870 

1890 

1910 

1930 

1970-2050 


Table  All  -  continued 
_  Notes 


1750-1780 :  Special  values  are  Inserted  here  to 

replace  those  Just  calculated  for  the  line  k»ki 
FRAGMENT  WEIGHT  DISTRIBUTION  INPUT 
Fragment  weight  distribution  input.  In  the  example 
(Table  AI),  the  Mott  distribution  is  used  and  the 
weight  distribution  is  calculated,  from  Eqs.  (82)  and 
(83),  in  21200. 

CASINO  DIMENSION  HUNTOUT 

Print  out  coordinates  hKc  for  K*0  to  Kl+1,  3>0  to 
Ll+1.  Notice  that  this  includes  the  outer  boundary  of 
the  metal,  which  will  be  used  only  to  calculate  the 
masses  associated  wiuh  the  mass  points.  To  bypass  this 
print  out  insert  1805  GO  TO  1875. 

INITIALIZE  Z3H,  R3„ 

Initial  values  of  the  slide  points  Z3„,R3,,  are  taken  as 
,  respectively.  This  is  done  in  11040-11080. 
SET  Z  FOR  LIGHTING 

Set  Z  3  •  This  18  the  value  of  H^where  the  HE 

column  starts.  In  the  burn  fraction  calculation  for 
the  plane  detonation  case,  the  distance  that  the 
detonation  front  has  traveled  Is  measured  from  this 
line. 

INITIAL  CELL  DENSITIES 

The  initial  densities,  of  the  fuze  material  or  the 
undetonated  HE,  are  stored  for  the  individual  cells 

This  refers  to  the  cell  with  the  lower  left 
corner  at  (k, l). 

COORDINATES  OF  ZONE  CENTERS:  ZONE  HASSES 
Calculate  coordinates  of  zone  centers  by  averaging  the 
coordinates  of  the  corners.  Calculate  the  scaled 
masses, wt  associated  with  the  cells.  This  is  done 
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Number 


2070-2200 


2220-2300 


2320-2360 


2380-2870 
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with  7300-7380,  which  locates  the  cell  corners,  (Fig.Al), 
and  7390-7720,  the  subroutine  which  cuts  the  cell  into 
two  triangles  and  calculates  their  areas,  volumes,  and 
scaled  masses.  In  2030, W5  is  the  sum  of  the  scaled 
masses  of  the  two  triangles,  calculated  in  7390-7720. 

INTERIOR  GRID  POINT  MASSES 

Get  scaled  masses  associated  with  the  interior  grid 
points.  The  centers  of  the  adjacent  cells  are  located  in 
2070-2160,  and  the  triangle  subroutine  7390-7720  is 
used. 

INTERIOR  AXIS  POINT  MASSES 

Scaled  masses  associated  with  the  axis  grid  points 
are  found  by  calculating  the  scaled  masses  of  the  two 
quarter  cells  at  the  points.  Subroutines  starting  in 
8090  and  8180  (see  Fig.Al)  are  used  to  locate  the 
corners  of  the  quarter  cells. 

MASS  CORRECTION  AT  THE  FUZE-HE  BOUNDARY 
In  the  calculation  of  the  scaled  masses  associated  with 
the  interior  grid  points  (2070-2200),  the  densities 
of  the  cells  with  (K,L)  at  the  lower  left  were  used. 

For  grid  points  on  the  line  K=K4+1  (K4>0,0<I£L1-1), 
this  is  now  corrected  to  take  into  account  different 
initial  densities  on  both  sides  of  the  line.  Here  R2 
is  the  fuze  material  density  and  R3  is  the  HE  density. 

It  is  assumed  that  the  cells  on  both  sides  of  the  line 
K*K4+1  are  of  the  same  length  in  the  if  direction.  Other¬ 
wise,  these  masses  should  be  gotten  by  calculating  the 
masses  of  the  quarter  cells  on  both  sides  of  the  line. 

Scaled  gas  masses  associated  with  the  metal  mass 

points  are  calculated  (see  Fig. 5).  If  the  boundary 
point  Is  not  a  corner  point,  these  are  the  sums  of  the 
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scaled  masses  of  the  two  adjacent  quarter  cells.  If 
the  point  is  a  corner  point,  (1,0),  (1,L1),  (K1SL1), 
or  ( K1 , 0 ) ,  W3,ik  is  the  scaled  mass  associated  with  the 
adjoining  quarter  cell.  The  subroutines  in  7290-8M0 
are  used. 

2880-3600  The  scaled  masses  associated  with  the  metal  mass  points 

are  calculated  here,  with  the  subroutines  in  7290-8M0. 
The  metal  is  split  up  as  in  Fig. 5.  In  the  example  used 
here,  the  105  mn.  projectile,  special  coding  for  the 
base  end  was  done,  so  that  the  lines  dividing  the 
masses  would  be  perpendicular  to  the  lines  Joining 
them,  as  shown  in  Fig. 3.  This  special  coding,  done 
In  9110-9730,  is  called  with  3560-GOSUB  9110.  An 
alternate  arrangement,  useful  if  the  base  is  perpendi¬ 
cular  to  the  axis,  is  shown  in  Fig. A3.  This  Is  done 
with  the  statement  3560-GOSUB  8910. 

PRINTOUT  OF  GRID  POINT  AND  ZONE  CENTER 
COORDINATES,  AND  SCALED  MASSES 

3610-3750  Printed  out  here  are  the  coordinates  of  the  zone 

centers,  21*.,.  ,  the  scaled  zone  masses 

the  grid  point  coordinates  and  the  scaled 

masses  associated  with  the  grid  points  and  metal  mass 
points  .  This  printout  can  be  bypassed  by  inserting 
the  statement  3625  -  GO  TO  3760. 

C/M  RATIO  AND  SCALE  FACTOR 

3770-3850  The  C/M  ratio  ar.d  ~X  ,  defined  in  Eq.  (80)  are  calculated 

and  printed  out  here,  for  the  metal  at  the  cross  sections 
where  K*1  to  Kl.  These  are  used  to  get  from  Fig. 10  or 
a  sxm'.lar  relationship,  the  values  of  R  to  assign  to 
the  mass  points  on  the  side  wall  of  the  projectile.  It 
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is  usually  convenient  to  stop  the  program  at  this  point, 
with  the  statement  3855  -  GO  TO  30GG0.  After*  the  values 
of  R  for  the  various  mass  points,  determined  from  the 
calculated  X  values,  are  put  in  as  data  for  the  fragment 
prediction  scheme,  which  starts  at  20000,  the  statement 
3855  is  removed  and  the  entire  program  is  run  from  the 
beginning. 

INITIALIZE  N,  T,  AND  FLOW  VARIABLES 

3870  Set  the  cycle  number  N»0. 

3880  Set  the  initial  time  equal  to  that  for  the  detonation 

front  to  traverse  half  of  the  first  HE  cell. 

3900-4000  Initialise  and£2K>u.  Usually,  to  start,  vl„t=  1  , 

E2  K  =  o  in  the  fuze  material  region,  and  E  2^  -  £  1  (=  p„  El) 
for  the  HE  cells. 

4020-4070  Initialize  the  velocity  components  u„>t  and  Vw><> 

4080-4140  Initialize  the  intermediate  slide  points  H4h>R4„  , 

and  the  marker  points  K9„. 

4145  Convert  dimensions  of  Dl,  the  arena  radius,  from  ft 

to  cm. 

ENERGY  CHECK  FOR  CYCLE  ZERO 

4160  The  total  energy  in  the  system,  at  cycle  zero,  before 

the  computation  starts.  Is  found  and  printed  out  here 
with  the  energy  check  subroutine  that  is  used  later  at 
the  end  of  eacn  cycle. 


Statement 

Number 


V 


4i8C 

4190 

4200 

4210 

4220-4330 


4340 

4360-4480 

4500 

4530 


FLUID  DYNAMICS  -  MAIN  ROUTINE 


Advance  tine  cycle  number. 

Exit  If  time  exceeds  T4  microseconds. 

Exit  if  cycle  number  exceeds  N6 
Advance  time  (see  Eq.  (17)). 

KJ  for  limited  computation  is  calculated  (Eq.  (18)). 
This  is  done  to  save  computing  time,  by  not  calculating 
over  portions  of  the  grid  where  nothing  is  happening. 
For  the  cycle,  K3  is  the  maximum  K  for  which  new  values 
of  the  variables  located  at  the  cell  centers  are 
calculated.  If  K3=K1-1,  the  motions  of  the  end  points 
(K1,L)  are  also  calculated. 

The  new  cycle  number  N  and  the  new  time  T*  are  printed. 


NEW  POSITIONS 

Equations  (1$^  are  used  here. 

New  slide  point  positions  i3HRZ>K  and  related  variables 
are  calculated  Ir.  10200-10650 

NEW  VALUES  OF  VARIABLES  JCATED  AT  THE  CELL  CENTERS 
This  is  the  start  of  a  double  loop  which  ends  at  4960. 
New  values  of  ,  qC’*,  .  £2^  ,  C2;v 

are  calculated  In  this  loop. 
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45 50 


4560 


ii  570-H620 


4630 


4640 


4650 


4660 


4670-4700 


4?30 


4740 


4750 
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Saves  old  Pifck(*  PZ^-Oi^J  temporarily. 


Saves  old  vi*  temporarily. 


New  cell  areas  and  relative  specific  volumes, Ai.>t  and 
Vl"tJare  calculated  with  the  same  triangle  subroutine 
used  for  the  initialization  (7300).  If  L  <L1~1, 
subroutine  7300  is  used  to  get  the  cell  corners.  If 
L*L1-1,  1C660-10750  is  used.  This  takes  into  account 
the  fact  that  the  upper  corners  of  the  cells,  for  this 
case,  are  (2f3^.R3w)  and  (Z3„^  ,  . 


=  Vi..,  -  vi,; 


.  "**/t 

v7  =  (vi)  =rdvi/dT) 


ve  *  vi 


* (Vi/v) 


Ql„  is  calculated  with  Eqs.  (29)  and  (30). 

BURN  FRACTION 

If  all  F*  now  equal  one,  no  burn  fractions  need  be 
calculated. 


This  is  the  plane  detonation  case.  For  each  X  the  burn 
fraction  is  the  same  for  all  L.  This  i?  calculated 
when  L*0 ,  in  11740-11910. 


P 2.c  ,  E2"t  and  C2*w  are  calculated  in  the  aquation  of 
state  subroutine  (8450-8900). 
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i*770-^940 

4950 

4990-5060 

5070-5110 

5130 

5150 

5380-5530 
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The  time  step  for  the  (K,L)  cell,  T3*\  is  calculated 
Kith  Eqs.  (33) -(35).  Notice  that  iZ3fc,R3m)  and 
(Z3„.t.«3w  )  are  the  upper  corners  of  the  cells  when 
L-Ll-1. 

Save  old  time  step  temporarily. 

T5  is  the  minimum  of  the  cell  time  steps  T3*t 

The  new  time  steps  Tl"*1  and  T2****  ere  found  with 
Eqs.  (56)  and  (27). 

H4H}ReHF0R  SLIDING 

Intermediate  slide  point  positions  are  found 

ir.  11090-11390,  while  the  old  velocities 
are  still  available. 

VELOCITY  0?  INTERIOR  GRIP  POINTS 
New  interior  grid  point  velocities 
are  calculated  in  10830-11020,  with  the  scheme  in 
Eqs.  (5CM54). 

VELOCITIES  OF  INTERIOR  AXIS  POINTS 
The  new  velocity  components  of  the  interior  axis  points 
are  gotten  with  Eqs.  (56)- (58).  Subroutines  used  In 
the  initialization  are  used  here  to  calculate  the 
volumes  of  the  two  quarter  cells  associated  with  the 
axis  grid  points.  Extrapolation  of  the  pressure  to  the 
axis,  in  the  radial  direction,  is  not  done  because 
dv/dT*oon  the  axis.  This  implies,  from  (9),  that 
3P2/dR*o  on  the  axis.  Hence,  the  pressure  can  be 
assumed  uniform  between  R*Oand  R*iX/2.. 

VELOCITIES  OF  METAL  MASS  POINTS 
u.  Is  the  axial  velocity  component  of  the  axis  mass 
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5600-5810 


5830-5840 

5870-5930 

5940 

5950-6040 
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point  on  the  left  boundary.  The  pressure  associated  with 
the  point  is  taken,  by  extrapolation,  to  be -(3P2“„ - P2^#)/2. 
The  difference  form  of  Eq.(10)  is  used,  with  Xt**T(R,A/2f  • 


Velocity  components  of  the  mass  points  on  the  line 
K=1  (L^0)  are  found  by  using  Eqs.  (64)  and  (65),  with 
approprjate  subscripts.  For  example,  for  K*l,  L=L1 

A,  *  *■(  R*  -  R?),  where  R,  Rl*(Rtu  +  RM».,V2. 

and 

A,'2ir[(^u-ZtiU.i)/2][(3Rtill  +  Rl>tM)/4] 

For  the  corner  mass  point  (1,L1)  the  corner  cell  pressure 
is  used,  i.e.,  in  Eq.  (10)  PS  =  -Pz  ^  is  used  in 
calculating  the  axial  component,  and  P2lilt.j  is  used 
in  calculating  the  radial  component.  The  opposite 
signs  are  needed  so  that  the  accelerations  will  have 
the  correct  signs. 

Assign  pressures  in  the  cells  (1,L1-1)  and  (K3,L1-1) 
to  the  slide  points  (Z3,  ,rs4)  and 
respectively . 

Get  pressures  P 3.  at  the  slide  points  Z3H>  P3m  f 
(i$K«K3)by  extrapolation  from  the  interior,  with 
Eq.  (59). 

Start  of  loop  for  velocity  components  of  mass  points 
(KjLl)  for  K  €  K3. 

Finds  closest  slide  points  on  both  sides  of  the  mass 
point,  with  Eq.  (6l).  If  there  is  no  sign  change. 
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6050-6070 

608C-6190 

6200 

6210-6220 

6240-6520 

6560 

6570 

6580 
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diagnostics  are  printed  out  In  9750-9810  and  the 
program  3tops. 

The  pressure  at  the  mass  point  is  found  by  interpolating 
between  adjacent  slide  points,  using  Eqs.  (62)  and  (63). 

Calculate  mass  point  velocity  components,  for  2«tc<K3, 
with  Eqs.  (64)  and  (65). 

End  of  loop  which  starts  at  5940. 

If  K3 <ki-i  velocity  calculations  for  the  mass  points  on 
K*K1  are  bypassed. 

Velocity  component  calculations  for  the  mass  points  on 
the  right  boundary , K*K1 .  The  formulas,  with  appropriate 
subscripts,  are  the  same  as  those  used  for  the  mass 
points  on  the  left  boundary ,K*1. 

FLUID  DYNAMICS  -  PRINT  ROUTINE 
Values  of  the  flow  variables  will  be  printed  at  the 
end  of  the  first  cycle,  before  the  second  cycle  is 
calculated,  unless  this  instruction  is  deleted. 

The  fluid  dynamics  print  routine  and  the  fragment 
prediction  routine  are  entered  at  Intervals  of  T7  usee 
(actually,  at  the  first  cycles  past  these  times).  This 
is  done  N5  times.  The  counter  N7  is  advanced  by  one  in 
7100  each  time  the  fragment  prediction  scheme  is  used. 

Print  routine  and  fragment  prediction  bypass. 


Statement 

Number 


6590 


6600 


6610 
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If  the  cycle  number  N  is  less  than  N^  the  fragment 
prediction  routine  is  bypassed.  Usually  the  fragment 
prediction  scheme  is  not  used  until  all  the  mass  points 
have  started  to  move. 

The  fragment  prediction  is  dene  here,  starting  in  20020. 

Prints  N  and  t" ,  TlN*,/*'>  tz”  in  ^isec.  Notice  that  N 
and  T  are  for  the  cycle  just  completed.  The  time  step 
Tl  is  for  use  in  the  next  cycle,  and  Tz" was 

used  to  calculate  the  new  velocity  components. 


6660 


The  flow  variable  printout  can  be  bypassed  on  every 
cycle,  if  one  is  just  interested  in  the  fragment 
prediction  results,  by  Inserting  the  statement 
6650  GO  TO  7080. 

Start  of  flow  variable  printout  loop.  BASIC  programming 
now  allows  for  a  maximum  of  five  numbers  per  line. 


6670  Prints  K,  the  burn  fraction  (this  is  the  plane 

detonation  case),  and  the  metal  velocity,  i.e., 

(U* *  v*)1  *  for  the  mass  point  at  (K,L1),  in  ft/sec. 

6680-6710  Prints  L,  2  and  R  coordinates,  and  velocity  components. 

Units  are  cm  and  cm/jjsec. 

6730-6760  Prints  L,  P2r#<.  (mbar),  (mbar),  VlHu[(cc/gram)*p0] 

and  E2"(k  [  (r.bar-cc/gram)  *  p0  ] 


6770-6800 


Prints  k,  L,T3„iU  (usee),  and  [ (cm/usec) 2 ] , 


End  of  loop. 
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6830-6880 

6920-6950 

6970-7030 

7040-7070 


7090 

7100 
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If  KS^Ki-l  the  printout  for  the  line  K*K1  is  bypassed. 

Prints  L,  Z  and  R  coordinates,  and  velocity  components 
for  the  right  end  mass  points. 

Prints  positions  of  slide  points  Z3*,R3r.  and  associated 
pressures  P3« .  Recall  that  the  adjoining  cell  values 
of  P 2Hl  are  used  for  the  corner  points  (5830-5840)  and 
the  other  values  are  gotten  by  extrapolating  t>2,,Avaluus 
in  (5870-5930). 

Prints  velocities  (u*  +  vM  *  of  end  mass  points  in 
ft/sec.  The  velocities  are  calculated  in  11920-11960. 

Prints  out  values  of  variables  H4..R4*  ,  SB*  , 
used  in  the  slide  routine.  Note  that  Z4„ , R4„  are 
calculated  for  use  during  the  next  cycle.  The  variable 
Ka„  denotes  the  K  value  of  first  mass  point  to  the 
left  of  the  Kth  slide  point.  If  (?3*,R3,,)  and 
(Z*^,  coincide,  K4*=k-1  The  variable  K4*  Is 

used  in  11090-11390  to  locate  the  line  segment  along 
which  the  point  (zs*,^*)  moves  to  (Z4M,R4H).  In 
another  version  of  the  program  which  allows  for  free 
gas  expansion  out  the  ends  of  a  tube,  it  is  used  to 
tell  if  a  slide  point  is  Inside  or  outside  the  tube. 

If  it  is  outside  K4,=  Oor  Kl. 

Does  energy  check  if  N=1  and  returns  to  the  beginning 
of  the  mair.  routine. 

Advances  N7  (see  statement  6570). 
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7110 

7130 


7160 

7190-7280 

7290 

7300-7380 

7390-7720 


Table  All  -  continued 
Notes 


Exit  if  T  >  T7-  N5. 

Does  energy  check  (lr.  -.'5-11730)  every  cycle. 
Provision  to  stop  the  program  if  the  total  energy 
varies  by  more  than  a  prescribed  amount  can  be  inserted 
at  this  point.  For  example,  put 
A095  E9=E 

7132  If  ABS  ( E-E9 )  <  1 . 1  E9  then  71*»0 

7133  Print  "STOP  ON  ENERGY  CHECK" 

7135  GO  TO  30000 

Return  to  start  of  main  routine. 


SUBROUTINES  AND  TRANSFERS 

INITIAL  DENSITIES  OF  CELLS 
Initialize  HE  and  fuze  cell  densities  R2,iL.  It  is 
assumed  here  that  the  HE  column  starts  at  K^K^+l. 
Special  initial  densities  can  be  assigned  to  individual 
cells  in  7271-7279. 

MASS,  VOLUME,  AMD  AREA 

These  routines  are  used  in  the  initialization  and  in 
the  calculation  of  vi"t  and  (z«k(H3). 

Locates  corners  of  cell  with  (K,L)  at  lower  left 
(see  Fig.Al). 

Calculates  area,  scaled  volume  and  scaled  mass  of 
quadrilateral  with  corners  (G<,  ,  t=  1  to  4  /  by 

dividing  it  into  two  triangles  and  using  Eqs.  (1)  and 
(2).  At  large  expansions  the  corner  cells  (l,ul-l)  and 
(K1-1,L1-1)  tend  to  become  concave  if  the  masses 
associated  with  the  corner  points  are  large  relative 
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7730-7810 


7320-7900 


7910-7990 


8000-8080 


8090-8170 


8160-8260 


8270-  j50 


8360-8440 


8460-8480 


Table  All  -  continued 
Notes 


to  the  masses  of  the  adjacent  points.  The  subdivision 
into  triangles  is  done  differently  for  the  left  and 
right  sides  (see  Fig.A2),  in  order  to  get  the  true 
areas  of  these  corner  cells. 

Locates  corners  of  half  cell  with  (K,L)  at  lower  left 
and  (K+1,L)  at  lower  right  (Fig.Al). 

Locates  corners  of  half  cell  with  (K,L)  at  upper  left 
and  ('K+1,L)  at  upper  right  (Fig.Al). 

Locates  corners  of  half  cell  with  (X,L)  at  lower  right 
and  (K,L+1)  at  upper  right  (Fig.Al). 

Locates  corners  of  half  cell  with  (K,L)  at  lower  left 
and  (K,L+1)  at  upper  left  (Fig.Al). 

Locates  corners  of  quarter  cell  with  (K,L)  at  lower  left 
(Fig.Al). 

Locates  corners  of  quarter  cell  with  (K,L)  at  lower 
right  (Fig.Al). 

Locates  corners  of  quarter  cell  with  (K,L)  at  upper 
right  (Fig.Al). 

Locates  corners  of  quarter  cell  with  (K»L)  at  upper 
left  (Fig.Al). 

EQUA1T0N  OF  STATE  SUBROUTINE 
If  N*0,  the  equation  of  state  constants  are  entered 
and  printed  out  with  8760-8900.  Notice  t/at  this 
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subroutine  must  be  called  before  the  cell  and  grid 
point  masses  are  computed,  because  it  contains  the 
value  of  R3,  the  solid  HE  density.  (This  is  done  in 
1310.) 

8490  If  1«K<K4  then  P2*4.>E2"ai  C2*l  are  calculated  for  the 

inert  fuze  material  between  the  lines  K“1  and  KcK4+l. 

This  Is  done  in  10150-10190. 

8500  Set  N2*0.  This  is  the  counter  for  the  iteration  which 

solves  for  the  pressure  and  energy. 

8510-8530  The  equation  of  state  used  here  Is  Eq.  (11),  with  C4 

and  C5  from  Eq.  (13). 

8540  Saves  the  old  energy  E2^  temporarily,  as  E6. 

8550-8620  Solution  of  Eqs.  (6)-(8),  by  iteration.  The  iteration 

Is  done  only  once.  The  iteration  can  be  done  by  the 
method  in  the  HEMP  Code2,  by  inserting  the  statement 
8535  CrO  TO  9830. 

863O-866O  Calculates  the  sound  speed  squared,  C2M#1.  using  Eq.  (12). 

8680-8740  Diagnostics  if  the  calculated  value  of  C2^  is 

negative.  This  will  cause  the  program  to  stop  when 
it  attempts  to  take  a  negative  square  root  In  the  time 
step  calculation  (4940). 

8760-8840  Put  detonation  product  equation  of  state  constants 

here.  The  constant  Q3  (=  1 -vicy)  where  vit,  is  the 
relative  specific  volume  at  the  Chapman  Jouguet  state 
is  used  in  the  burn  fraction  calculation  (Eq.  (25)). 
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8850-8880 


8910-9100 


9100-9730 


9750-9820 


9830-9900 


9910-9950 


Printout  of  equation  of  state  constants  the  first  time 
the  routine  is  called  (1310). 

SPECIAL  CODING  FOR  BASE  END  MASSES 
This  routine  reduces  the  mass  of  the  mass  point  at  the 
base  corner  (K1,L1)  and  adds  the  mass  removed  to  the 
adjoining  mass  points  (K1-1,L1)  and  (Kl,Ll-l).  The 
metal  then  associated  with  the  three  mass  points  is 
shown  in  Fig. A3.  To  use  this  routine  insert  3560 
GOSUB  8910. 

Special  coding  for  the  masses  associated  with  the  mass 
points  can  be  done  here.  In  this  particular  case,  the 
105  mm  projectile  base  (Fig. 3'  is  split  up  so  that  the 
separation  lines  are  perpendicular  to  the  lines  joining 
the  mass  points,  and  the  associated  masses  are 
calculated.  The  routine  is  called  in  3560,  with 
GOSUB  9100. 

MASS  FOINT  PRESSURE  -  DIAGNOSTIC  PRINTOUT 
Diagnostic  printout  if  J>K1-1  in  the  mass  point  pressur'* 
routine,  i.e.,  if  the  two  adjacent  slide  points  are  not 
found.  It  is  called  in  6030. 

ALTERNATE  PRESSURE  ENERGY  ITERATION 
Alternate  method,  used  in  HEMP2  for  the  pressure,  energy 
iteration  in  the  equation  of  state  routine  (8450-8900). 
Insert  with  8535  GO  TO  9830. 

INPUT  METAL  BOUNDARY  POINTS 
Input  of  metal  boundary  points,  one  K  line  at  a  time. 
Called  with  1630  GOSUB  9910. 


9960-10060  Input  of  metal  boundary  points,  every  other  K  line 

between  K*1  and  K*K1.  The  Intermediate  points  are  put 
in  automatically,  with  10020-10050.  Called  with  1630 
GOSUB  9960. 
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10070-10130 


10140-10190 


10200-10650 


Input  special  values  for  metal  boundary  grid  points. 
This  is  needed,  for  example,  when  the  base  is  curved, 
as  in  the  105  mm  projectile.  Called  with  1650 
OOSUB  10070. 

EQUATION  OF  STATE  OF  FUZE  CELLS 
Equation  of  state  for  cells  with  1*K<K4.  Calculate 
P2*  E 2*1 1  C2*  here.  Called  with  8490. 

SLIDE  ROUTINE  (S3*  ,  R5„  .  Kfl„,SSL) 

The  new  slide  points  Z3a,R3,  and  the  related  values 
K4*  and  S8h  are  found  here  (see  Eqs.  (20)- (24)).  The 
values  of  S8,.  and  58,^  in  10630  and  10640  are 
fictitious,  inserted  for  use  in  the  mass  point  pressure 
calculation  (5960,5970). 


10660-10750  Locates  cell  corners  when  the  upper  corners  are  slide 
points  (L*L1-1) .  Called  in  4580. 


IIO30-IIO8O  Initialize  slide  points  H3„ji?3h  by  setting  them  equal 
to  >  respectively.  Called  in  1890. 

11090-11390  Intermediate  slide  points  Z4KlR4„  (see  Fig. 6)  are 
calculated  with  Eqs.  ( 38 )— ( 49 ) .  Notice  that  the 
pressures  in  the  (K,L1-1)  cells  are  used  to  move  the 
two  quarter  cells  associated  with  a  slide  point.  This 
is  done  because  the  assumption  that  the  boundary  is 
being  held  fixed  while  the  gas  slides  implies  that  the 
normal  pressure  gradient  is  zero  there.  Called  in  5130. 

ENERGY  CHECK 

11410-11730  The  gas  internal  energy,  the  kinetic  energies  of  the 

gas  and  the  metal  mass  points,  and  the  total  energy  in 
the  system  are  summed  with  Eqs.  (66)- (63)  and  printed 
out.  Here 
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11740-11910 


11920-11960 


11970 


Notes 


E3  *  internal  energy  in  the  gas  (mbar-cc), 

E4  =  kinetic  energy  in  the  gas  (mbar-cc), 

E5  ■  kinetic  energy  in  the  metal  (mbar-cc), 

E  *  total  energy  in  the  system  «•  E3+E4+E-5. 
Called  before  the  first  cycle, in  4090  (this  is  the 
total  energy  released  by  the  HE),and  at  the  end  of 
every  cycle,  in  7130. 

BURN  FRACTION  (PLANE  DETONATION) 

The  burn  fraction  FK  is  calculated  from  Eq.  (25). 
Called  in  4740. 

Velocities  of  the  end  mass  points,L*0  to  Ll-l,are 
calculated  from  the  components  uKL  and  and 

printed  out  in  ft/sec. 

EXIT. 

(Continued  on  the  following  page) 
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20040-20090 


20100 


20170-20260 
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Notes 


TOTAL  WEIGHT  OF  FRAGMENTS  IN  EACH  POLAR  ZONE 
Sets  all  W7(J2,M5)=0.  The  variable  W7(J2,M5)  is  the 
scaled  mass  (in  grams)  of  the  fragments  in  polar  zone 
J2,  which  come  from  the  mass  point  M5.  The  J2  nd 
polar  zone  lies  between  polar  angles  5(J2-l)<a<5J2. 

Thus  polar  zone  number  36  lies  between  175£ot<l80, 
and  polar  zone  number  37  contains  only  the  180° 
direction.  Where  it  is  convenient  the  material  in 
polar  zone  37  is  included  in  polar  zone  36. 

The  scaled  mass  of  metal  associated  with  each  mass 
point,  W6„r  and  the  angles  A'*,,,.,,,  (see  Eq.  (73))  are 
found  in  20100-20500.  Fart  of  the  computation  of  the 
angles  is  done  in  subroutine  21390,  which  is  called 
from  20250,  20350,  and  20450. 

The  angles  A<lnJ4JA  (see  Eq.  (73))  and  the  velocities 
of  the  mass  points  is'ws  <L1  are  calculated  here.  The 
point  Z5,  R5  is  the  current  position  of  the  midpoint 
of  the  line  segment  joining  the  mass  points  M5  and 
M5+1,  taken  from  the  fluid  dynamics  results.  The 
variables  U5  and  V5  are  the  current  velocity  components 
of  the  mass  point  M5+1.  The  variable  Q3„g  is  the 
current  velocity  of  the  mass  point  M5.  The  computation 
of  the  angle  Ac\^vl/X  is  explained  under  21390-21690. 
Note  that  arM5,1/x  is  the  polar  angle  associated  with 
the  midpoint  of  the  line  segment  joining  mass  points 
M5  and  M5+1.  The  metal  mass  associated  with  the  mass 
point  M5  is  assumed  to  be  distributed  uniformly  between 
the  angles  and 
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20270-20360  The  current  velocities  Q3_  and  the  polar  angles 

“5 

are  calculated  for  the  mass  points  on  the 
line  L=L1. 

20370-20460  The  current  velocities  Q3„f  and  the  polar  angles 

are  calculated  for  the  mass  points  corresponding 
to  K=K1,  L*1  to  LI— 1. 


20470  is  set  equal  to  zero.  This  is  used  in  calculating 

the  spatial  distribution  of  the  mass  associated  with 
the  mass  point  M5=l. 


20480-20490  The  scaled  mass  and  velocity  associated  with  the  last 
mass  point,  M5=K1+2L1  are  calculated  here. 


20500  For  the  last  mass  point,  A%a4K4»,/Jt  is  set  equal  to 

180°. 


20520-20610  Varlcjs  quantities  used  to  get  fragment  weight  distri¬ 
butions  in  the  individual  polar  zones  are  set  equal 
to  zero. 


20650-20910  Here,  the  scaled  mass  associated  with  each  mass  point 
M5  is  distributed  uniformly  over  the  polar  region 
between  the  angles  A%r_,/x  and  A^„laa/i  .  If 
is  larger  than  the  fragments  paths  cross,  but 

the  treatment  is  the  same.  In  either  case  C7  is  the 
smaller  and  C8  is  the  larger  of  the  two  angles  defining 
the  polar  region  over  which  W4,mir  the  mass  from  the 
point  M5  is  to  be  distributed.  Also  J2  and  J3  are  the 
numbers  of  the  polar  zones  containing  C7  and  C8 
respectively.  If  C7  or  C8  is  on  the  boundary  between 
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zones  J2  and  J3  it  is  assigned  to  J3.  If  J2  is  equal 
to  J3,  all  of  the  mass  from  the  point  M5  is  assigned 
to  this  zone,  i.e.,  w7„  MJ  ,the  scaled  mass  of  metal 
from  mass  point  M5  in  polar  zone  J2,is  set  equal  to 
If  J2  is  not  equal  to  J3,  the  appropriate 
scaled  masses  are  assigned  to  zones  J2  and  J3  in 
20770  and  20730.  Appropriate  scaled  masses  are  assigned 
to  intervening  zones  in  20790-20810,  with  J4  from  20750. 
The  numbers  N7Jfc  3  of  fragments  of  weight  greater  than 
n,  and  the  scaled  masses  wsJt  J  of  fragments  in  each 
of  the  J9  weight  ranges  in  each  of  the  polar  zones 
are  accumulated  here.  Also  the  weights  of  fragments 
(in  grains)  from  each  of  the  mass  points  in  each  of 
the  polar  zones  are  printed  out. 

209*70-20970  The  total  weight  of  metal  is  summed  here  and  printed 
out, in  pounds  and  in  grains. 

210.10-21140  For  each  polar  zone,  the  total  fragment  weight  in  the 
zone  (in  grains),  the  percent  of  the  total  fragment 
weight  in  the  zone  and  the  average  velocity  of  fragments 
in  the  zone  (in  ft/sec)  are  calculated  and  printed  out. 
Here,  W8„  is  the  scaled  mass  of  fragments  in  the  polar 
zone  and  the  average  velocity, in  cra/jusec,is  w<?„ /W8„  . 
where  2«+kj 

w<*n  =]Zlw7«,M,-Q3Mt  .  The  average  velocity  in 
ft/sec  is  YJa  (see  21120). 

21150  At  this  point  the  following  subroutines  can  be  called 

(for  details  see  the  notes  under  the  subroutine 
statement  numbers).  Subroutine  21710  must  be  called 
before  either  22270  or  22480,  since  quantities  needed 
in  the  latter  two  subroutines  are  calculated  there. 
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21710-Caleulates  and  prints  out,  for  each  polar  zone, 
the  average  weight  of  fragments  in  each  of  the  J9 
weight  ranges  (in  grains)  and  the  number  of  fragments 
in  each  of  these  weight  ranges. 

22270-Values  of  the  average  fragment  velocity,  fragments 
per  steradian,  and  average  fragment  weight  (in  grains) 
are  printed  out  for  each  polar  zone.  This  form  is 
useful  as  input  for  a  particular  lethal  area  program 
(AMSAA).  The  calculated  values  can  also  be  put  on 
tape,  from  which  cards,  in  a  prescribed  FORTRAN  format, 
can  be  punched  directly. 

22480-Stored  values  are  put  on  tape  for  use  by  a 
FORTRAN  routine  (Table  AVI)  which  has  the  computer 
punch  a  card  deck  suitable  for  input  to  the  JMEM  lethal 
area  program. 

21190  End  of  fragment  prediction  subroutine  which  starts  at 

20020  and  is  called  at  6600,  in  the  fluid  dynamics 
program. 

WEIGHT  DISTRIBUTION  INPUT 

21200-21230  The  left  endpoints  t“%  of  the  J9  fragment  weight  intervals 
(in  grains)  are  read  in.  The  data  are  in  22710. 

21240-21380  In  this  case  the  Mott  distribution  (Eq.  (8l))  is  used. 

For  each  mass  point  M5,  the  quantity  RMS  }  the  average 
weight  in  grains  of  fragments  weighing  more  than  one 
grain,  is  read  in  (the  data  are  in  22720-22740).  Then 
the  corresponding  y.  is  calculated  with  Eq.  (83)  ar-l 

and  W2MS  y  are  calculated  with  Eqs.  (8l)  and  (84), 
respectively.  Notice  that  m,  and  are  in  grains, 

is  in  fragments  per  gram,  and  W2M< ,  is  dimensionless. 
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21390-21690 


21720-21750 


21770-21930 


Table  All  -  continued 


Not  os 


POLAR  ANGLE  SUBROUTINE 

The  polar  angle  A*„%„/x  (Eq.  (73))  is  calculated  here. 
The  subroutine  is  called  in  20250,  20350,  and  20*150. 

To  solve  Eqs.(72)  one  needs  tan  j3MS+1/l .  This  is  found 
as  follows:  Note  that 

tan  %H  =  tan  {(£„,+£„,**)/ 2} 

=  [l  -  Cos(pMf  +  /s  1  n  ( 

+  ) 

_  1  ~  (c Q-S1  Cos pMttl  --sin  $ t n 

•S  l rx  pMS  cos  -  cos  6-irt  )3W4t 

where  tan  AV  and  tan  pM,+1  =  vs/uS  ,  with 

U5  and  V5  (defined  in  20170-20*460)  the  velocity 
components  of  the  mass  point  M5+1.  Then  tan  1/A4, 

where  A*4=(B1+B2)/B5,  and  Bl,  B2,  and  B5  are  defined 
in  21*400-214*40.  Equations  (72),  with  the  second  part 
now  written  A4(R t-«’5)-(24-E5)=0  is  solved 
simultaneously  in  21470-21520.  Special  cases  are 
treated  at  the  end  of  the  subroutine.  It  was  assumed 
here  that  the  arena  center  (ZQ  l?0)coincides  with  the 
center  of  the  weapon.  Hence,  if  the  radial  velocity 
component  is  zero,  the  angle  is  taken  to  be 

0°  for  Z5*Z0.and  130°  for  Z5  >2 o.  The  scaled  metal 
masses  associated  with  the  mass  points  are  assigned 

in  21670. 

NUMBER  AND  AVERAGE  WEIGHT  OF  FRAGMENTS  IN  EACH 
WEIGHT  RANGE  IN  POLAR  ZONES 
Assigned  values  of  n  for  each  mass  point  are  printed 
out.  These  are  read  in  at  21250  (the  data  are  in 
22720-22740). 

The  quantities  M3Jb;r  ,  ,  H87<> ,  Y%t  are  calculated 
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Statement 

Number 


21950-21990 


22000-22260 


22270-22470 


22480-22650 


Table  All  -  continued 
Notes 


here.  Recall  that  N7,fc;  is  the  cumulative  number  of 
fragments  of  weight  greater  than  M,  in  the  J6th 
polar  zone.  The  actual  number  of  fragments  in  the 
Jth  weight  range  in  the  polar  zone,  N3yfcjJ. ,  is  found 
by  difference.  Similarly,  XI,  the  scaled  weight  of 
fragments  in  the  weight  range  in  the  polar  zone,  is 
gotten  from  the  W5JC  3  .  The  average  weight  of  fragments 
in  the  Jth  weight  range,  in  the  J6th  polar  zone, 
is  then  the  weight  divided  by  the  number  (see  21890). 

In  21900,  N8,4,  the  total  number  of  fragments  in  the 
J6th  polar  zone  is  accumulated.  The  number  of  fragments 
per  steradian  In  the  polar  zone,  Y<*JfcJ  is  calculated  in 
21920. 

The  fraction  of  the  number  of  fragments  in  the  J6th  polar 
zone,  in  the  Jth  weight  group,  ,  is  calculated  here, 

for  all  J  and  J6. 

The  array  N3Jti,  which  consists  of  the  calculated 
number  of  fragments  in  each  weight  group  for  each 
polar  zone,  is  printed  out. 

AMSAA  LETHAL  AREA  PROGRAM  INPUT 
Quantities  needed  for  the  AMSAA  lethal  area  program 
are  printed  out  here.  These  are  (see  21890), 

O'.,,,  (see  21970),  and  for  each  polar  zone  the  average 
fragment  velocity  (there  is  no  correction  for  air  drag), 
the  fragments  per  steradian  (see  21920),  and  the 

average  fragment  weight  (in  grains). 

JMEM  LETHAL  AREA  PROGRAM  INPUT 
Quantities  needed  for  the  JMEM  lethal  area  program  are 
stored  on  tape  for  use  by  a  FORTRAN  program  (Table  AVI) 
which  prepares  a  card  deck  in  the  JMEM  format.  The 
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Statement 

Number 


NOLTR  7*1-77 
Table  All  -  continued 


Notes 


quantity  YJfc  (see  21120)  is  the  average  fragment 
velocity  in  the  polar  zone  (ft/sec).  Yl^  is  the 
average  fragment  velocity  at  the  intersection  of  the 
polar  zones  J6  and  J6+1.  Also,  Mfe,w  and  N3,tiJ 
(see  21770-21930)  are  used  here. 

DATA 

22670-22700  Casing  dimensions  (called  in  1630). 


22710 


Beginnings  of  fragment  weight  range  intervals 
(called  in  21220). 


22720-227*10 


Average  weight  of  fragments  weighing  more  than  one 
grain, ,  for  each  mass  point  (called  in  21250). 


8000-8080 


K.L+l  \ 

f2 - BT" 


U _ 4j 

M- _ 

8180-8260 


K-1,  L+l  K,L+1 

r - •*, 

I  J 


( _ _ U 

K-1,  L  K,  L 

8360-8440 


K+1,L-1 


MSg^sssBga&tss 


>£$3 ?M*fr*#'*s>J>,'«»'S  jj,  «^£,  ,. 
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Table  A  III 

List  of  Variables  -  Fragment  Prediction  Code 


Variable  and 
Dimensions 

Finite 

Difference 

Notation* 

R(Kl+2 ,Ll+2 ) 

H 

Z(Kl+2,Ll+2) 

2k. 

R1 (Kl-1 , LI ) 

n 

Ri.. 

Z1(K1-1,L1) 

Z1M«. 

W(K1,L1+1 ) 

WK. 

W1 ( Kl-1 , LI ) 

wk- 

W3(K1,L1+1) 

W3„iU 

U(K1,L1+1) 

u"+t/t 

WM,l 

V(K1,L1+1 ) 

VMtW 

VI (Kl-1, LI) 

H 

Vl„. 

V5(K1-1,L1) 

v-M,‘ 

*  H,4. 

Q1(K1-1,L1) 

h-1/x 

T3(K1-1,L1) 

T3”t 

C2 ( Kl-1 , LI ) 

Description 


Radial  coordinate  of  grid  point  or 
metal  mass  point  (cm). 


Axial  coordinate  of  grid  point  or 
metal  mass  point  (cm). 


Radial  coordinate  of  cell  center  (cm). 

Axial  coordinate  of  cell  center  (cm). 

Scaled  mass  (mass/2n)  associated  with 
interior  grid  point  or  metal  mass  point 
( grams ) . 

Scaled  mass  (mass/^tr )  associated  with 
computation  cell  (grams). 

Scaled  mass  (mass/Tn )  of  gas  associated 
with  boundary  grid  point  (grams). 

Axial  velocity  component  of  interior 
grid  point  or  metal  mass  point  (cm/'/isec), 

Radial  velocity  component  of  interior 
grid  point  or  metal  mass  point  (cm/^usec). 

Relative  specific  volume  (cc/cc). 

Old  relative  specific  volume  (cc/cc). 

Artificial  viscosity  (mbar). 

Cell  time  step  fyisec). 

Sound  speed  squared  (cm/^isec)  . 


At  end  of  computation  cycle. 
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Table  AIII 

-  continued 

Variable  and 
Dimensions 

Finite 

Difference 

Notation* 

Description 

P2 (Kl-1 ,L1 ) 

P2h„ 

Pressure  +  artificial  viscosity  (PI  +Q1) 
(mbar) . 

E2(K1-1,L1) 

Specific  internal  energy  per  original 
cc  (mbar-cc/cc) . 

A1(K1-1,L1) 

aC 

Cell  area  (cm2). 

R2(K1-1,L1) 

Initial  density  in  cell  (grams/cc). 

P(K1-1) 

f; 

Burn  fraction  -plane  detonation  front 
perpendicular  to  axis. 

R3(K1) 

R3Hk 

Radial  coordinate  of  slide  point  (cm). 

Z3(K1) 

23* 

Axial  coordinate  of  slide  point  (cm). 

P3(K1) 

P3k" 

Pressure  +  artificial  viscosity  at  slide 
point  (mbar). 

R4(K1) 

Intermediate  radial  coordinate  of  slide 

point  (cm). 

Z4(K1) 

Intermediate  axial  coordinate  of  slide 

point  (cm). 

K9(K1+1) 

H 

K  coordinate  of  first  mass  point  to  the 
left  of  the  Kth  slide  point. 

S8(K1) 

S6k 

Negative  reciprocal  of  slope  of  K  line 
to  slide  point. 

U7(K1) 

Axial  component  of  velocity  of  slide 
point  (cm/^tsec). 

V7(K1) 

Radial  component  of  velocity  of  slide 
point  (cm/^sec). 

N 

N 

Cycle  number. 
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Variable  and 
Dimensions 


T 

T1 

T2 

Z 


P 


N3(37,J9) 


N7(37 , J9> 


W2(2L1+K1,JS) 

W5(37,J9) 


Mb(37  »J9) 


NOLTR  7^-77 

Table  AIII  -  continued 


Finite 

Difference 

Notation 


vi 

T2 

m 

H 


H  *itx 


<♦♦1,0 


P 

N^Jfc.T 


NT, 


T*,y 


N(2L1+K1,J9)  N, 


*9>9 


W7(37,2L1+K1)  wt 


J4,n> 


W2 


M*,J 


W5 


Jfc.T 


M4> 


y*,j 


Description 


Time. (ysec ) 

Time  step,  (ysec) 
Time  step,  (ysec) 


Axial  coordinate  of  initial  position  of 
detonation  front  for  plane  detonation. 


H-  3/x. 


Temporary  storage  for  PZ_  “Ql«, 

Number  of  fragments  in  the  Juh  weight 
range  in  the  J6th  polar  zone  (0°!£rC<5o 
for  J6=l,  .  .  .  .,175°^A<180°  for  J6=36, 
=180°  for  J6=37). 


Cumulative  number  of  fragments  of  weight 
greater  than  M,  grains  in  the  J6th  polar 
zone. 


Cumulative  number  of  fragments  per  gram 
of  weight  greater  than  M,  from  the  metal 
mar.s  point  M5. 


Scaled  mass  (mass/2n),  in  grams,  of 
fragments  from  the  metal  mass  point  M5> 
in  the  polar  zone  J6. 


Fraction  of  mass  of  fragments  from  mass 
point  M5,  of  weight  greater  than  m,. 


Scaled  mass  (mass/2ir),  in  grams,  of 
fragments  in  the  Jth  weight  range  in  the 
polar  zone  J6. 


Average  weight  of  fragments,  in  grains, 
in  the  Jth  weight  group  in  polar  zone  J6. 
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Table  AIII  -  continued 
Finite 

Variable  and  Difference 

Dimensions _ Notation _ Description 


Q6(37,J9) 

M(J9)  m, 

M1(2L1+K1) 


V/6  (2L1+K1) 

A9(2L1+K1) 

Q3C2L1+K1)  Q3„, 

W8(36)  W8„ 

W9(36)  wq„ 

Y(36)  Y„ 

N8(37)  N8,t 

Y8(37)  Y8rt 

Y9(37)  Y<*,4 


Fraction  in  the  Jth  weight  group  of  the 
number  of  fragments  weighing  more  than 
one  grain  in  polar  zone  J6. 

Weight  (grains)  at  the  left  end  point 
of  the  Jth  fragment  weight  interval. 

Value  of  F»,  the  average  weight  of 
fragments  weighing  more  than  one  grain, 
assigned  to  metal  mass  point  MS- 

Total  scaled  mass  (mass/2ir)  of  fragments 
from  the  metal  mass  point  M5- 

Polar  angle  cf  direction  of  fragments 
from  the  midpoint  of  the  line  joining 
metal  mass  points  M5  and  M5+1  (degrees). 

Velocity  of  metal  mass  point  M5 
(cm/^isec) . 

Total  scaled  mass  (mass/2v),  in  grams, 
of  fragments  in  polar  zone  J2. 

Sum  of  velocities  times  masses  contributed 
to  polar  zone  J2. 

Average  fragment  velocity  in  polar  zone 
J2  (ft/sec). 

Total  number  of  fragments  weighing  more 
than  one  grain  in  polar  zone  J6. 

Number  of  steradians  in  the  J6th  polar 
zone. 

Number  of  fragments,  weighing  more  than 
one  grain,  per  steradian  in  the  J6th 
polar  zone. 
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1110 

J9 

Number  of  fragment  weight  ranges. 

1120 

T7 

The  fluid  dynamics  print  routine  and  fragment 
prediction  routine  are  entered  every  T7 

microseconds. 

1130 

T4 

Maximum  time. 

1140 

T1 

Initial  time  step. 

1150 

E8 

Factor  for  maximum  time  step  increase. 

1160 

N4 

The  fragment  prediction  routine  is  bypassed 
for  cycle  numbers  less  than  N4. 

1170 

N5 

Number  of  times  fluid  dynamics  print  routine 
and  fragment  prediction  routines  are  entered. 

1180 

N6 

Maximum  cycle  number. 

1190 

N7 

Counter  for  N5. 

1200 

K4 

The  inert  fuze  material  ends  and  the  HE  starts 

at  K»K4+1. 

1210 

R2 

Fuze  material  density  (grams/cc). 

1220 

R4 

Metal  casing  density  (grams/cc). 

1230 

C3 

Constant  in  artificial  viscosity  formula. 

1240 

E7 

Cut  off. 

1250 

D1 

Arena  radius,  in  ft. 

1270 

R0 

Radial  coordinate  of  arena  center  (cm). 

1280 

ZO 

Axial  coordinate  of  arena  center  (cm). 

1290 

LI 

Maximum  L  in  computation  grid. 
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Table  AIV  -  continued 

Statement 

Number _ Symbol _ Description _ 


1300 

K1 

Maximum  K  in  computation  grid. 

8760 

D 

Detonation  velocity. 

8770 

D2 

Equation  of  state  constant. 

8780 

R3 

Density  of  undetonated  HE  (j*.). 

8790 

D3 

Equation  of  state  constant. 

8800 

D4 

Equation  of  state  constant. 

8810 

El 

Energy  released  jy  HE  (mbar-cc/original  cc). 

8820 

PI 

Equation  of  state  constant. 

8830 

F2 

Equation  of  state  constant. 

88J*0 

Q3 

l-viet  (where  vic,  is  the  relative  specific 
volume  at  the  Chapman  Jouguet  state. 

1520 

Z(K1,0) 

Initial  axis  position  of  inside  of  metal  at 
the  base  (cm). 

1530 

Z(0,0) 

Initial  axis  position  of  outside  of  metal  at 
the  nose  (cm). 

151*0 

Z(1,0) 

Initial  axis  position  of  inside  of  metal  at 
the  nose  (cm). 

1550 

Z(K1+1, 0) 

Initial  axis  position  of  outside  of  metal  at 
the  base  (cm). 

1610 

Z(K,0) 

Input  special  (K,0)  values  (cm). 

1620 

R(K,L1), 

R(K,L1+1) 

Ir.^ut  R(K,L1),  R(K,LI+1)  coordinates  of  metal 
casing  (In  cm).  This  can  be  done  with  sub¬ 
routines  9910  or  9960.  The  data  is  entered 
in  22670-22700. 

1650 

R(K,L1), 

R(K,L1+1) 

Input  special  R(K,L1)  amd  R(K,L1+1)  values,  if 
necessary. 

1650 
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Statement 

Number 


1750 


3560 


21210 


21220 


Table  AIV  -  continued 

Symbol _ Description _ 

Z(KI,L)  Input  special  (K1,L)  values,  if  necessary; 

for  example,  when  the  inside  of  the  base  is 
curved. 

W(K,L)  Special  coding  for  base  end  metal  masses. 

Subroutine  8910  is  useful  when  the  inside  of 
the  base  is  straight.  It  makes  the  change 
shown  in  Fig. A3. 

Subroutine  9110  can  be  used,  with  appropriate 
values  of  the  coordinates,  when  the  inside  of 
the  base  is  curved.  Here,  the  lines  separating 
the  areas  assigned  to  the  mass  points,  as  in 
Fig. 3. 

M(J)  Read  in  left  endpoints  of  the  J9  weight 

Intervals  in  grains.  The  data  are  inserted 
in  22710. 

M1(M5)  Read  in  the  assigned  value  of  FI,  the  average 

weight  of  fragments  weighing  more  than  one 
grain,  for  each  metal  mass  point  M5.  The  data 
are  inserted  in  22720-22740.  At  the  outset, 
if  these  data  are  not  available,  the  beginning 
of  the  program  can  be  run  with  the  statement 
3855  -  GO  TO  30000.  This  will  stop  the  program 
after  the  values  of  X  (Eq.(80))  have  been 
calculated.  These  values  can  be  used  to  get 
the  values  of  from  Fig.  10  or  a 

similar  plot.  After  these  values  of  are 

inserted  in  22720-22740,  remove  the  EXIT  card 
3855  and  start  over.  If  desired  some  of  the 
initial  printout  and  the  X  computation  can  now 
be  bypassed  with  1815  -  GO  TO  1880,  3620  - 
GO  TO  3860. 
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Statement 

Number 


Table  AV 

Output  Statements  -  Fragment  Prediction  Code 
_ Description _ 


1310 


1320-1370 


1820-1870 

3630-3680 


3700-3750 


3800 


43UO 

6610 


6660-7070 


HE  density  (“R3  here),  internal  energy  El(mbar-cc/ 
orig.cc)  and  equation  of  state  constants. 


Kl,  LI,  Kh,  J9,  T7,  T4,  Tl,  E8,  N5,  N6,  E7,  R2,  R4, 
R0,  Z0,  D1 

(3ee  Table  A1V  for  descriptions.) 


Case  dimension  input  -  (bypass  with  1815  GO  TO  1880) 


Initial  coordinates  of  zone  centers  (Z1Ka>R1„a)  and 
scaled  masses  (masses/2ir)  wi,^  associated  with  zones 
(bypass  with  3625  -  GO  TO  3690). 


Initial  coordinates  of  grid  points  and  scaled  masses 
(masses/27r)  wR(t  assigned  to  grid  points  and  metal  mass 
points  (bypass  with  3695  -  GO  TO  3760). 


Initial  axial  position  , H*0  , initial  inside  radius  (cm), 
C/M,  and  X  (  in*^)  (see  Eq.  (80)).  (Bypass  both  the 
computation  and  the  printout  with  3675  -  GO  TO  3860.) 


N,  T,  and  K3  are  printed  out  here  for  each  cycle. 


N,  T,  Tl,  and  T2  are  printed  out  whenever  the  fluid 
dynamics  print  routine  is  entered. 


Fluid  dynamics  print  routine..  Prints  F*  t  metal  velocity 
(ft/sec),  ZB<t  (cm)  ,  RKt  (cm)  ,  U.^cm/ysec)  (cm/ysec) 

P2„,(=  Pl  +  Ql)  (mbar)  ,  Qi„  (mbar)  ,  VL  t(cc/cc) 

E2*>L  (mbar-cc/orig.cc)  , T3KL  (ysec)  ,  C2„ik  ^(cm/ysec)  j  , 

R3k  (cm)  ,  H3„  (cm)  ,  P3R  (mbar)  ,  KS*  ,  R4„  (cm)  ,  H4K  (cm)  ,  S8„. 
Bypass  with  66^0  -  GO  TO  7080. 
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.Statement 

Number 


8670-87^0 

8850-8880 

9750-9810 

10520 

11690 

11930-11950 

20890 

20990 

21130 


mm 


NOLTR  7*1-77 
Table  AV  -  continued 
Description _ 


Diagnostic  printout  if  C2K#L<o. 

HE  and  equation  of  state  data.  Called  in  1310. 

Diagnostic  printout  if  J>Ki-lin  the  mass  point  pressure 
routine. 

Diagnostic  printout  before  exit,  if  J=K1  in  slide 
routine. 

Energy  check  printout.  Gas  internal  energy,  gas  kinetic 
energy,  metal  kinetic  energy,  and  total  energy  are 
printed  out.  These  should  be  compared  with  their 
values  at  T*0,  called  for  in  *4090. 

Metal  mass  point  velocities  (ft/sec)  on  the  ends 
(lines  K*1  and  K=Ki,  for  L*3  to  Ll-1).  These  are 
called  for  in  7020). 

Mass  point  number  M5,  polar  zone  J6,  weight  of  metal  from 
mass  point  M5  in  polar  zone  J6  (grains). 

Total  metal  weight,  in  pounds  and  grains. 

Polar  zone  J2,  weight  of  metal  in  polar  zone  J2  (grains),, 
percent  of  total  weight  of  metal  in  polar  zone  J2, 
average  velocity  in  polar  zone  J2  (ft/sec). 
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Statement 

Number 

21680 

21730-21750 

22030-22260 

22290-22470 


22480-22660 


Table  AV  -  continued 

Description _ - 

Mass  point  number  M5>  corresponding  K  and  L,  total 
scaled  mass  (mass/2ir,  in  grams)  associated  with  M5, 
polar  angle  for  fragments  from  midpoint  of 

line  segment  joining  M5  and  M5+1. 

Input  data  R(=m„t)vs  M5  (R  is  given  in  grains). 

Number  of  fragments,  in  each  weight,  range,  in  each 
polar  zone. 

AMSAA  lethal  area  input.  The  average  weighs  in  grains 
and  the  fraction  of  the  total  number  of  fragments  in 
the  polar  zone  are  printed  out  for  each  fragment  weight 
range  in  each  polar  zone.  Also  Y(J6),  the  initial 
velocity  (ft/sec)  for  the  polar  zone  (this  is  the 
average  fragment  velocity  in  the  zone  (see  statement 
no.  21120),  Y9(J6)  the  number  of  fragments  per  steradian 
in  the  polar  zone,  and  the  average  fragment  weight 
(grains)  in  the  polar  zone.  The  latter  quantity  is  not 
called  for  by  the  AMSAA  lethal  area  program. 

Calculated  quantities  needed  by  the  JMEM  lethal  area 
program  are  stored  on  tape,  for  use  by  a  FORTRAN 
program  (Table  AVI)  which  prepares  the  appropriate 
card  deck. 


table  a-vi  fortran  routine  for  preparation  of  card  deck  for  input  to  jmem  lethal 

AREA  PROGRAM  * 


PROGRAM  PCARD<  TAPE  1, PUNCH, TAP £6=PUNCH,1NPUT=65,0UTPUT  =  65) 
DIMENSION  A ( 7) 

DATA  FMT/8H(7F10.4)/ 

REWIND  1 
10  CONTINUE 

C  CALLS  A  SPECIAL  VERSION  OF  GCARD  TO  READ  FROM  TAPEl  INSTEAD  OF  INPUT 
CALL  GCARD (1.7* All) ,A<2) • A ( 3 ) » A ( 4 ) , A ( 5 ) » A ( 6 ) » A ( 7 ) ) 
IF(A(6).E0.1.E36)  CALL  GCARD { 1 ,2 . A ( 6 ) ,A ( 7 > ) 

PRINT  100, A ( 7 ) 

WR I TE (6, 100 )  A 
100  FORMAT { 7F10.4 1 
105  NDATA=A ( 7 ) 

IF< (NDATA.EQ.O ) .OR. <A(7) .EQ.l .E36) >  GO  TO  200 
CALL  BASICR ( 7*NDA TA.FMT ) 

CALL  BASICR ( 7,NDATA ,FMT ! 

GO  TO  10 

200  PRINT  201, A 

201  FORMAT ( *  ERROR  ON  INPUT  *,7F10.4> 

GO  TO  10 

END 

SUBROUTINE  BASICR ( NP.NDATA.FMT ) 

DIMENSION  A ( 19 ) 

A  ( 19 ) =1 .E36 
I NUM=0 
K=1 

10  IF ( INUM.GE.NDATA)  RETURN 

CALL  GCARD ( 1,12,A(K) ,A(K+1) ,A(K+2) ,AlK+3) ,A(K+4) ,A(K+5J ,A(K+6), 
+A(K+7) ,A(K+8),A(K+9),A(K+10) ,A(K+11) ) 

IE=!C+11 
DO  110  I=K  » I E 

IF(A( I ).EQ.1.E36)  GO  TO  120 
110  CONTINUE 
I *K+12 
120  1=1-1 

I NUM= I NUM+I 

IF ( I.LT.NP)  GO  TO  140 

JF{ ( INUM.GE.NDATA) .AND. ( I. EQ.NPJ)  GO  TO  140 
125  WR I TE( 6  ,FMT }(A(K),K=1,NP) 

IF ( I .EQ.NP )  GO  TO  150 
IS=NP+ 1 
DO  130  K=IS.I 
A ( K-NP ) =A ( K ) 

130  CONTINUE 
I=I-NP 

IF(I.GE.NP)  GO  TO  125 

C  HAVE  LESS  THAN  NP  DATA  FROM  LAST  CARO  AND  LEFTOVERS 
135  »I  +  1 

IF( INUM.EQ.NDATA)  GO  TO  145 

INUM=INUM-I 

GO  TO  10 

C  HAVE  LESS  THAN  NP  ON  ONE  CARD 
140  IF( INUM.LT.NDATA)  GO  TO  135 
145  IE=MIN0( I.NDATA+I-INUM) 

WRITEIS.FMTI ( A ( < ) ,K=1 , I E ) 

RETURN 
150  !«0 

GO  TO  135 
END 

•This  routine  was  written  by  C. 3. Wilson  (NOL/CODE332) 
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